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Preface 



This book describes in 2 parts experimental data with simple explanations (Part I) 
and itinerant electron theories (Part II) about magnetism and its related properties 
of 3J-intermetallic compounds. Unlike 3J-metal alloys and oxides, theoretical as 
well as experimental studies on SJ-intermetallic compounds such as 3J-pnictides 
and chalcogenides, on which we focus in this book, seem unfortunately delayed. 
The objective of this book is to motivate active studies in this field in the future. 

We discuss in detail magnetic and related properties of the SJ-transition-metal 
pnictides and chalcogenides, which include the intermetallic compounds expressed 
as MX and M2X, and their mixed compounds Mi-JA'^X, MXi_3;X^ and M2-xM!^X, 
where M (M') is a 3d element and X (X') a pnicogen (P, As, Sb, and Bi) or a 
chalcogen (S, Se, and Te). Most of the MX-type compounds crystallize either in 
the hexagonal NiAs-type structure or in the orthorhombic MnP-type structure which 
is regarded as a distorted NiAs-type structure. Crystallographic phase transition bet- 
ween the NiAs- and the MnP-types occurs in some of MX-type compounds when the 
temperature changes. The M2X-type compounds crystallize in the tetragonal Cu2Sb- 
type structure. As described in detail in this book, many of the compounds mentioned 
above exhibit very interesting magnetic and crystallographic phase transitions caused 
by various means such as change of temperature, applications of external magnetic 
field or pressure, and change of the composition x in the case of mixed compounds. 
This variety of phase transitions seems to have motivated many researchers to in- 
vestigate the group of compounds mentioned above. We, one experimentalist and 
three theorists, had been collaborating for a long time to clarify the mechanism of 
phase transitions as well as the peculiar magnetic properties appearing in the group of 
compounds mentioned above. Considerable progress in understanding the properties 
of the compounds mentioned above has been achieved; however, there still remain 
various phenomena not clarified yet. Besides the fundamental problems, those of po- 
tential applications are also considered for the phenomena observed in some of the 
compounds mentioned above. 

In Part I the crystallographic and magnetic properties are described together with 
simple phenomenological analyses the consideration of applications is given in a few 
sections. In Chap. 1 we show the method to prepare a poly crystalline sample, crystal 
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structures and the phase diagram of a typical MX-type compound. In Chap. 2 we sur- 
vey basic magnetic data such as magnetic transition temperatures, magnetic moments 
and magnetic orders, etc. for the compounds of MX- and M2X-types. Chapters, the 
main part in Part I, presents detailed descriptions of magnetic and related properties 
for respective compounds. In Sects. 3.1-3.6 properties of the compounds containing 
Mn or Cr, which are the most interesting from both the fundamental and applica- 
tion view points, are described in detail, together with phenomenological theory to 
explain magnetic properties, and application of MnAs as a magnetic refrigerant is 
also included in these sections. Properties of the M2X-type compounds are shown in 
Sect. 3.7. Many pages of this section are used for the Cr-modified Mn2Sb. 

Part II, where theoretical results are shown, is divided further in to two chap- 
ters. In Chap. 4 theoretical results of the NiAs-(MnP)-type compounds are shown. 
The magnetism of materials arises from 3d or 4/ electrons. There are two differ- 
ent models to describe the magnetic electrons: the localized model and the itinerant 
electron model. Theories shown in Part II are based on the itinerant electron pic- 
ture. In Sects. 4.1 and 4.2 we show and discuss in detail electronic band structures 
calculated by the first principle method for paramagnetic, ferromagnetic or antiferro- 
magnetic states of the MX-type compounds in the 3(i-pnictides and chalcogenides. 
In Sects. 4.3-4.7 we theoretically investigate various experimental data on the ba- 
sis of the calculated electronic band structures. In Sect. 4.3 we clarify origin of the 
anomalous behavior of paramagnetic susceptibility and thermal expansion of MnAs 
by taking the spin fluctuation into account. Correspondence between Fermi surfaces 
and de Haas-van Alphen oscillations in NiAs and NiSb is discussed in Sect. 4.4. 
Pressure induced magnetic transition in CrTe is discussed in Sect. 4.5. In Sect. 4.6 
we present a theory to search for the most stable magnetic ordering in FeAs, where 
instability of the paramagnetic state is investigated from the wavevector-dependent 
susceptibility tensor. In the final Sect. 4.7 structural transformation from the NiAs- 
type to the MnP-type structures is discussed. On the basis of the calculated band 
structures, we investigate the instability of the NiAs-type structure against MnP-type 
lattice distortion by taking the electron-lattice interaction into account. 

In Chap. 5 we show the theoretical results of the Cu2Sb-type compounds. Crystal 
structure and magnetic ordering are shown in Sect. 5.1. In Sect. 5.2 band structures 
of the Cu2Sb-type compounds calculated by the first principle method are shown, 
followed by a comparison between experimental and theoretical results in magnetic 
ordering and optical spectra. 

This book is the result of a longlasting close collaboration between experimen- 
talists and theorists, which is also the tradition of the study of magnetism in Japan. 
It would be our great pleasure if this book proves helpful to graduate students, 
researchers and engineers in the field of science and technology of magnetism. 



Miyagi, Osaka 
October 2009 



Kazuko Motizuki 
Hideaki Ido 
Tadaei Itoh 
Masato Morifuji 
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Experimental 



1 



Basic Properties of 3J-Pnictides 



1.1 Main Compounds of 3rf-Pnictide 

Among various intermetallic compounds between 3d transition metal M and 
pnicogen X (X = P, As, Sb, or Bi), compounds expressed as MX and M2X are 
of special interest due to their rich variety of magnetic properties. The MX com- 
pounds crystallize basically in the hexagonal NiAs-type (B81) structure, however, 
some of them show distortion into the orthorhombic MnP-type (B31) structure. 
The M2X compounds take the tetragonal Cu2Sb-type structure. The MX and M2X 
compounds are listed in Table 1.1. 

In this table, the MX compounds denoted by boldface letters take the NiAs- or 
MnP-type structure. The M2X compounds written in boldface take the Cu2Sb-type 
structure. It should be pointed out that the ternary compounds MnAlGe and MnGaGe 
also have the Cu2Sb-type structure. Some compounds take nonstoichiometric com- 
position such as Mi+;cX {x > 0) although they are written as MX in Table 1.1 [1]. For 
example, MnSb actually forms a compound Mni+^^^Sb with 0 < x < 0.2. Some other 
NiAs-type compounds also have similar tendency [1]. As the excessive Mn atoms, 
which occupy interstitial positions, have a strong effect on magnetic properties, a 
special attention should be paid for preparation of stoichiometric MnSb compound 
[2] . For the NiAs-type structure, we can regard that M atoms with smaller size than 
that of X atoms occupy interstitial positions of the hexagonal close-packed struc- 
ture formed first by the X atoms. Therefore, this may permit existence of excessive 
M atoms. Among the compounds shown in Table 1.1, only Cr-, Mn-, and Fe- com- 
pounds take magnetic ordering. 



1.2 Crystal Growth by Sintering Method 

Crystal growth is the starting point of experimental studies of solid state physics. 
Although there has been rapid progress in equipment and technique, we show here 
a conventional method to prepare a polycrystalline sample. This method is called 
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Table 1.1. Compounds between 3 J metal and pnicogen 





P 


As 


Sb 


Bi 


Ti 


TiP 


TiAs 


TiSb 


- 


V 


VP 


VAs 


VSb 




Cr 


CrP 


CrAs, Cr2As 


CrSb 




Mn 


MnP, Mn2P 


MnAs, Mn2As 


MnSb, Mn2Sb 


MnBi 


Fe 


FeP, Fe2P 


FeAs, Fe2As 


FeSb 




Co 


CoP, C02P 


CoAs 


CoSb 




Ni 


NiP 


NiAs 


NiSb 


NiBi 


Cu 






Cu2Sb 





The compounds written in boldface have the NiAs-, MnP-, or 
Cu2Sb-type structure 



a sintering method, in which mixture of powders of constituent elements settled in 
an evacuated silica tube is sintered in a furnace. This method is applicable for most 
intermetallic compounds. 

Usually particle size of metal powder is about 100 |im, but larger or smaller par- 
ticle sizes can be acceptable. Mn can be ground easily in a mortar of agate; however, 
other 3d metals are very difficult to pulverize. It is recommended to obtain the pow- 
dered metal from specialty store. Pnicogen is ground easily, but special care is nec- 
essary for phosphorus and arsenic. As for phosphorus, we use red phosphorus (not 
white phosphorus), which may explode while grinding. One can purchase powder of 
red phosphorus. Since As is toxic (AS2O3 is deadly poisonous), one must wear mask 
not to inhale the fine particles and be careful so that particles do not scatter. Tidy- 
ing up after experiments is also necessary. Purity of all ingredient should be more 
than 99.9%. 

We show procedure of sintering for the case of MnAs. First, we wash electrolytic 
manganese in dilute nitric acid and then in water. Next we grind it in a mortar of 
agate. High purity bulk of As taken out of a vacuum glass container is crushed into 
powder. Powders of Mn and As are weighted in desired proportion and mixed in 
a mortar carefully so that the powder does not scatter. These processes should be 
finished as quickly as possible to avoid oxidation of the powdered metals. Then, the 
mixture of Mn and As (total weight of 5-8 g) is sealed in an evacuated quartz tube. 
An example of shape of the quartz tube is shown in Fig. 1.1. 

The neck of the quartz tube should be carefully prepared. If this part is too thin, 
keeping the tube in high vacuum may fail, and ingredient powder cannot go inside 
through the neck with too small internal diameter. 

The mixture sealed in the evacuated quartz tube is heated in an electric furnace. 
Since As sublimates at 613°C and has large vapor pressure at high temperature, fur- 
nace temperature must be raised slowly spending a few days up to 500°C in order to 
avoid explosion. Then the temperature is kept constant for another few days so that 
rough reaction of the mixture finishes. In the final stage, the temperature is raised 
up to about 800°C (you can heat up rapidly this time) and kept for a week. We can 
take out the reacted product after cooling down. As the sintered crystal thus prepared 
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ingredient 



evacuation 



c 



outer 0 4 ~ 5mm 
inner 0 2-^ 3mm 



Fig. 1.1. A sample tube of quartz. After evacuating up to about 10~^ Torr, the neck of the 
quartz tube is fused and cut by a gas burner. The dimensions shown in this figure are typi- 
cal ones 



may be inhomogeneous, the sintering at 800°C for a week is desirable to be repeated. 
Namely, the reaction product is pulverized and sealed again in an evacuated quartz 
tube. One can heat the sample rapidly this time. Longer the period of sintering, better 
the quality of the crystal. 

As phosphorus also has risk of explosion, the procedure mentioned earlier must 
be taken for compounds containing P. Other compounds also can be prepared in the 
similar way. The diagram of binary alloys [3] will be helpful for determining the 
sintering temperature. 



1.3 Crystal Structure of NiAs-, MnP-, and Cu2Sb-Type Structure 

1.3.1 NiAs-Type Structure 

As the MnP-type structure is regarded as distorted NiAs-type structure with small 
shifts of atomic positions, these two structures are similar to each other. We explain 
below these crystal structures connecting them with each other. 

Figure 1.2 shows the hexagonal NiAs-type structure. Pnicogen atoms form a lat- 
tice like the hexagonal close-packed structure, and 3d metal atoms a simple hexag- 
onal lattice. Table 1.2 shows lattice constants of the MnX compounds at room 
temperature. 

As seen in Table 1.2, with change of pnicogen from Bi to P, the lattice constants 
decrease. In particular, decrease of the a-axis is remarkable. It is noted that a-axis of 
MnP is 73% of that of MnBi, corresponding to the fact that radius of P atom is 75% 
of that of Bi. Hence we see that the a-axis is almost proportional to radius of the 
pnicogen atom. Such tendency is also found in the compounds CrP, CrAs, and CrSb. 

1.3.2 MnP-Type Structure 

The MnP-type structure is obtained by introducing an orthorhombic distortion as 
well as slight shifts of atom positions for the hexagonal NiAs-type structure. As 
distortion and shifts are small, these two structures are basically similar. However, 
even the slight difference affects magnetic properties through change of electronic 
structure. 



6 1 Basic Properties of 3J-Pnictides 




Table 1.2. Lattice constants a and c (room temperature) of MnX [4] and radius of the X atom 





(MnP) 


MnAs 


MnSb 


MnBi 


a (A) 


3.17 


3.68 


4.15 


4.34 


c(A) 


5.26 


5.72 


5.78 


5.98 


Atom radius of X (A) 


1.28 


1.39 


1.59 


1.70 



Note that MnP takes the orthorhombic MnP-type structure 
(a and c are lattice constants corresponding to those of the 
hexagonal NiAs-type structure) 




M: (0,0,0); (0,0, 1/2) 
0,0 X : (2/3, 1/3, 1/4), (1/3, 2/3, 3/4) 



Fig. 1.3. Relation between the NiAs-type and the MnP-type structure projected onto the 
c-plane. (a) The NiAs-type structure (The thick lines denote the unit cell. See also Fig. 1.2.) 
(b) Distortion into the orthorhombic lattice and displacement of M atoms along the Z^-axis are 
shown. The displacement is written as ub using a displacement parameter u. (c) Displacement 
of X atoms in the /^-direction. The displacement is written as wb (w: 3. displacement param- 
eter) whose origin is 1/3 or 1/6 of the Z^-axis. Note that the displacements and distortion are 
enlarged in order for the illustration to be legible. See also Fig. 4.38 



We explain the difference between the MnP-type and the NiAs-type structure. 
Figure 1.3 shows atomic positions projected onto the c-plane for the NiAs-type and 
MnP-type structures, where displacements of atoms in the c-plane are also illustrated. 
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In Fig. 1.3, crystallographic parameters u and w indicate atomic displacement 
along the /7-axis. The parameter w is one order smaller than u. Distortion from the 
NiAs-type to MnP-type structure is expressed with a parameter 5, which denotes 



Atoms also shift along the c-axis as shown in Fig. 1.4. 

In summary, the MnP-type structure is characterized by the five crystallographic 
parameters: the displacement parameters of metal (u and x) and anion (v and w) and 
the distortion parameter (5). The NiAs-type structure is thus expressed by setting 
u = v = x = w = 5 = 0. The crystallographic parameters are listed in Table 1.3 for 
MnP and CrAs, which take the MnP-type structure. 

Other MnP-type compounds have similar values of the parameters. Since w 
and X are generally smaller than u and v, only v, and 5 can be determined in 



Fig. 1.4. Unit cell of the MnP-type structure. Displacements of M atoms and X atoms in the 
c (A) -direction are given by xc and vc, respectively. The parameter x is smaller than v by one 
order. In the expression of lattice parameters of the MnP-type structure, (A, 5, C) is often used 
in addition to {a,b,c). Note that displacements are enlarged, and displacement of X atom in 
the /7-direction (w) is not shown to avoid confusion (see Fig. 1.3) 

Table 1.3. Crystallographic parameters of MnP and CrAs in the MnP-type structure (room 
temperature) [4] 



deviation of the /?-axis from >/3a-axis : 



5 = {b-V3a)/V3a. 




a(B) 



9 © : M atom 

O O ' Xatom 



5 



u 



V 



X 



w 



MnP 0.077 0.05 0.06 0.005 0.01 
CrAs 0.036 0.05 0.05 0.007 0.006 
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Table 1.4. Positions of atoms in the MnP-type structure expressed in the coordinates based on 
A-, and C-axes 

M {-x,0,u),{x,l/2,l/2-u), (l/2-x,0,-^), (1/2 + x, 1/2, 1/2 + ^) 

X (1/4 - V, 1/2, 5/6 -w), (l/4 + v,0,l/3-w), 

(3/4 - V, 1/2, 1/6 + w), (3/4 + v, 0, 2/3 + w) 




Fig. 1.5. Unit cell of the Cu2Sb-type structure. 3d transition metal occupies the Cu(I)- and 
Cu(II)-sites. In MnAlGe, Al atoms occupy the Cu(II)-sites and Ge occupies the position of Sb. 
The crystallographic parameters u and v are these of Mn2Sb 

reliable accuracy by X-ray diffraction for powdered samples. Table 1.4 shows posi- 
tion of atoms (see Figs. 1.3 and 1.4). We note that the coordinates are based on the 
A-, B-, and C-axes (c-, a-, and /7-axes) shown in Fig. 1.4. 

In some literatures, origin of coordinate for the atomic position is shifted to 
(1/4, 1/4) in the ab (5C)-plane shown in Fig. 1.3 (or 1.4). In this case, the positions 
shown in Table 1 .4 must also be shifted. There have been several (about four) man- 
ners of expression for the lattice parameters of the MnP-type structure. The (a, /?, c)- 
system used in Figs. 1.3 and 1.4 is useful because the a- and c-axes of the NiAs-type 
lattice correspond to the a- and c-axes of the MnP-type structure. We use mainly 
(a^b^c) system and sometimes (A,5,C) system in this book. 

1.3.3 Cu2Sb-Type Structure 

MnAlGe, MnGaGe, and M2X denoted by boldface letters in Table 1 . 1 take this crys- 
tal structure, whose unit cell is shown in Fig. 1.5. In this structure, distance between 
the Mn(II) atoms (denoted by the thick dotted line) is large. For example, the distance 
is 3.94 A for Mn2Sb. Nearest distance between Mn(I)-Mn(II) and Mn(I)-Mn(I) is 
2.81 and 2.88 A, respectively. In the case of MnAlGe, which is obtained by substi- 
tuting Mn(II) and Sb of Mn2Sb by Al and Ge, respectively, Mn-Mn distance along 
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the c-axis is very large and equal to c (= 5.93 A). This is also the case of MnGaGe. 
Crystallographic characteristics mentioned above affects magnetic properties, which 
we discuss in Chap. 3. 



1.4 Nonstoichiometric Compounds 

One of typical nonstoichiometric MX compounds is Mni+j^Sb. Figure 1.6 is a phase 
diagram of Mn-Sb binary system [3]. The NiAs-type structure exists at room tem- 
perature for Sb composition between 44% and 49%. In other words, Mni+j^Sb with 
0.04 < X < 0.25 has the NiAs-type structure. Curie temperature Tq of Mni+j^Sb 
decreases linearly from 314°C (51at% Mn) to 90°C (56at% Mn) with increasing 
Mn composition. This diagram shows that stoichiometric (i.e., 50at% Mn) MnSb 
does not exist. MnSb is considered to have about 1% excessive Mn atoms occu- 
pying partially the interstitial sites ((2/3, 1/3, 1/4) and (1/3, 2/3, 3/4)) shown in 
Fig. 1.7. It has been reported that an excessive Mn atom at the interstitial site is 
surrounded by six nearest Mn atoms and does not have magnetic moment, while 
an Sb atom polarizes with small moment 0.2/iB antiparallel to the total moment [5]. 
For Mni+;cSb (0.04 <x< 0.25), Tq and spontaneous magnetization rapidly decrease 
with increasing x. Magnetic anisotropy constants also vary with x [6] . 

Some other MX compounds that forms the NiAs- (or MnP-) type structure in 
nonstoichiometric composition region are shown in Table 1.5. 



1000°C 



800 



U 
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200 
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Fig. 1.6. A phase diagram of the Mn-Sb binary system. The label 2P denotes two phases 
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O Mn 

9 Excessive Mn 

(partially occupied) 

Osb 



Fig. 1.7. Interstitial positions of the NiAs-type structure at which excessive Mn atoms of 
Mni+;(;Sb locate are denoted by the filled circles 



Table 1.5. Composition region where MX-type compound exists (room temperature) 
VSb Single phase only for x = 0.4 [7] 

CrSb 0 <x< 0.04. Upper limit of x increases with temp [1,3] 
MnSb See Fig. 1.6 

FeSb 0.2<x< 0.3 [8], 0.1 < x < 0.3 [3] Region incr. with temp [8] 

Exists only for excessive Fe 
CoSb 0<x< 0.04. Region incr. with temp [1]. -0.04 < x < 0.15 [3] 
NiSb -0.06 <x< 0.02. Region incr. with temp [1]. -0.08 < x < 0.10 [3] 
VAs 1:1 compound 
VP 1:1 compound 
CrP 1:1 compound 
CrAs 1:1 compound 
MnP 1:1 compound 
MnAs 1:1 compound 

MnBi MnBi for T < 340°C. Mni.ogBi for T > 340°C 

Phase separation into Mn and Mn-Bi hquid for T > 446°C [3, 8] 

The letter x in this table means composition in Mi+^^X. Generally, the region 
expands as temperature rises 
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Overview of Magnetic Properties of NiAs-Type 
(MnP-Type) and Cu2Sb-Type Compounds 



2.1 Compounds That Have Magnetic Ordering Phase 

In this section, we look over structural and magnetic properties of the MX and M2X 
compounds. In Table 2. 1, we list the NiAs- (or MnP-) type and the Cu2Sb-type com- 
pounds that show magnetic ordering. Basic magnetic properties of typical MX and 
M2X compounds are shown in Tables 2.2-2.6. In the tables, only representative data 
are shown when there are several data. Refer to Chap. 3 for more detailed magnetic 
properties. 

As shown in Tables 2.1-2.4, magnetic ordering occurs only in compounds con- 
taining Cr, Mn or Fe. The compounds with the MnP- type structure show the double 
helical magnetic ordering as shown in Table 2.6. 

As a typical example of the double helical structure, we explain the magnetic 
structure of CrAs in Fig. 2.1. 

We summarize the double helical magnetic orderings of the MnP-type com- 
pounds in Table 2.6. MnP has a long wavelength and thus is close to a ferromag- 
net, and is known to show metamagnetic magnetization process. Magnetization of 
MnP saturates at magnetic field about 4kOe when the applied field is along the 
a-axis. There are no magnetization measurements for the other compounds listed 
in Table 2.6. There is only an unpublished data by us for powdered sample of CrAs 
measured in pulsed fields up to 35 T at liquid nitrogen temperature. According to 
this report, magnetization shows no saturation and reveals anomaly at lOT. CrAs is 
an interesting compound with various unsettled problems. We will discuss some of 
these problems in Chap. 3. 

MnX (X = As, Sb, or Bi) are ferromagnets. CrAs, CrSb, MnAs, and MnBi show 
a first order phase transition at Tq (or 7n). See Chap. 3 for the detail of these phase 
transitions. 

The Cu2Sb (M2X)-type compounds have rich variety of magnetic structures 
as shown in Fig. 2.2 [1]. In particular, ferrimagnetic Mn2Sb is interesting; Slight 
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Table 2.1. The shape of the signs indicates crystal structure at room temperature 





\m 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


p 


MP 




0 


0 


0 


0 


0 




As 


MAs 
M^As 




0 


0 
□ 


o 
□ 


0 

□ 


0 


o 


Sb 


MSb 
M^Sb 


O 


o 


o 


o 


o 


o 


o 


Bi 


MBi 








o 









Q: NiAs-type Q- MnP-type □: Cu^Sb-type 



Compounds that show magnetic ordering are denoted by the gray sign 



Table 2.2. Magnetic properties of CrX [1] 



MX 


CrP 


CrAs 


CrSb 


Magnetic structure 


X has a peak ~ 200 K 


Double helix 


Antiferro 


7n(K) 




250 (1st) 


718 (1st) 


PA (MB/Cr) 


;^(RT^3)xlO-6 


1.67 


3.0 



Unit of magnetic susceptibility x is emug~ . See also Chap. 3. and 
Pa are Neel temperature and magnetic moment per atom. (1st) means 
first order transition 



Table 2.3. Magnetic properties of MnX [1] 



MX 


MnP 


MnAs 


MnSb 


MnBi 


Magnetic structure 


Double helix (T < 47 K) 


Ferro 


Ferro 


Ferro 




Ferro (47K<r<291K) 








7c (K) 


291 


318 (1st) 


587 


628 (1st) 


PA (MB/Mn) 


1.3 


3.4 


3.6 


3.8 



See also Chap. 3. Tc is Curie temperature. (1st) means first order transition 

substitution of Cr etc. for Mn induces a first order phase transition from an anti- 
ferromagnetic to a ferrimagnetic state with increasing temperature. We will discuss 
this phenomenon in Chap. 3. The various magnetic ordering shown in Fig. 2.2 is also 
discussed in Part II of this book. 
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Table 2.4. Magnetic properties of FeX [1] 



MX 


FeP 


FeAs 


Fei+;,Sb 


Magnetic structure 


Double helix 


Double helix 


Triangular in c-plane 


7n(K) 


125 


77 


105 ~ 21 1 (depends on x) 


PA (Mb/Fc) 


0.41 


0.51 


0.88 



Table 2.5. Magnetic properties of M2X (See [1] and Chap. 3) 



M2X 


Cr2As 


Mn2As Fe2As 


Mn2Sb 


Magnetic structure 


Antiferro 


Antiferro Antiferro 


Ferri 


7n(K) 


393 


573 353 


550 


PA (Mb/M(I),M(II)) 


0.40, 1.34 


2.2,4.1 1.28,2.05 


2.13,3.87 


Table 2.6. Double helical structure 






MnP 


CrAs FeP FeAs 




On 


16° 


-126° 169° 140° 




2n/k{= 


A) 8.2/? 


2.61b 5M 2.61b 




See Fig 


2.1 for 012. 


b is length of the Z^-axis. 





Data are taken at 4 K (12 K for FeAs) (see [1] 
and Chap. 3) 



be ©6 




• • 1 X^27r/k = 2.67 b 



Fig. 2.1. Double helical magnetic ordering of CrAs. (a) Cr atoms projected onto the aZ^-plane 
(see Fig. 1.3b). The spins are parallel within the plane 1, 2, 3,. . . normal to the Z^-axis. (b) As 
the plane changes 1, 3, 5,. . . (or 2, 4, 6,. . .), the direction of the spin rotates by 67.4° with 
the relative angle 126° between adjacent (i.e., between 1 and 2) planes. The wavelength is 
A = 2.61b and /icr = 1-67/iB [2]. Slightly different data are reported in [3]. See also Fig. 4.33 
in Part II 

2.2 Compounds without Magnetic Ordering 

The compounds denoted by unshaded symbols in Table 2.1, which do not include 
Cr, Mn, or Fe, are paramagnetic or diamagnetic. We summarize properties of these 
compounds that show no magnetic ordering in Table 2.7. 



14 



2 Overview of Magnetic Properties of Ni As-Type and Cu2Sb-Type Compounds 



0«M(I) 

® M (II) 



unit 
cell 




CrjAs 
(antiferro) 



MnsAs FgoAs 
(antiferro) (antiferro) 



(>► 

Mn^Sb 
(ferri) 



Fig. 2.2. Magnetic orders of the Cu2Sb-type compounds. M(I) and M(II) (see Fig. 1.5) have 
different moments. As (Sb) is not shown [1] 



Table 2.7. Properties MX compounds that have no magnetic ordering 



MX 


Struct. 


Magnetism (unit of x is emu g ^) 




TiP 


TiP 






TiAs 


TiP 


PauU para. ;^ ^ 1.5 x IQ-^ (T = 4.2 K) 


[4] 


TiSb 


NiAs 


Pauh para 


[5] 


VP 


NiAs 


X slightly decreases with temp, x — ^x 10~^ (RT) 


[6] 


VAs 


MnP 


X slightly increases with temp. ;^ ^ 1.5 x 10~^ (RT) 


[6] 


VSb 


NiAs 


Single phase in Vi.4Sb. Pauli para, x — 1-5 x 10~^ 


[7] 


CrP 


MnP 


X has a broad valley at 200 K. - 3 x 10"^ (RT) 


[8] 


CoP 


MnP 


PauU para. ;^ ^ 1.5 x 10"^ (T > 100 K) 


[9] 


CoAs 


NiAs 


X has a broad peak at 200 K. - 2 x 10"^ (RT) 


[10] 


CoSb 


NiAs 


Pauli para. ;^ ^ 1 x 10~^ 


[4] 


NiP 


NiP 






NiAs 


NiAs 


PauU para. ;^ ^ 2.5 x 10"^ (T > 300 K) 


[11] 


NiSb 


NiAs 


Diamag 


[12] 


NiBi 


NiAs 


Superconducting, Tc = 4.25 K 


[13] 



IVIX compounds that have no magnetic ordering are of little interest in mag- 
netism. However, we can investigate genuine electronic structures that are helpful 
to understand properties of other IVIX compounds. In Table 2.7, some of such com- 
pounds with weak temperature dependence of susceptibility are denoted by "Pauli 
para". The susceptibility x of the order of 10~^ to 10~^ emug~^ may be considered 
by not only 3d states but also diamagnetism of anion, and possibly by some kind of 
magnetic impurities. Some compounds such as VSb has a nonstoichiometric phase 
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in which the composition deviates from 1:1. Magnetic properties are often affected 
by nonstoichiometry. As already described in Sect. 1.4, ferromagnetic Mni+^^^Sb is a 
typical example of nonstoichiometric compounds [7, 14, 15]. 
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Properties of the Compounds with NiAs-Type 
(MnP-Type) and Cu2Sb-Type Structures 



3.1 MnP and Related Compounds 

A magnetic phase diagram of MnP with respect to magnetic field and tempera- 
ture is shown in Fig. 3.1. This diagram was obtained from magnetization measure- 
ments for a single crystal MnP [1]. It has been known that in the temperature region 
r < 47 K, MnP shows a double helical ordering with relatively long wavelength 
of 8.2/? propagating along the Zp-axis [6]. Above 47 K, MnP is a ferromagnet with 
7c = 291 K. For 420K < T < 660 K, susceptibility obeys the Curie-Weiss law with 
jUeff = 2.36/iB {'^S = 1.56) and Op = 344 K. This value of /igff is close to saturation 
moment 1.3/iB observed at T = OK [7]. The X~^-T curve shows upward convex in 
the range 660 K <T< 1400 K [8]. Studies on mixed compounds with the type of 
Mni-JAjP (M = other 3d metal) also have been reported [9-12]. 



3.2 MnAs and Related Compounds 

3.2.1 Magnetic Transition of MnAs and the Bean-Rodbell Theory 

MnAs is especially interesting among the MX compounds. In this section we show 
magnetic and crystallographic properties of MnAs and give a simple explanation for 
them on the basis of the phenomenological theory developed by Bean and Rodbell. 
We also mention some interesting properties of mixed compounds between MnAs 
and MnP (and MnSb).^ 

MnAs has long been of interest since the study carried out by Guillaud et al. 
[13]. Temperature dependence of magnetization and inverse susceptibility are shown 
schematically in Fig. 3.2 [14]. With rising temperature, spontaneous magnetization 
disappears discontinuously at Tq (up) = 3 1 8 K. This first order transition is accom- 
panied with structural transformation from the Ni As-type to the MnP-type structure. 
Susceptibility shows a peak at 7^ = 398 K and obeys the Curie-Weiss law for T > Ti. 



See Part II of this book for the itinerant theory on magnetism of MnAs. 
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fan struct. \ 


ferro-Hke structure 






double hcHx 
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T =291.5 



0 T 100 200 300 

f (K) 

Fig. 3.1. Magnetic phase diagram of MnP [1-4]. See also [5] for more detailed diagram 
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X ' e„==260K 

Pa = 3.56//5 

6'^-288K 

T(K) 

T, ( = 398K) 
Te (ext.) 

Fig. 3.2. Spontaneous magnetization CTs and inverse susceptibility of MnAs plotted as 
functions of temperature 



As we will see later, 7^ is a structural transition temperature from the MnP-type to 
the Ni As-type structure. Anomalous temperature dependence of paramagnetic sus- 
ceptibility in the temperature range 7c < T < 7^ is explained from the viewpoint of 
itinerant electrons (see Sect. 4.3). 

Figure 3.3 shows temperature dependence of the crystallographic parameters 
[15]. The (2-axis and the unit cell volume show remarkable discontinuous changes 
of 1.2% and 2.1%, respectively, at 7c(up), while the c-axis does not show visible 
change at 7c (up). These facts suggest that there is a correlation between the a-axis 
and ferromagnetism of IMnAs. Thermal expansion coefficient has a large value in 
the IMnP-type phase just above 7c. This is due to the rapid structural transformation 
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(A) 




100 200 300 400 500 300 To 350 400 

T (K) T (K) 



Fig. 3.3. (a) Temperature dependence of lattice constants of MnAs. (a, b, c are lattice constants 
of the hexagonal lattice.) a is the thermal expansion coefficient, (b) Temperature dependence 
of the parameters u, v, and 5. Xmax means the temperature where the ;^~^-r curve shows the 
minimum as shown in Fig. 3.2 (Reproduced from [15]) 



from the MnP-type to the NiAs-type in the small range of temperature (see Part II for 
detail). The crystallographic parameters u, v, and 5 (see Table 1 .3) take the maximum 
value just above 7c (up), approaching zero at 7^ on further increasing temperature as 
indicated in Fig. 3.3b. 

MnAs shows a metamagnetic magnetization process in the temperature range 
Tc<T <Tt, corresponding to a paramagnetic-ferromagnetic transition accompanied 
with a crystallographic change from the MnP-type to the NiAs-type structures. As 
the magnetization process is similar to that of MnAsi-^Px, which we will discuss 
later, we show in this section only temperature change of magnetization in if = 2 T 
and if = 40 T in Fig. 3.4 [16]. The curve in if = 40 T suggests that, if the NiAs-type 
structure below Tq were maintained up to the temperature region just above Tq of 
the spontaneous magnetization, Tq could be about 450 K. It is noted that this roughly 
extrapolated Curie temperature Tq — 450 K is considerably higher than Op o:^ 270 K 
estimated in the high-temperature NiAs-type region, which is characteristic of MnAs 
as 7c ^ Op is satisfied for MnSb and MnP. 

Various characteristic properties of MnAs mentioned above are well explained 
by a phenomenological theory based on an assumption of strong volume dependence 
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Fig. 3.4. Temperature dependence of magnetization of MnAs in two magnetic fields H = 2 
and 40 T. Tq means 7c(up) (Reproduced from [16]) 



of exchange interaction between Mn atoms [17, 18]. Namely, the first order transition 
of IVLnAs (such as the discontinuous change of crystal volume and spontaneous mag- 
netization, and their hysteresis for temperature variation), metamagnetic phase tran- 
sition in the temperature region just above Tq (e.g., magnetization process and tem- 
perature dependence of hysteresis), and pressure dependence of Tq can be generally 
explained. We have to note, however, that 3d electrons in IVLnAs from which magne- 
tization arises have itinerant character. Therefore, it is necessary to treat 3d electrons 
in IMnAs as itinerant electrons as we will show in Part II of this book. In this section 
we describe the phenomenological theory to show how the various peculiar proper- 
ties of IVLnAs are explained. 

First, we assume that the Curie temperature depends on crystal volume (or lattice 
constants) [17]. This assumption is reasonable as Tq depends on applied pressure 
in most ferromagnets. The crystal volume is affected by various factors, for exam- 
ple, thermal expansion, pressure, applied magnetic field, and magnetovolume effect 
(magnetostriction). Structural phase transition also changes crystal volume. Consid- 
ering such factors. Bean and Rodbell [17] set a relation 

Tc{V) = To[l^p{V-Vo)/Voi (3.1) 

where Tq is assumed to be a linear function of crystal volume V. This may be justified 
as long as volume change is small. In (3.1), Vq is the volume at T = 0 K in the absence 
of magnetovolume effects. Tq is an expected Curie temperature when V is assumed 
to be kept in the constant Vq. Volume dependence of Tq is given by the coefficient j3, 
which takes a large value for IVLnAs. 

We consider a localized model where a magnetic material consists of N magnetic 
atoms with spin j. Under an applied magnetic field H and pressure P, free energy of 
this material is written in the form 
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G =-(3/2) + \)]Nk^TcG^ - HM,G + {\/2K) [(V - yo)/yo]' 

-Tk^NS^P[{y-yo)/VQ], (3.2) 

where a(= Ms(r)/Ms(0)) is a magnetization normalized by Ms(0) specific magne- 
tization at OK. The first term denotes exchange energy written by using Tq? The 
second and the third terms are Zeeman energy and elastic energy {K is the com- 
pressibility), respectively. The fourth term is the entropy one and k^NS {k^ is the 
Boltzmann constant) is the entropy, where S is summation of magnetic contribution 
So {a J) and lattice contribution S\^. The last term denotes energy due to applied 
pressure. 

When G is smaller enough than 1, ^'cr of the magnetic entropy term k^NSa{(y^ j) 
can be expanded in the following form: 

Sa{G, j) = Sa{0, j) + aa^ + + . . (3.3) 

where the coefficients a,b, ... are given by 

a=-(3/2)[7/a+l)], b=-{9m[{lj+lf-l]/[l{j+l)]\ (3.4) 

See Appendix for the detailed expression of (3.4). This expansion is useful to in- 
vestigate magnetic phase transition. We note that the exact expression of SaioJ) is 
given by 



Scj{OJ)-Sa{(Jj) = aa- 




(3.5) 



where Bj{a) is the Brillouin function with angular momentum j. For the discussion 
in this section, we use (3.3) and (3.6), which is derived from (3.5) as given by 

dSa/da = -a. (3.6) 

Substituting (3.1) into (3.2), and using a condition of minimum free energy 
dG/dV = 0,wQ have a relation 

{V -Vo)/Vo = {3/2)[f /j{j^l)]NkBKTol3a^-PK. (3.7) 

This equation means that magnetovolume effect (or exchange striction) is propor- 
tional to a^ifP = 0. 

First, we consider a case of P = 0. Inserting (3.7) into (3.2), we obtain an expres- 
sion of G, which has been minimized with respect to V. For this free energy, we use 
a condition dG/da = 0 and (3.6), and then we have a relation among spontaneous 
magnetization a(= Ms(r)/Ms(0)), temperature T, and a magnetic field H as 

{gj^B/kB){H/To) = 2aa[l + {2b/a)rio^] + {T/To)a • • • , (3.8) 



^ Exchange energy is written as — ( 1 / 2)AM^ with the coefficient A of molecular field. The 
first term in (3.2) is derived from this expression. 
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where /Ib is the Bohr magneton and g is the g factor (=2 in the present case). In this 
equation T] is defined by 

7] = -{a^/2b)NkBKTol5^. (3.9) 

We note that rf is proportional to As we will show below, rf is an important 
parameter contained in the coefficient of the O"^ term of the free energy expansion. 
For T] > 1, a first order transition occurs, while order of transition is second for T] < 1 
(see also Appendix). The coefficients a and b are given in (3.4), and a is a function 
of a (i.e., a = Bj{a)). 

Using (3.8) with a given value for j and setting // = 0, we can calculate tem- 
perature variation of spontaneous relative magnetization a = a{T /Tq) for various 
values of t] . We can also calculate magnetization curve by fixing temperature and 
changing H. By considering that magnetization of MnAs at T = 0 K is 3.4/iB, we set 
j = 3/2 (with g = 2). Figure 3.5 shows temperature dependence of a calculated from 
(3.8) for various values of rf. Magnetization curves calculated for t] = 2 at various 
temperatures are plotted in Fig. 3.6. 

We see from Fig. 3.5 that 7] = 1 is the critical value between the first and the sec- 
ond order transition. For T] = 2, relative spontaneous magnetization O" is 0.81 at ^ = 1 
and is 0.6 at ^ = 1.08 (i.e. just below 7c). These results agree well with experimental 
values shown in Fig. 3.4. We also see that magnitude of hysteresis of spontaneous 
magnetization is A{T/To) = 0.08. This value is also close to the experimental value 
A(r/7b) = 0.06 (i.e., 7c (up) = 318 K and 7c (down) = 7b = 301 K). Therefore, tem- 
perature dependence of spontaneous magnetization of MnAs is well reproduced with 
the parameter rf =2. 



1.5 




T/T, 

Fig. 3.5. Relation between magnetization (cr) and normalized temperature (T/Tq) for various 
values of r]. The curves were calculated on the basis of (3.8). The curve with r/ = 2 is in good 
agreement with experimental curve for MnAs 



3.2 MnAs and Related Compounds 23 

1 



0.8 

0.6 

O 

0.4 

0.2 



-0.01 0 0.01 0.02 0.03 0.04 

H{T)IT, 

Fig. 3.6. Calculated magnetization curves of MnAs at various temperatures {T /T{) = t) with 
r/ = 2 fixed. We see that spontaneous magnetization is cr = 0.81 at ^ = 1. The curve at ? = 1.09 
shows metamagnetic behavior clearly but the curves for 1.134 do not show hysteresis as 
well as magnetization jump 
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Fig. 3.7. Temperature dependence of //c(up) and //c(down) of metamagnetic phase transition 
calculated for MnAs using rf =2 

It is known that MnAs shows metamagnetic behavior at temperatures higher than 
^c(up) [19]. Figure 3.6 is the result calculated from (3.8) with T] = 2. The magnetiza- 
tion curves for 1 < T/Tq < 1 . 134 show a jump at the transition field Hq. These curves 
also show hysteresis as indicated by the vertical dotted lines for t = 1.09. Figure 3.7 
shows the temperature variation of transition fields for the process of increasing and 
decreasing field evaluated from Fig. 3.6. 
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Setting 7b = 7c(down) = 301 K, the critical point {T /Tq = 1.13 and Hc/Tq = 
0.029) in Fig. 3.7 corresponds to T = 340 K and // = 8.7 T. These values agree well 
with experimental results T = 345 K and H = 9T found by Zieba [20] . Magnetiza- 
tion curves measured by Grazhdankina [19] show that the critical point locates near 
350 K, which is also in good agreement with the calculated results mentioned above. 
We note that Fig. 3.7 is also applicable to metamagnetic behavior of MnAsi-^Px 
above 7n shown in the next section. 

In summary, anomalous magnetic properties of MnAs is well explained by a 
phenomenological theory in which exchange interaction between Mn is sensitive to 
crystal volume (that is, T] or j3 is large). 

Next we apply this theory to the magneto volume effect of MnAs. We start 
from (3.7). According to experiments, the jump of magnetization at the tempera- 
ture T/Tq = 1.08 is about 0.6. This value is reproduced with T] = 2. By inserting 
the value of T/Tq = 1.08 and AV /V = 0.021 (obtained from experiments) shown in 
Fig. 3.3 into (3.7), we have 

A^y^B^:roj3 = 0.0648. (3.10) 

Also, a relation 

ri = 2 = 2.21 xNkBKTop^ (3.11) 

is known from temperature dependence of magnetization. From (3.10) and (3. 1 1), we 
obtain j3 = 14 and ^ = 3.8 x 10~^^ (in CGS). Temperature dependence of magnetiza- 
tion as well as magneto volume effect atT = rc(up) (AV /V = 0.021) is reproduced 
using these values. However, the experimental value of AV /V is a summation of 
magnetovolume effect and volume change due to structural transformation from the 
Ni As-type to the MnP-type. To investigate pure magnetovolume effect, it is neces- 
sary to investigate MnAsi-^Sb^, which does not show the structural transformation 
for X > 0.1 [21, 22]. We note that the magnetovolume effect of MnAso.QSbo.i is well 
explained with the critical value T] = 1 . 

Finally, we investigate the effect of pressure on the magnetic transition temper- 
ature 7b of MnAs. In the procedure to derive (3.8), we set if = 0 and keep P finite 
instead of setting P = 0. Then we have a relation 

0 = 2aa[{l - pKP) + {2b/a)ria^] + (r/7b)a. (3.12) 

This equation shows a relation between a and T/Tq under pressure P. We set T] = 2 
and 7b = 301 K as we did for MnAs. The value of is evaluated to be fiK = 
5.32 X 10~^^ from (3.10) and (3.11). Using these values of parameters, we calculated 
temperature dependence of a under various pressures as shown in Fig. 3.8. From 
this figure, 7c (up) and 7c (down) are obtained as a function of pressure as plotted in 
Fig. 3.9. 

The calculated results shown in Fig. 3.9 agree qualitatively with the pressure de- 
pendence of 7c (up) and 7c (down) in the low pressure region in Fig. 3.10. We see 
from this figure that MnAs under pressure transforms easily to the MnP-type struc- 
ture. We also see that this transformation is of the first order accompanied with large 
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Fig. 3.8. a of MnAs under pressure as a function of temperature calculated using Kft = 5.3 x 
lO-^HCGS) 




P (kBar) 

Fig. 3.9. Pressure dependence of first order transition temperatures 7c (up) and rc(down) of 
MnAs evaluated from Fig. 3.8 



hysteresis for pressure. For example, if we apply pressure larger than P = 3.1 kbar 
(0.31 GP) at r = lOOK, the crystal remains in the MnP-type even after the pressure 
is removed. This phase diagram is similar to that of MnAsi-^Px (Fig. 3.12), which 
we will discuss in the next section. These results suggest that NiAs-type structure of 
ferromagnetic MnAs is quite unstable; it seems to be maintained by occurrence of 
the ferromagnetic order. 

The model we have shown in this section is based on assumptions that Mn atom 
in MnAs has a local spin 3/2 and that the exchange interaction depends strongly 
and linearly on the crystal volume. This phenomenological model has explained 
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Fig. 3.10. A phase diagram (pressure-temperature) of MnAs (Reproduced from [23]) 

successfully characteristic properties of MnAs. This model, however, has defects. 
One is that it is based on the localized spin model. Another is that the structural 
transformation at 7c (up) is neglected. In Part II of this book, we will show that the 
anomalous magnetic properties and structural transformation of MnAs are well ex- 
plained from the viewpoint of itinerant electrons and band structures. 

3.2.2 Various Type of Phase Transition of MnAsi_^P^ 

There are various interesting mixed compounds of MnAs such as Mni-j^Cr^As 
[24-27], MnAsi_;,P;c [15,28-30], and MnAsi_;,Sb;c (see Sect. 3.2.3) [5]. As we have 
noted, MnAs is unstable against the MnP-type lattice distortion. Substituting P for 
As realizes the MnP-type structure. On the other hand, by substitution of Sb for As, 
the NiAs-type structure is stabilized. We describe MnAsi-j^Px in this section. 

First we show a relation between magnetic moment and the a-axis of MnX 
(X = P, As, Sb, or Bi) in Fig. 3.11. In this figure, MnP and MnAso.88Po.i2 talces 
the MnP-type structure. (We note that the a-axis of the MnP-type structure in this 
figure corresponds to that of the NiAs-type structure.) We see that magnetic moment 
is sensitive to P- substitution, but the moment does not change largely as long as the 
(2-axis is longer than that of MnAs. 

Figure 3.12 shows a phase diagram of the MnAs-MnP system, which is de- 
picted using data in [3 1-34] . In this figure, the hatched region indicates the region 
in which the NiAs-type structure realizes. The compound in the vicinity of MnAs 
(x < 0.03) is a ferromagnet, which shows a first order transition at 7^. In the region 
of the MnP-phase, the double helical ordering realizes although thorough investiga- 
tion has not been performed yet. The first investigation of these mixed compounds 
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Fig. 3.11. Relation between magnetic moment and the a-axis of MnX compounds. The NiAs- 
type and the MnP-type structure are denoted by filled and open symbols, respectively 
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Fig. 3.12. A phase diagram (magnetic state and crystal structure) of the MnAs-MnP mixed 
crystal 



were carried out by Goodenough et al. [32]. In Fig. 3.13, we show M-T curves of 
MnAso.9Po.i. The Neel temperature is about 230 K. At the temperatures T < lOOK 
the double helical ordering realizes. By applying a magnetic field, the double heli- 
cal structure transforms into ferromagnetic ordering with moment 1.6/iB/Mn [35]. 
In the temperature range lOOK < T < 7n, magnetic structure is unknown though 
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Fig. 3.13. Magnetic moment [35] and inverse susceptibiUty of MnAso.9Po.i plotted against 
temperature, /igff is effective Bohr magneton (Reproduced from [32]) 
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Fig. 3.14. Magnetic moment of Mn in MnAsi-j^Px- The data for x = 0, 0.1, 0.25, and 1 were 
obtained from saturation moment measurement [14,35-37]. Others were taken from results of 
neutron diffraction measurements [28,31] 

ferromagnetic behavior of a magnetization curve has been observed. IMagnetic mo- 
ment of IMnAsi-^Px was measured by neutron diffraction [28, 31] and saturation 
moment measurements [16, 35]. The moment slightly decreases with increasing P 
component for x > 0.03 as shown in Fig. 3.14. 
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As we see from Fig. 3.13, behavior of x~ above 7n is similar to that of MnAs 
above 7c (up). In the temperature range where decreases with rising temperature, 
the MnP-type structure is changing into the NiAs-type (the parameters u, v, and 5 
approach zero) [15, 36], and the NiAs-type structure realizes above 7^. Moment 
estimated from the Curie-Weiss-like behavior of is /ieff = 2.09 iib{pa = ^^Mb = 
1.3/iB). This is close to the saturation moment observed at low temperatures shown 
in Fig. 3.14. In the range T > 7^, /igff is 4.87/iB, indicating that magnetic moment is 
as large as the moment of MnAs. 

Figure 3.13 suggests that transition between a low spin state and a high spin 
state occurs accompanied with the structural transformation. This transition was dis- 
cussed first by Goodenough et al. [32]. There is an experiment that indicates relation 
between susceptibility and crystal structures [15, 36]. Such problems have been ex- 
plained from the viewpoint of itinerant electrons by Motizuki and her co-workers as 
described in Part II of this book. 

Next we describe field-induced structural transformation of MnAso.QPo.i. We 
show temperature dependence of moment of MnAso.QPo.i in high magnetic fields in 
Fig. 3.15 [16]. As shown in Fig. 3.16, metamagnetic behavior is seen above 7n. This 
is similar to the behavior of moment of MnAs just above 7c (up). The data shown in 
Figs. 3.15 and 3.16 were measured at the high magnetic facility of Osaka University. 
Behavior similar to Figs. 3.15 and 3.16 have been observed for Mno.9Cro.9As [37]. 

In a moderate magnetic field, MnAso.9Po.i has the double helical ordering with 
low spin state below 7n ^ 230 K. In high-fields, metamagnetic transition to the 
state similar to MnAs occurs for T > 7n as shown in Figs. 3.15 and 3.16. This 
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Fig. 3.15. Temperature change of magnetization of MnAso.9Po.i under high magnetic fields. 
Data of MnAs are also plotted for comparison (Reproduced from [16]) 
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Fig. 3.16. Magnetization process (a) and transition fields (b) of metamagnetic behavior of 
MnAso.9Po.1- The open (filled) symbols are for increasing (decreasing) field (Reproduced 
from [16]) 



Magnetization 



MnAs B81 




J 




1 B8. 


MnAsci.*tP(u 


r 1 

B31 





















MnAs^ 






•V4OT 






' ^ 












M11AS0.9P0J ^ 


...i Ay^ 

■4 


13% 








1.6T * 




r 


40T 







rv = 230 300 Tc=S18 400 450 T (K) 



Fig. 3.17. Schematic figure for magnetization-temperature curves of MnAs and MnAso.gPo.i 
under high magnetic fields. Temperature range for the Ni As-type and MnP-type crystal struc- 
ture is also shown 

field-induced phase transition is expected to be a transition from paramagnetic state 
(MnP-type structure, the low spin state) to ferromagnetic state (NiAs-type, the high 
spin state). This transition is similar to the metamagnetic transition of MnAs just 
above 7c, as we have noted. Figure 3.17 shows a comparison between MnAs and 
MnAsQ.QPo.i. 

The fact that the magnetization of MnAso.QPo.i is almost the same as that of 
MnAs in the temperature region between 7n (=250 K) and Tq (=318K), as well 




as very large magnetostriction, which will be explored next, indicates that the 
metamagnetic transition of MnAso.QPo.i above accompanies the structural trans- 
formation. In MnAso.93Po.07. very large magneto volume effect AV /V = 0.15 accom- 
panied with the metamagnetic transition has been observed at T = 300 K, as shown 
in Fig. 3.18 [38]. 

The magneto volume effect shown in Fig. 3.18 is extraordinarily large but rea- 
sonably understood. Temperature dependence of lattice constants of MnAsi-^Px is 
plotted in Fig. 3.19 [15]. Unit cell volume V{= abc) for x = 0.1 at 300 K is 121 
while that of MnAs in the NiAs-type phase is 137 at 300 K (Fig. 3.3). Thus, there 
is 13% difference (as also shown in Fig. 3.17), which suggests that the large magne- 
to volume effect shown in Fig. 3.18 is reasonable. Therefore, the metamagnetic tran- 
sition shown in Fig. 3.16 is accompanied with the structural transformation (between 
the MnP-type and NiAs-type) with extraordinary volume change more than 10%. 
Metamagnetic transition accompanying large magnetostriction is possible to occur, 
because increase of elastic energy due to strain is small in the present case (this is 
because the very large thermal expansion coefficient of MnAsi-^Px results in a small 
elastic coefficient). In Fig. 3.16, we see that transition field increases with decreasing 
temperature, which is closely related to the decrease of thermal expansion coefficient 
as seen in Fig. 3.19 [15]. The thermal expansion coefficient of MnAso.QPo.i is small 
in the temperature range T <T^ and the transition field of the metamagnetic transi- 
tion shown in Fig. 3.16 shows rapid increase. At T = 231 K, metamagnetic transition 
does not occur even in 400 kOe. Results shown in Fig. 3.16 suggests that no meta- 
magnetic transition occurs for T <T^. These facts seems to be important hints to 
understand the field-induced metamagnetic transition, which occurs only above 7n. 
The phenomenological theory of metamagnetism of MnAs just above Tq shown in 
Sect. 3.2.1 may be applicable to MnAso.QPo.i- Divergence of transition field at 7n 
(Fig. 3.16) and metamagnetic transition of MnAso.QPo.i accompanied with structural 
change are complicated and interesting phenomena. Theoretical studies based on the 
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Fig. 3.19. (a), (b) Lattice parameters of MnAsi-j^Px- ^ and c are the lattice constants of the 
hexagonal lattice, (c) Unit cell volume V = abc of MnAso.9Po.i and MnAs is schematically 
plotted as a function of temperature (Reproduced from [15]) 



itinerant electron picture, as well as experiments in an ultra high field, are necessary 
to clarify the mechanism of metamagnetic transition in MnAsi-^Px- 

3.2.3 Anomalous Behavior of MnAsi_^Sbx 

The NiAs-type structure is stabilized by substituting Sb for As of MnAs. Figure 3.20 
shows a magnetic and structural phase diagram of the MnAs-MnSb system [39,40]. 
The shaded area that indicates the IVLnP-phase disappears for x :^ 0. 1 of IVLnAsi-^Sbx. 
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Fig. 3.20. A phase diagram of MnAsi-j^Sbx. The shaded area denotes the MnP-type structure. 
The compound takes the NiAs-type structure in the other area. Ferromagnetism reaUzes for all 
composition. For x < 0.1, a first order transition occurs at Tq 



For X > 0.1, magnetic phase transition at Tq is of the second order. Tq takes the 
minimum value for x = 0.35. It is also interesting that Tq and Op are significantly 
different in composition region 0. 1 < x < 0.7. Magnetic moment per Mn shows weak 
dependence on x [41] but almost constant in the range 0.15 <x< 0.7, where the 
minimum of 7c (x = 0.35) occurs. We will discuss the minimum of Tq later. 

As shown in Fig. 3.21, the a-axis increases drastically with Sb composition. The 
difference of the a-axis between MnAs and MnSb is 13%, while variation of the 
c-axis is much smaller (about 1% between MnAs and MnSb). The small jumps of 
lattice parameters at the compositions designated by the vertical dotted lines are ex- 
plained by the composition dependence of Tq shown in Fig. 3.20. If Tq is higher 
than the room temperature, the lattice parameters have an effect of spontaneous 
magnetostriction as shown in Fig. 3.22. On the other hand, when Tq is lower than 
room temperature (0.06 < x < 0.08), there is no effect of the magnetostriction. 
As shown in Fig. 3.22, magnetostriction of MnAso.ySbo.s at T = 80K is +1% for 
the a-axis and —0.7% for the c-axis. This corresponds to volume magnetostriction 
AV/y = +1.5x 10-^ 

The experimental result that Tq takes a minimum for the composition (Fig. 3.20) 
is phenomenologically explained from the viewpoint of localized model similar to 
the model of Bean and Rodbell shown in Sect. 3.2.1. 
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Fig. 3.21. Lattice constants of MnAsi_;cSbx at room temperature [39] 
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Fig. 3.22. Temperature dependence of lattice constants of MnAsojSbo.s- Spontaneous mag- 
netostriction arises below 227 K (=7c) [22] 

We assume that Tq is a function of exchange parameters Jc and Ja, which change 
linearly with lattice parameters c and a around cq and ao, respectively, as follows: 



Ja{a) =J{ao)^A{a-ao), 
Jc{c)=J{co)^B{c-co), 



(3.13) 
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where A and B are constants. Therefore, change of Tq due to lattice change is 
written as 

dTciJcM = {dTc/dJa){dJa/da)da^{dTc/dJc){dJc/dc)dc (3.14) 

and 

Tc{a, c) - Tc{ao,co) = {dTc/dJa)A{a - gq) + {dTc/dJc)B{c - cq) 

= a(a-ao)+j3(c-co), (3.15) 

where {dTc/ dJa)A and {dTc/ dJc)B are replaced by a and j3 , respectively. We insert 
the experimental lattice constants and Tq of MnSb for Tc{ao^co). By using (3.13) and 
the composition dependence of lattice constants shown in Fig. 3.21, we can calculate 
7c of MnAsi-xSbx as plotted in Fig. 3.23. We note that for x > 0.6 (where Tq is 
higher than room temperature) we extrapolated the lattice constant curve (Fig. 3.21) 
in the region x < 0.5 so as to exclude the effect of ferromagnetism. 

Results shown in Fig. 3.23 were obtained by assuming linear change of exchange 
parameters Jc and Ja against corresponding lattice constants. In calculation we set 
the coefficients a = 1360 K/A and j3 = — 3000K/A. With these parameters we 
successfully reproduced the minimum of Tq against composition of MnAsi-^Sb^ 
qualitatively. It is apparent that anomalous composition dependence of the c-axis 
gives rise to the composition dependence of 7c, though there is no theoretical foun- 
dation for the anomalous composition dependence of the c-axis. We also see from 
Fig. 3.22 that expansion of the a-axis as well as the shrink of c-axis brings about in- 
crease of 7c. This corresponds to the sign of a and j3. If elastic constants along 
the a- and c-directions are available, it is possible to evaluate quantitatively the 
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Fig. 3.23. Tq of MnAsi_;cSbx calculated using Ja and Jc whose temperature change was eval- 
uated from experimental values of the lattice constants. The minimum of Tq is caused by 
anomalous variation of the c-axis for the composition (Fig. 3.21 and also see text) 
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composition dependence of Tq as well as the temperature dependence of sponta- 
neous magnetostriction and magnetic moment on the basis of the Bean-Rodbell the- 
ory [21]. Lattice constants of MnSb and CrSb under high pressure up to 1.8 GPa 
have been measured by Nagasaki [42,43]; however, a more precise measurement 
using state-of-the-art equipment is expected. 

Next we describe metamagnetic behavior of MnAsi-^Sb^. The metamagnetism 
of MnAs just above Tq arises from the fact that ferromagnetic-paramagnetic transi- 
tion at 7c is of the first order (?] > 1). In MnAsi-^Sb^, similar metamagnetic tran- 
sition induced by applied magnetic field occurs only in the composition region of 
X < 0.08 (see Fig. 3.20). On the other hand, for x > 0.1, metamagnetic transition 
does not occur; however, metamagnetic-like behavior remains near x = 0.1, which 
means the parameter r\ in the argument in Sect. 3.2.1 is slightly smaller than 1. As 
an example, metamagnetic-like behavior of MnAsojSbo.s just above Tq is shown in 
Fig. 3.24 [22]. 

This field-induced metamagnetic-like transition does not accompany the struc- 
tural transformation between the MnP-type and Ni As-type structures. As shown in 
Appendix, metamagnetic magnetization process is generally given by the expression 

H = Aia +A2a^ ^A^a^ + . . ., (3.16) 

where H isa. magnetic field, a magnetization, and Ai etc. are expansion coefficients. 
The form of (3.16) is common to the Bean-Rodbell model and the itinerant electron 
model developed by Yamada [44], though physical meaning of the coefficients is 
different. The experimental results shown in Fig. 3.24 are well reproduced by (3.16) 




Fig. 3.24. Magnetization curves of MnAsojSbo.s- The curves are measured at T = 4.2, 77, 
180, 200, 210, 220, 225, 230, 240, 250, 255, 260, 270 K from larger magnetization. Inset: 
magnetic fields at inflection points are plotted as a function of temperature. From this figure 
we see that Tc = 230 K (Reproduced from [22]) 
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with temperature-dependent Ai, A2, and A3. However, both the above two models 
cannot explain successfully the temperature change of A2 [22]. We also note that 
there is some allowance for the experimental values of Ai, A2, and A3. 

3.2.4 Effect of High Pressure on MnAsi_^Sb^ 

Pressure effects on Tq of MnSb and MnAso.88Sbo.12 [45] and changes of magnetic 
states of MnAso.7Sbo.3 [46] have been investigated under high pressure. Tq of MnSb 
decreases by 200 K (from 590 to 390 K) by applying pressure 5GPa [47]. Under 
this pressure, the a- and c-axes decrease about 2% (from 4.13 to 4.03 A) and 4% 
(from 5.78 to 5.55 A), respectively [43]. This corresponds to the compressibility 
K = -{l/V){dV/dP) = 8.6 X 10-^^ Pa'^ (= 8.6 x lO'^^ cm^ dyn'^), which is an 
ordinary value for this kind of intermetallic compounds [42] . The elastic constant of 
MnSb is anisotropic and the compressibility is much larger than that of a material 
which has closed packed structure (e.g., four times larger than that of Cu and eight 
times larger than that of Ni). Unlike MnAso.7Sbo.3 shown in Fig. 3.22, the temper- 
ature variation curve of c-axis of MnSb shows no anomaly at Tq while the a-axis 
shows small positive magnetostriction [48]. This fact indicates that the ferromag- 
netism is stabilized by the increase of a-axis. If we are based on the localized model, 
the positive magnetostriction means that increase of a-axis results in increase of an 
exchange interaction parameter Ja connecting Mn spins within the c-plane. This is 
consistent to the behavior of MnAsi-^Sb^ with x > 0.6 shown in Figs. 3.20 and 3.21, 
where Tq increases as a-axis increases (c-axis is almost constant) with increasing 
composition x. Therefore, the decrease of Tq of MnSb by application of pressure 
is mainly caused by the decrease of a-axis. It is apparent that there is a close cor- 
relation between magnetism and lattice parameters. To clarify the mechanism from 
microscopic electronic structure is a problem left for the future. 

Next, we discuss pressure effect on Tq of MnAso.88Sbo.12. Unlike MnAs, this 
compound shows a second order transition at Tq. As shown in Fig. 3.25, Tq de- 
creases and temperature interval Ar of hysteresis for magnetization increases with 
increasing pressure. The curve for 1.2 GPa suggests a structural change caused by 
the pressure from the NiAs- to the MnP-types. 

MnAso.88Sbo.12 has an Sb composition x just above the critical composition x = 
0.1, as seen in Fig. 3.20; therefore, shows at Tq the second order transition having 
the characteristic near the first order transition. The hysteresis generated by applied 
pressure shown in Fig. 3.25 means that pressure changes the second order transition 
to the first order transition at Tq. (The large hysteresis shown in Fig. 3.25 indicates the 
first order transition though the order is unclear especially under pressures P = 0.8 
and 1.0 GPa). 

We will show that these experimental results are explained by the Bean-Rodbell 
model. We assume t] = 1 (see (3.11) and (3.12)) as this compound locates near 
the border between the first order and the second order transition. Then, by using 
parameters Tq = 300 K and Kfi = 5.3 x 10~^^ (CGS), which are the same as those 
as MnAs, we can calculate by making use of (3.12) the magnetization curves plotted 
in Fig. 3.26. 



38 3 Properties of the Compounds with Ni As-Type and Cu2Sb-Type Structures 




T(K) 

Fig. 3.25. Temperature dependence of magnetization (in arbitrary scale) of MnAso.ggSbo.n 
under high pressure. The labels 1, 2, 3, 4, 5, 6, and 7 denote the pressure 0, 0.2, 0.4, 0.6, 0.8, 
1.0, and 1.2 GPa, respectively (Reproduced from [45]) 
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Fig. 3.26. Magnetization (cr) vs. temperature curves of MnAso.ggSbo.n = 1^ = 300 K, 
PK = 5.3 X 10-11 (CGS)) under atmospheric and high pressure (P = 1 GPa). 15 K width of 
hysteresis is seen for P = I GPa 

Calculated results shown in Fig. 3.26 exhibits that pressure alters order of tran- 
sition from the second to the first, and the calculated values of A < 15 K is in fairly 
good agreement with the experimental results of A 19K (Fig. 3.25), even with- 
out parameter adjustment. It is noted that the molecular field theory based on the 
localized picture (i.e., the Bean-Rodbell model) is successful in explaining various 
properties of itinerant electron system such as present compounds. 
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Fig. 3.27. Magnetization-temperature curves of MnAso.ySbo.s under atmospheric and high 
pressure (1.2 GPa) (Reproduced from [46]) 

Next, we show pressure effect on MnAsojSbo.s in magnetic field up to 7 T 
and in quasi-hydrostatic pressure up to 1.2 GPa [46]. We have already noted that 
this material exhibits a second order transition at Tq ; however, this compound shows 
metamagnetic-like behavior just above Tq as shown in Fig. 3.24. Figure 3.27 shows 
magnetization-temperature curves underpressures 0 and 1.2 GPa. We see that Tq de- 
creases by 70 K due to 1.2 GPa pressure. (7c of lVLnAso.88Sbo.i2 decreases by 105 K). 
It is not clear in Fig. 3.27 whether the order of transition at Tq changes by pressure 
from the second to the first. However, temperature dependence of magnetization un- 
der 8 and 10 GPa shown in Fig. 3.25 suggests that phase transition under 1.2 GPa 
shown in Fig. 3.27 is of the first order It is possible to clarify whether the transition 
is of the first order or of the second by measuring the hysteresis A or temperature 
dependence of spontaneous magnetization under pressure. 

Finally, we show pressure dependence of magnetization curves at T = 4.2 K 
in Fig. 3.28. It is clearly seen that spontaneous magnetization at 4.2 K decreases 
largely by pressure; the spontaneous moment per IVLn atom M(/iB ) is found to be ex- 
pressed by M = 3 . 1 85 - 0. 126P- 0.052P^ where P means the pressure in the unit of 
GPa. A relation between spontaneous magnetic moment Ms and Tq under pressure, 
d InMs / JP = (2/3 ) d In 7c / dP, has been obtained for a weak itinerant ferromagnet by 
Takahashi [49] ; however, this relation seems to be inapplicable to a strong ferromag- 
net ]V[nAso.7Sbo.3. Pressure dependence of spontaneous magnetization at T = 4.2 K 
is explained qualitatively by the energy band calculation for changing lattice param- 
eters [46]. Note that there are many problems left for the future. 

3.2.5 Mni_^Cr^As, Mni_^Ti^As, etc. 



As IVLnAs is unstable against the IVLnP-type deformation, a mixed compound 
Mni-jCvx^^ (CrAs has the IVLnP-type structure) has the IVLnP-type structure in 
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Fig. 3.28. Pressure effect on the magnetization curve of MnAsojSbo.s at T = 4.2 K (Repro- 
duced from [46]) 

low temperature region even when x is small, and undergoes the structural trans- 
formation to the Ni As-type when temperature rises. These structural properties of 
Mni-jcCr^As are similar to those of MnAsi-^Px [50]; however, magnetic properties 
are somewhat different. Existence of two kind of 3d atoms in JMni-jcCr^As makes 
the magnetic properties a little complex. 

By substituting a nonmagnetic atom Ti for Mn of IMnAs and IMnSb, density 
of magnetic atom in the compound decreases and the NiAs-type structure is sta- 
bilized, except the vicinity of TiAs (TiSb). Tq of these mixed compounds shows 
rapid decrease with Ti composition, approaching OK for x = 1. IMni-j^TixSh be- 
comes spin glass in low temperature for the compositions around x = 0.8 [51]. Sub- 
tracting Pauli paramagnetic contribution (Xp = 1 x 10~^emug~^), susceptibility of 
]V[ni_;cTixAs obeys the Curie-Weiss law [14]. Form this Curie-Weiss type suscepti- 
bility, the paramagnetic Curie temperature Op and the effective Bohr magneton per 
Mn atom, /igff = ^Mb 1), were evaluated. By setting the g-factor 2 in /igff , we 
can calculate gSji^ • spin) = Pa (atomic moment). Figure 3.29 illustrates 0p, 7c, 
and Pa as a function of Ti composition. ]V[ni_;cTixSb also shows the similar behavior 
[14]. It is notable that Pa is almost independent on Ti composition. There has been 
report on temperature dependence of lattice constants and pressure effect on Tq of 
]V[ni__;cTixAs (0 < X < 0.9) [21]. An itinerant electron theory is desired to explain the 
properties shown above. 

3.2.6 Magnetic Refrigeration Using MnAs and the Related Compounds 



If an applied magnetic field largely reduces entropy of a material, this material 
can be a candidate for a refrigerant of magnetic refrigeration. Enhancement of 
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Fig. 3.29. Tq, 0p, and Pa(= g^/ie/Mn, see text) of Mni_;cTixAs. Mni_;cTixSb shows similar 
behavior (Reproduced from [14]) 

magnetization due to applied magnetic field means increase of magnetic order or 
reduction of magnetic entropy. For a ferromagnet, large increase of magnetization 
due to applied field occurs in the temperature region around 7c; just above and be- 
low 7c for the second order transition and only just above 7c for the first order 
transition (cf. Figs. 3.4 and 3.6). A ferromagnet that shows a first order transition 
(e.g., MnAs) shows metamagnetic transition just above 7c. Metamagnetic transition 
field He increases with rising temperature (see Fig. 3.7). In Fig. 3.30, magnetiza- 
tion curves above and below 7c are show schematically for a ferromagnet with the 
second or the first order transitions at 7c. 

If we measure or calculate magnetization curves at various temperatures, it is 
possible to evaluate corresponding magnetic entropy using the expansion formula 
of a (3.3).^ In an adiabatic condition where total entropy is constant, reduction of 



Temperature dependence of decrease of magnetic entropy A5m by application of external 
field H is calculated from difference of specific heats in // = 0 and H = H. ASm is also 
calculated by integrating the Maxwell relation {dS/dH)T = {dM/dT)H of thermodynam- 
ics as 

rH 

ASm= / (dM/dT)HdH, 
Jo 
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Fig. 3.30. Schematic magnetization curves above and below Tq 
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Fig. 3.31. Temperature dependence of magnetic entropy (a schematic figure for the case of a 
second order transition). See text for Ti 

magnetic entropy due to applied field enhances entropy of lattice vibration, bringing 
about increase of temperature by AT. This is called magnetocaloric effect. A material 
that has large AT is suitable for magnetic refrigerant. 

We explain briefly how AT arises below [52]. Figure 3.31 shows schematic tem- 
perature dependence of magnetic entropy Sm with and without magnetic field. 

Total entropy Sj of a magnetic material is Sm + 5'l lattice entropy). Assum- 
ing that the lattice specific heat Cl is constant in the range 7i < T (7i is a temperature 
just below 7c), 5'l is calculated as 

SLiT)-Si^{n)= rcLdr/r = CLin(r/ri). (3.17) 

This equation means that the lattice entropy is a monotonically increasing function 
of temperature. Taking account of temperature dependence of in Fig. 3.31, we 
can express Sj schematically as shown in Fig. 3.32. 

where M is magnetization under a magnetic field. We need to measure magnetization in 
various magnetic fields to obtain ASm- By integrating dM/dT with respect to H at respec- 
tive temperatures, we can calculate temperature variation of ASm- We note that M in this 
argument means magnetization in a state in which all the atomic moments align in the 
direction of H, that is, a state without effect of anisotropy and magnetic domain. 
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Fig. 3.32. Schematic curves of total entropy St vs. temperature T inH = 0 and H = H 

If we apply H adiabatically (i.e., is unchanged) to a magnetic material at 
temperature T = T\ (// = 0), temperature rises up to T = T2 by Ar. The resulting 
state at r = 72 (// = H) is equivalent to a state in which H is isothermally applied 
to a material at T = 72 (77 = 0). In the former adiabatic process, entropy of spins 
decreases (because spins align along the applied field) and lattice entropy increases, 
and then temperature of the material increases. On the other hand, dissipation of heat 
keeps the temperature unchanged in the latter isothermal process. 

The process of magnetic refrigeration is explained next. By applying 77, the state 
(72,77 = 0), which is taken as an initial state here, changes to the state (72,77 = 77). 
In this process, generated heat is released to outside. Adiabatic removal of magnetic 
field brings the state (72,77) to (ri,77 = 0). In this process, the temperature of the 
material lowers from 72 to T\ , namely the material absorbs heat from outside (i.e., air 
around). Repetition of the similar process is the magnetic refrigeration. The magnetic 
material must show large magnetocaloric effect within a finite temperature range 
because the temperature of refrigerant decreases during the refrigeration process. 
A magnetic material that has large Ar in a wide temperature range is favorable. From 
Fig. 3.32 we see that dS^/dT = dS^/dT ASu/AT for T > Tq. Since dS^/dT = 
Cl/^ as seen in (3.17), we have 



Therefore, AT is proportional to reduction of magnetic entropy due to applied field. 
For T <Tc- AT, dSj/dT = dS^/dT + dSu/dT = C^/T + Cu/T is approximately 
satisfied and then 



where Cm (77 = 0) is the magnetic specific heat in 77 = 0. From (3.18) and (3.19) 
we see that Ar for T < Tc is expected to be smaller than Ar for T > Tq because of 
the difference between denominators of (3.18) and (3.19). Wada et al. investigated 
MnAs and MnAsi-^Sb^ as magnetic refrigerant [53-55]. We show some of their data 
in Fig. 3.33. 



AT - TASm/Cl - 7cA5m/Cl. 



(3.18) 



Ar^7cA%/(CL + CM(77 = 0)), 



(3.19) 
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Fig. 3.33. (a) Magnetocaloric effect AT of MnAs evaluated indirectly from temperature 
variations of total entropy *St and magnetization, (b) Decrease of magnetic entropy — A^m 
of MnAsi-j^Sbjc measured by applying magnetic field 5 T (Reproduced from [53]) 
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Fig. 3.34. Schematic figure to show the mechanism of refrigeration by pressure for 
MnAsi_j(:Sbjc with small x. This figure is for a first order transition at Tq 



In Fig. 3.33b, A5'm at if = 5 T was evaluated from temperature change of mag- 
netization in various magnetic fields, whereas in (a), ASu was determined indirectly 
from the data shown in (b) and temperature dependence of entropy in zero field [55]. 
IMethod of direct measurement of the magnetocaloric effect and experimental data 
for Gd metal have been reported together with analysis based on the molecular field 
theory in [56]. 

We have discussed in Sect. 3.2.4 that Tq of IVLnAsi-^Sbx with small Sb content 
decreases remarkably by application of pressure. It may be possible to apply the 
compounds for refrigeration using pressure instead of magnetic field for the refrig- 
eration process. As seen from Fig. 3.34b, removal of pressure isothermally reduces 
magnetic entropy. This is equivalent to application of magnetic field isothermally. 
It is interesting to clarify which is more practical. 
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3.3 MnSb and MnBi 

We have already shown basic properties of these ferromagnetic compounds in 
Sects. 1.4 and 2.1. As MnSb has already been mentioned in Sect. 3.2.4, we explain 
MnBi in this section. As shown in Table 1.2 and Fig. 3.11, MnBi takes the NiAs- 
type structure with lattice constants larger than MnSb, and magnetic moment per Mn 
atom is 3.8/iB- MnBi shows a first order transition with temperature hysteresis 15 K 
at 7c = 628 K [57]. Figure 3.35 shows temperature dependence of spontaneous mag- 
netic moment, where data not only for low temperature phase (LTP) MnBi but also 
for Mni.ogBi (quenched high temperature phase: QHTP) are plotted. The LTP shows 
a first order transition at 7c, which is accompanied with crystallographic transition to 
the nonstoichiometric Mni.ogBi. The QHTP has smaller moment and lower 7c than 
those of the LTP MnBi. These properties are similar to those of Mni+j^Sb; however, 
X is limited to 0.08 for Mni+;cBi. 

The QHTP MnBi separates into Mn and Mn-Bi liquid at the peritectic tempera- 
ture of 446° C. As MnBi has large magnetic anisotropy with magnetic moment par- 
allel to the c-axis (anisotropy constant K c^2x 10^ ergcm"^ = 2 x 10^ Jm~^) [57], 
there were attempts for applications to permanent magnet [58]. However, this ma- 
terial is not suitable for a magnet because of the relatively small magnetization and 
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difficulty in preparation of compound. This material may be utilized for other type of 
application such as magnetic memory. There has been no theoretical study to explain 
the peculiar properties of MnBi. 



3.4 CrAs and Related Compounds 
3.4.1 Anomalous Magnetic Transition of CrAs 

In contrast with the ferromagnetic Mn-pnictides with MnX-type, the compounds 
with CrX (X = As or Sb) type take generally antiferromagnetic ordering. As mag- 
netic properties of CrAs are especially interesting, we describe mainly CrAs and its 
related compounds. 

As explained in Sect. 2.1, CrAs takes the double helical magnetic ordering below 
7n = 272 K [59], and undergoes a first order transition at 7n to paramagnetic state in 
high temperature region above 7n. We show susceptibility of CrAs and CrAso.QPo.i in 
Fig. 3.36 (Ido, unpublished). Similar data for CrAsi-^Px have already been reported 
by Selte [60]. 

As seen in Fig. 3.36, susceptibility of CrAs increases with increasing temperature 
without anomaly at 7n, suggesting strong itinerant character. It is known that CrAs 
with the MnP-type crystal structure takes the NiAs-type structure in the temperature 
region higher than 7^ = HOOK, and the magnetic susceptibility obeys the Curie- 
Weiss law in this temperature range, which is not included in Fig. 3.36 [50]. For 
CrAsi-xSbx (x < 0.6), the Curie-Weiss law in the region of T > 7^ is clearly observed 
as Ti decreases due to substitution of Sb for As [61]. For example, CrAso.QSbo.i has 
transition temperature 7^ = 780K and shows the Curie-Weiss behavior of suscep- 
tibility in r > 780K. From the Curie-Weiss law, moment of Cr in CrAso.QSbo.i 
is evaluated from the effective Bohr magneton to be 2.8/iB, which is equal to the 
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Fig. 3.36. Susceptibility of CrAs and CrAso.9Po.1- The magnetic order occurs in CrAs below 
Tn, while no magnetic order occurs in CrAso.QPo.i down to 0 K 
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moment of Cr in CrSb determined by neutron diffraction measurements (see also 
Sect. 3.5). It is interesting that susceptibility of CrAs in the NiAs-type phase obeys 
approximately the Curie-Weiss law expected for Cr^+ (i.e., 3/iB)- This phenomenon 
will be shown and discussed later in Fig. 3.42. 7n of CrAs was determined by neu- 
tron diffraction measurements (Fig. 3.37) [59, 62]. IVLagnetic phase transition at 
is of the first order accompanied with discontinuous change of lattice constants. We 
will discuss in detail later the magneto volume effect of CrAs. 

It has been found from specific heat measurements that latent heat accompanied 
with magnetic phase transition is 5.3 x 10^ Jmol~^ (=0.98calg~^) [63,64]. Temper- 
ature dependence of electric resistivity is similar to that of a gap-type antiferromagnet 
such as Cr and Cri-^^clVLn^ (x < 0.25), whose resistivity increases first with decreas- 
ing temperature below 7n and decreases next on further cooling [63] . It is difficult to 
measure resistivity of CrAs accurately because of cracks arising due to large volume 
change at 7n. This is similar to the situation in IVLnAs. 

In Fig. 3.38, magnetic moment per Cr atom for some Cr pnictides is plotted as a 
function of length of the a-axis. This figure is similar to Fig. 3.11, where IMn-moment 
is plotted against a-axis length for the IVLn-pnictides. We note again that the a-axis 
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Fig. 3.38. Relation between moment per a Cr atom and the a-axis at 4.2 K. Filled and open 
symbols denote the MnP-type and NiAs-type structures, respectively. This figure has been 
drawn by making use of the data in [65] and [66, 67]. 



for the orthorhombic IMnP-type structure in Fig. 3.38 corresponds to the a-axis of the 
hexagonal NiAs-type structure. The reason why we take only a-axis as a parameter 
for jU(Cr) in Fig. 3.38 is that the c-axis variation from 5.45 A of CrSb to 5.36 A of 
CrP is negligibly small compared with the large a-axis variation from 4. 13 A of CrSb 
to 3.11 A of CrP. 

Figure 3.38 shows apparently that magnetism of CrAs is very sensitive to the 
change of a-axis, which will be discussed in detail in the next section. It may need to 
point out that there is little change of /i(Cr) around a = 3.75 A where the structural 
change occurs. This may mean that the distortion to the IMnP-type structure does not 
affect so much. This is sharply different from IMnAso.QPo.i, whose moment exhibits 
large change accompanied with the structural change as shown in Fig. 3.14. 

3.4.2 Critical Lattice Constant of CrAsi_^Px and Cri_^M^As 
(M = Mn, Ni, etc.) 

We see the anomalous magnetic properties of CrAs more clearly in the behavior of 
mixed compounds. A typical example is a Cr As-based compound in which small 
fraction of As is replaced by P. Both electron number and crystal structure remain 
unchanged by the P substitution for As. According to [60], the double helical or- 
dering with the first order transition at 7n ^ 250 K disappears for CrAsi-^Px with 
X > 0.03 [60], which is also roughly expected from Fig. 3.38. This is a very interest- 
ing phenomenon whose mechanism has not been clarified yet. 

We show and discuss experimental data of magneto volume effect of CrAsi-^Px 
in the following part of this section. Figure 3.39 shows temperature change of lattice 
constants of some of the mixed compounds derived from CrAs [68]. As we noted 
in Sect. 1.3, the A-, B-, and C-axes here correspond to the c-, a-, and /7-axes of the 
hexagonal NiAs-type structure, respectively. 
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Fig. 3.39. Temperature change of lattice constants of CrAs and related compounds. Note that 
the B- and C-axes correspond to the a and Z^-axes of the hexagonal NiAs-type structure, re- 
spectively (Reproduced from [68]) 



In Fig. 3.39 data for the A(c)-axis, where (c) means the A-axis corresponds to the 
c-axis, are not plotted because the discontinuous change of A (c) -axis at 7n is much 
smaller than those of B{a)- and C(/?)-axes. For example, in the case of CrAso.98Po.02, 
the discontinuous change of the B- and C-axes are 5.5% and 1.1%, respectively, 
while that of the A(c)-axis is only 0.4%. Therefore, as generally seen in Fig. 3.39, 
magneto volume effect takes place mainly in the 5 (a) -direction. It is noted that the 
change of the 5-axis is exceptionally large. We may regard the change in the C{b)- 
direction at 7n, which has an opposite sign to that of 5 (a) -axis, is only induced by 
the large change of the 5(a)-axis. To obtain the normal thermal expansion curve 
of the 5(a)-axis, temperature dependence of the 5(a)-axis of CrAso.9Po.i, which 
shows no magnetic ordering at T = 0 K, was measured for the temperature from 80 
to 330 K. As shown by the dotted curves in Fig. 3.40, the part of the curve below 80 K 
was extrapolated on the basis of the specific heat (H. Ido and T. Kamimura, unpub- 
lished) and the Griineisen relation [69], which says that specific heat is proportional 
to a thermal expansion coefficient. By applying the normal thermal expansion curve 
presumed for CrAso.QPo.i, magneto volume effect at T = OK has been estimated 
for CrAs and CrAso.98Po.02 to be 4.6% (0.16 A) and 5.6% (0.19 A), respectively, 
as shown in Fig. 3.40. These values are extremely large. As no magnetic ordering 
occurs in CrAso.9Po.i, this compound shows no magneto volume effect at T = OK. 
We see from Fig. 3.39 that magnetic ordering of CrAs and CrAso.98Po.02 is stabilized 
by the large lattice distortion in the 5 (a) -direction. We can also say that increase of 
elastic energy due to magnetostriction is larger than decrease of magnetic energy in 
CrAso.9Po.i ; therefore, this compound has no magnetic ordering. Therefore, in CrAs 
and the related compounds there seems to be a close relation between the occurrence 
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Fig. 3.40. Spontaneous magnetostriction of CrAsi_;cPx and critical 5-axis Bq (Reproduced 
from [68]) 



of magnetic order and the length of B{a)-dod^. As indicated in Fig. 3.40, the criti- 
cal length of 5 ((2) -axis, Bq, is about 3.37 A, namely magnetic ordering occurs in the 
compounds with B larger than Bq, and no magnetic ordering occurs if 5 < ^c- CrAs 
with B larger than Bq has moment 1.1 jl^ and 7n = 250 K. Figure 3.40 indicates that 
magnetic ordering of CrAs is expected to disappear by more than 5.6% compression 
of the 5(a)-axis. Figure 3.40 also suggests that, if spontaneous magnetostriction were 
fixed at the value just below the first order transition point 7n, the Neel temperature 
would be more than 350 K. 

As seen from Fig. 3.39, 7n (a first order transition) increases due to Mn sub- 
stitution for Cr of CrAs. In Fig. 3.41a, 7n determined by thermal expansion curve 
is plotted against composition of various substitution elements. We see that de- 
creases (increases) when the 5 (a) -axis shrinks (increases) due to the substitution. 
On the other hand, in Fig. 3.41b, 7n is plotted not for composition but for the length 
of 5(a)-axis just above 7n (=250) K, which is almost close to the 5-axis length at 
r = 0 K in the absence of magnetostriction. 

As seen from Fig. 3.40, the important point we see from Fig. 3.41b is that there 
is a critical length of the 5-axis Bq = 3.38 A, which is independent of the kind of 
substituting atom. To construct a magnetic phase diagram, the crystallographic tran- 
sition temperatures Tt between the MnP-type and the NiAs-type, which are taken 
from Fig. 3.42, are also shown in Fig. 3.41b. 

In Fig. 3.41b, the labels "Pauli para" and "Curie-Weiss type para" indicate the re- 
gions that magnetic susceptibility behaves as those in Fig. 3.36 and as those above 7^ 
in Fig. 3.42, respectively [61]. In Fig. 3.42, we show moment of a Cr atom /icr (evalu- 
ated from the Curie-Weiss law of susceptibility), G)p, etc. As CrAso.QSbo.i has higher 
7n and lower Tt compared with CrAs (see Fig. 3.36 and [50]), the Curie-Weiss-type 
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behavior of susceptibility appears more clearly especially for the compound with 
X = 0.3. We will discuss again the results shown in Fig. 3.42 in Sect. 3.4.4. 

3.4.3 Phenomenological Theory for the First Order Transition of CrAs 

The phenomenological theory for the first order transition of MnAs (see Sect. 3.2.1 
and Appendix) is applicable to CrAs. We note that physical meaning of the terms of 
free energy expansion is different. When magnetization of space variation is written 
as M(r) = Mq Qxp{iqr) (q = In /2.67b along the /7-axis in the case of CrAs. See 
Fig. 2.1), we may assume molecular field AMq within the molecular field theory. 
Then we have magnetic energy = — {1/2) AM^. By setting an assumption that 
the molecular field coefficient A depends linearly on crystal volume as A = Ao[l + 
j3 (y — Vb) /Vb)] . we have an expression 

£n, = -(i/2)Ao[i+i3(y-yo)/yo)]cT', (3.20) 

where we wrote Mq = a, which corresponds to the relative magnetization in the Bean 
and Rodbell theory for MnAs (see Sect. 3.2.1). As mentioned in Sect. 3.2.1, elastic 
energy due to lattice distortion is written as 

E,i = {l/2K)[{V-Vo)/{Vo)]^ (3.21) 

and magnetic entropy per unit mass of the crystal is given by 

S{a) = NkB [S{0) + aa^ + ba"^ + co"^ + . . .] . (3.22) 

By setting derivative of E = Em-\- Eq\ with respect to V zero, we have a relation 

(V - yo)/yo = {Kl5Ao/2)a\ (3.23) 

Inserting (3.23) into E (= £^ni + ^ei). we have an expression of free energy G, which 
is minimized with respect to V as 

G/NkBTo = -tS{0) - a{t - l)a^ - b{t -11)0"^- cta^ + . . ., (3.24) 

where Tq, t, and T] are expressed by the following equations (3.25) and (3.26), re- 
spectively, 

-{Ao/2aNkB) = To and T/To = t (3.25) 

and 

T] = -{a^/2NkBb)Kl3^To (> 0). (3.26) 

We see that (3.26) is the same as (A8) in Appendix; therefore, a first order tran- 
sition occurs when T] > 1 , namely if the parameter expressing volume dependence 
of magnetic energy of CrAs is large enough for T] > 1 to hold, the magnetization 
a {= Mq) disappears discontinuously at 7n. As (3.23) shows that the spontaneous 
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magnetization is proportional to j3, the experimental results shown in Fig. 3.39 
clearly show that magnetic energy is very sensitive to crystal volume (mainly to 
the 5 (a) -axis). Results of neutron diffraction measurements revealed that a of the 
double helical CrAs at 7n is about 0.7, which is close to the value of MnAs. Spon- 
taneous magnetostriction {V - Vo)/Vo at 7n is 0.022 for CrAs [68] and 0.021 for 
MnAs. If we express (3.23) as {V — Vq) /Vq = Aa^, the coefficient A of CrAs is very 
close to A of MnAs; however, there is contribution of structural transformation to the 
magnetostriction at Tq in the case of MnAs. Therefore, we may regard the value of A 
for CrAs as twice as that of MnAs. The magnetism accompanied with the first order 
transition at 7n of CrAs, which is closely correlated to the critical length Bq of the 
5-axis, is a very interesting problem, and is expected to be clarified in the future in 
the view point of itinerant electron picture. 

3.4.4 CrAsi_^Sb^ 

As noted in Sect. 2.1, CrSb takes the NiAs-type structure. In contrast to CrAsi-^Px, 
substitution of Sb for As in CrAs reduces 7^ (i.e., transition temperature between 
the MnP-type and the NiAs-type structure). See also Fig. 3.42. Figure 3.43 shows a 
magnetic and crystallographic phase diagram of the CrAs-CrSb system constructed 
from temperature dependence of lattice constants obtained by X-ray diffraction mea- 
surements [61]. 

The rise of in the region near CrAs seen in Fig. 3.43 is probably arising from 
the competition between magnetic energy and elastic energy due to magnetostric- 
tion; decrease of magnetic energy and increase of elastic energy is caused by the 
spontaneous magnetization in CrAs; however, increase of lattice constants due to the 
Sb-substitution reduces the increase of the elastic energy, resulting in stabilization 
of magnetic ordering and higher T^. By this mechanism the experimental results in 
the region near CrAs is explained. A first order transition occurs when there is com- 
petition between the increase of elastic energy and the decrease of magnetic energy. 
Therefore, if there is no such competition, order of transition becomes the second, 
which is realized in the Sb-composition x between 0.1 and 0.4 as seen in Fig. 3.43. 
7n seems to take a minimum around x 0:^0.5. This behavior is similar to that in the 
MnSb-MnAs system; however, may have more complex origin. 

Next, we explain 7/ for x = 0.6 as an example in Fig. 3.43. On heating the 
lattice constants shows at T/, a large jump as Ac = 0.15 A (Ac/c = 0.026), and 
Aa = — 0.07 A (Aa/a = —0.018). We note that we defined T/ as a middle point of 
temperature region about lOOK wide where a high-temperature phase and a low- 
temperature phase coexist. Corresponding to the jump of the lattice constants at T/, 
susceptibility shows steep increase at 7/ as shown in Fig. 3.42. The first order transi- 
tion at r/ is quite interesting but left unsolved. 

We pay our attention again to the temperature dependence of susceptibility of 
CrAso.9Sbo.i shown in Fig. 3.42 in Sect. 3.4.2. Behavior of susceptibility in the tem- 
perature region below 7^ is very similar to that of CrAs (see Fig. 3.36). On the other 



54 



3 Properties of the Compounds with Ni As-Type and Cu2Sb-Type Structures 




0 0.5 1.0 

CrAs X CrSb 



Fig. 3.43. A magnetic and crystallographic phase diagram of CrAs-CrSb system. Thick lines 
indicate that the first order transition takes place. T/ corresponds to a transition that accom- 
panies discrete change of the c- and a-axes; however, detail of this transition has not been 
revealed yet (Reproduced from [61]) 



hand, susceptibility for T >Tt obeys the Curie-Weiss law with Cr moment 2.8/iB 
(estimated using g = 2) and paramagnetic Curie temperature Op = — 77 K. These 
values agree well with the values of CrSb [70]. It is noted that CrAso.QSbo.i in the 
temperature region of T < 7^, where the crystal structure is with IVLnP-type, is similar 
to CrAs, while it is similar to CrSb in T > 7^, where the crystal structure is of the 
Ni As-type. 

Finally, in Fig. 3.44, we show magnetic moment per Cr atom in the CrAs-CrSb 
system determined by Kallel et al. from their neutron diffraction measurements at 
T = 4.2 K [65]. The discrepancy around x = 0.6 is caused by the differences in mag- 
netic order and crystal structure; on the CrSb side, the compounds have an anti- 
ferromagnetic order similar to that of CrSb with the NiAs-type structure, while on 
the CrAs side the compounds take the double helical ordering similar to that of CrAs 
with the IVLnP-type structure. However, change of magnetic moment at around x = 0.6 
is not large, which is in contrast to the case of IMnAsi-j^Px, where large difference of 
moment between the NiAs-phase and the IVLnP-phase occurs (see Fig. 3.14). 
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CrAs X CrSb 

Fig. 3.44. Magnetic moment of a Cr atom in CrAsi-j^Sbx obtained by neutron diffraction 
measurements. For x < 0.5, the double helical ordering and the MnP-type structure realize 
(i.e., similar to CrAs). For x > 0.5, antiferromagnetism and the NiAs-type structure realize 
(similar to CrSb) [65] 

3.5 CrSb 

CrSb, which takes the NiAs-type structure, is an antiferromagnet with 7n = 718K 
and a simple spin arrangement shown in Fig. 3.45. 

Magnetic moment of a Cr atom is equal or close to Cr^+ moment (3.0/iB) with 
g = 2. As shown in Fig. 3.38, the (2-axis of CrSb is much longer than that of CrAs, 
while difference of the c-axis between CrSb and CrAs is small (5.65 A for CrSb and 
5.45 A for CrAs). Such a large a-axis of CrSb brings about narrow width of the 3d 
band and enhancement of intra-atomic exchange interaction, consequently the large 
moment. As we mentioned in Sect. 3.4.2, temperature dependence of susceptibility 
of CrSb behaves like an antiferromagnet with two sublattice [71]. This behavior is 
quite different from that in the temperature region < T < Tt for CrAs. We can 
obtain Op = — 200K and Cr-moment 2.9/iB for CrSb from the Curie-Weiss type 
susceptibility together with assumption of ^ = 2. The value of 2.9/iB thus obtained 
is in good agreement with 3.0 and 2.7 ji^ by the neutron diffraction measurements 
shown in Fig. 3.45. As compared with CrAs, CrSb seems to have a strong localized 
character. This kind of tendency is also seen in IVLn-pnictides such as IVLnBi, IVLnSb, 
IVLnAs, and IVLnP (see Fig. 3.11). The relation between a-axis and magnetic moment 
shown in Figs. 3.38 and 3.11 seems to suggest a way of approach to understand the 
magnetism of 3 J-pnictides. There is relevant discussion in Sect. 4.2 of Part II. 

Resistivity of CrSb above 7n shows temperature dependence like a semiconduc- 
tor [72]; however, it is unsure whether such behavior of resistivity arises from an 
intrinsic property of CrSb or slight crystal decomposition at high temperatures such 
as Sb-precipitation as seen in IVLnSb (Sect. 1.4 and Fig. 1.6). 
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Fig. 3.45. Antiferromagnetic spin ordering of CrSb. Moment of a Cr atom has been reported 
as S.OjUb [65] or 2.1 jib [70] 



3.6 CrP 

CrP takes the IMnP-type structure. As shown in Fig. 3.38, the 5(a)-axis is shorter by 
0.5 A than that of CrAs, which broadens the 3d band. Therefore, density of states 
at the Fermi level is small and no magnetic ordering appears. Susceptibility of CrP 
has minimum at about 200 K [73]. IMagnitude and temperature dependence of the 
susceptibility are very similar to those of CrAso.QPo.i as shown in Fig. 3.36. 



3.7 Properties of Cu2Sb-Type Compounds 

We have already shown magnetic structures of these compounds in Fig. 2.2 in 
Sect. 2.1. Band structures are shown in Part II of this book. In this section, we mainly 
discuss the very interesting magnetic properties of ]V[n2-xCrxSb. 

3.7.1 Antiferromagnetic-Ferrimagnetic Transition of Mn2-xCr^Sb 
and Kittel's IModel 

It has been well known that compounds ]V[n2Sb whose Mn or Sb is slightly replaced 
by V, Cr, Co, Cu, etc. or Ge, As etc. undergo with increase in temperature a first order 
transition from antiferromagnetic state in low temperature side to ferrimagnetic state 
in high temperature side. Among these compounds, ]V[n2-xCrxSb shows the most 
typical transition with very sharp change of magnetization. Figure 3.46 shows tem- 
perature dependence of spontaneous magnetization measured for single crystalline 
samples, and Fig. 3.47 shows the lattice constants as a function of temperature [74]. 
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Fig. 3.46. Temperature changes of spontaneous magnetization of Mn2_xCr;cSb. Mn2Sb is 
ferrimagnetic. Other samples with x except for x = 0.013 show the antiferromagnetic- 
ferrimagnetic transition (Reproduced from [74]) 
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Fig. 3.47. Temperature change of lattice constants of Mn2_xCrjcSb (Reproduced from [74]) 



In Fig. 3.46, we see the first order antiferromagnetic-ferrimagnetic transition for 
the samples with x in the region 0.05 < x < 0.16, whereas we also see an interme- 
diate magnetic ordering with a small magnetization for the samples with x < 0.03. 
The small magnetizations of about 3 emug~^ is generally observed for all samples 
in the antiferromagnetic region due to small amount of impurity of ferromagnetic 
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Fig. 3.48. Schematic picture of the antiferromagnetic-ferrimagnetic transition of 
Mn2_jcCr;cSb. Mn atoms occupy two different sites Mn(I) and Mn(II). Moments at the Mn(I)- 
site and Mn(II)-site are antiparallel. Mi and M2 are summation of moments in three adjacent 
layers. When Mi and M2 are parallel (antiparallel), we refer it as ferrimagnetic (antiferromag- 
netic) state 



MnSb phase. As shown in Fig. 3.47, the c-axis shows a large jump at the first order 
while the small jump of the a-axis may be attributed to a kind of reaction to the large 
jump of the c-axis. 

Next we explain the change of spin alignment between the antiferromagnetic 
and ferrimagnetic states. The antiferromagnetic-ferrimagnetic transition can be ex- 
plained by the exchange inversion model developed by Kittel [75], where the sign 
of exchange parameter Jc shown in Fig. 3.48 changes its sign at a critical length of 
the c-axis. The Kittel model is a phenomenological model, and the stability of the 
ferrimagnetic and the antiferromagnetic states has also been discussed from the view- 
point of itinerant electrons by Motizuki and co-workers (see Part II). We see from 
Figs. 3.47 and 3.48 that the c-axis length is closely related to namely it is easily 
seen that the contraction of the c-axis brings about the change of Jc from positive to 
negative values. This situation is illustrated schematically in Fig. 3.49. 

We explain on the basis of the Kittel model for the antiferromagnetic- 
ferrimagnetic transition referring to Figs. 3.48 and 3.49. Among energies relevant to 
the transition, magnetic energy 

^ex = -AMi-M2 (3.27) 

is directly related to the exchange parameter Jc. In this equation. Mi and M2 are sub- 
lattice magnetizations as shown in Fig. 3.48, and A is a coefficient of the molecular 
field given by 

A = {2JcZ)/{NgiiB), (3.28) 

where we note that J^ as explained in Fig. 3.48, is an average exchange parameter of 
several different exchange parameters Ji connecting Mn-spins belonging to Mi and 
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Fig. 3.49. (a) Temperature change of the c-axis with hysteresis is schematically shown. The 
experimental data in the vicinity of Tt in the process of rising temperature are shown in 
Fig. 3.47. The dashed line (cq) means the normal thermal expansion curve without the sponta- 
neous magnetostriction (c — cq). (b) Jc (see Fig. 3.48) plotted schematically against the c-axis 
length. We assume that sign of Jc changes at the critical length Cc 



M2, and Z the number of Mn-Mn spin pairs. We consider here only two exchange 
parameters Ji and J2 indicated by the dashed lines in Fig. 3.48, and we have a re- 
lation JcZ = 4/1 + 2/2. For Z = 6 this relation becomes Jc = (2/3)/i + (l/3)/2- It 
is important that the coefficient A in (3.28) is proportional to J^ that is, the summa- 
tion of J I and J2 in the form mentioned above. (Note that Kittel did not mention A 
in detail.) Considering that sign of Jc or A changes at c = c^, we can expand A as 
follows, 

A{c)=A\cc){c-Cc), (3.29) 

where A' (cc) means the first derivative. Substituting (3.29) into (3.27) and replacing 
A'(cc) by B, we have 

^ex = -BM^ {c-Cc) cos e, (3.30) 

where 0 is a relative angle between Mi and M2, and M{= Mi = M2) is magnitude 
of magnetization of each sublattice. 

Next we consider elastic energy induced by the spontaneous magnetostriction 
given by the following equation: 

£ei = (1/2) [cn{el^ + e^yy) + css^^J + cne^cxeyy + ci^{eyyezz + ^zz^xx), 
= {\/2)[2cn{^a/af^c^^{^c/cY]^cn{^a/af 

+ 2ci3(Aa/a)(Ac/c), (3.31) 

where exx = eyy = l^a/ a= {a — ao)/aQ and 6^1 = ^c/c = {c — cq)/cq are shear (see 
Fig. 3.49a),andQ/sarestiffnessconstants[76].Bysetting — (Aa/a)/(Ac/c) =/?(> 0), 
the elastic energy is written as 



E,i = [C33/2 + (cii + cn)R^ - IRcu] {^c/cf = K{c - co)^ (3.32) 
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with K defined by 

K = [C33/2 + (cii + cn)R^ - 2Rcu]/cl (3.33) 

From Fig. 3.47 together with Fig. 3.49a, we see that Aa/a is smaller than Ac/c, so 
that R 0.3. Thus, neglecting small terms in (3.33), we find K is mainly contributed 
by the term Kc^^/2cq. Therefore, it is seen from (3.32) that the elastic energy in the 
present case is approximately proportional to square of magnetostriction along the 
c-axis. We note that different origin c — Cq has been used for Eqx in (3.30). Magnetic 
entropy is unchanged during the transition at 7^, because this is a transition between 
different magnetic orderings.^ We can thus omit the entropy term in free energy. 
As a result, we can write the free energy as 

F = ^ex + ^el = -Bhf{c - Cc) COS Q ^ K{c - c^f . (3.34) 

From a condition dF /dc = 0, we have 

c-co = {B/2K)M^cose. (3.35) 

Inserting (3.35) into (3.34), we have an expression of F which is minimum with 
respect to c as 

F{e) = -B{co-Cc)M^ cose -{B^M'^/4K)cos^e=acose-bcos^e, (3.36) 

where b>0 and a is negative (positive) for cq > Cc (cq < Cc). 
We can write as 

F/b = Z) cos e - cos^ e (3.37) 

with 

D = a/b = -{4K/BM^){co-Cc). (3.38) 

Using (3.37) in Fig. 3.50, we calculate and plot F/b with various values of Z) as a 
function of 6. Equation (3.38) shows Z) < 0 for cq > Cc and Z) > 0 for cq < c^. 

We see from this figure that 6 = n (antiferromagnetic state) corresponds to the 
smallest free energy when Z) > 0 (cq < c^) as denoted by the dotted circle. As in- 
dicated by the vertical arrow in the right hand side, D decreases as cq increases 
with rising temperature (this is normal thermal expansion plotted by the dashed 
line in Fig. 3.49a). At higher temperatures, D becomes negative. However, antifer- 
romagnetic state (6 = n) still remains stable as long as Z) > —2. It is apparent 
from Fig. 3.50 that a first order transition to ferrimagnetic state (6 = 0) occurs at 
D = —2. Similarly, the ferrimagnetic state in the process of decreasing temperature 
shows a transition to the antiferromagnetic state at Z) = 2. From (3.38), we see that 
co = Cc^ {BM^/2K) at Z) = -2 and cq = - {BM^/2K) at Z) = 2. This situation is 
summarized in Fig. 3.51. 



There is little change of entropy due to change of the molecular field. 
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Fig. 3.50. F/b (F: free energy) plotted as a function of 6 (an angle between the sublattice 
moments) for various values of D(= -{4K/BM^){co - Cc)) 



c, Co ferri c 




' 7^(K) 

Fig. 3.51. Normal thermal expansion (cq) and magnetostriction (c — cq) = ±{BM'^ /2K) in the 
vicinity of transition temperature. AT denotes the temperature hysteresis. Transition occurs 
when Co = Cc ± {BJVfl /2K) is satisfied (see text). The dotted lines indicates length of the c-axis 
to be measured 

From a relation Cc + {BM^/2K) = Cc[l + a(Ar/2)], Ar is written as 

AT = {B/ccK)M^/a, (3.39) 

where a is a thermal expansion coefficient of cq. We see from (3.39) that width of 
hysteresis is proportional to M^. 

We have explained qualitatively the first order antiferromagnetic-ferrimagnetic 
transition of Mn2-xCr;cSb in the vicinity of the first order transition temperature. 
Next, we compare the result mentioned above with experimental data. From (3.35), 
the jump of the c-axis at T = 7^ (up or down), Ac/c, is approximately written as 




A[(c - Co)/ Co] - A[(c - co)/cc] = {B/Kcc)m\ (3.40) 
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Fig. 3.52. Relation between the square of sublattice moment and the jump of the c-axis at 
the first order transition temperature Tt. (Drawn by using experimental data shown in Figs. 3.46 
and 3.47) 
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Fig. 3.53. Temperature hysteresis of magnetization for Mn2_xCrjcSb (x = 0.07). AT ^ 15 K 

and A(c — cq)/ Cc, which is simply expressed by Ac/c, at T = 7^ can be determined 
from experimental results shown in Figs. 3.46 and 3.47. We note that magnetization 
of a sublattice M has been assumed to be a half of the jump of magnetization at Tt. 

Figure 3.52 shows experimental result corresponding to (3.40). The linear rela- 
tionship between Ac/c and shown in Fig. 3.52 indicates that (3.35) is satisfied. 
From the slope of the dotted line in Fig. 3.52 and (3.35), we can evaluate {B/Kcc) as 
1.9 X 10~^. Applying this value to (3.39), we obtain 

Ar = 1.9x 10"^MVa. (3.41) 

Now we show temperature variation of magnetization in Fig. 3.53 to show an ex- 
ample of the hysteresis Ar. Note that it is somewhat difficult to determine the cor- 
rect value of AT from experiments, because the hysteresis depends on the rate of 
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temperature change and the sample quahty, etc. In addition, a large coexistence re- 
gion of low temperature antiferromagnetic phase and high temperature ferrimagnetic 
phase has been found by X-ray diffraction measurements. 

From Fig. 3.53, we see that M 15 K. Inserting AT = 15K and 
150(emu/g)^ (see Fig. 3.46) in (3.41), we can estimate as a 2.5 x 10~^ 
though direct measurement of a is not available. This value of a is larger than a of 
a similar compound Mn2Sbo.875Aso.125 [77], however, usual one in the intermetallic 
compounds. 

We note that there is no other model to explain the first order antiferromagnetic- 
ferrimagnetic transition of Mn2-xCrxSb. The model described above is based on the 
assumption that the exchange parameter Jc connecting different sublattices is sensi- 
tive to length of the c-axis and changes its sign at the critical length c^. As shown in 
Fig. 3.49, Jc is the summation of Ji and J2. If Ji and J2 have different sign and com- 
parable magnitude, sign of Jc may change due to even small change of the c-axis. It 
is noted that Ji and J2 mean exchange interactions between rather distant Mn atoms; 
the nearest distance is 5.05 A for Ji and 3.94 A for J2 as shown in Fig. 3.48. It is 
expected to investigate Ji and J2 on the basis of band structures. 

3.7.2 Magnetic Transition of Fe^-^Mn^ As (a 0^ 2) 

As for compounds written as Fe^-xMn^As, both Fe2As and Mn2As show antiferro- 
magnetic ordering. For a = 2.1 and 1.25 < x < 1.5, however, ferrimagnetic structure 
occurs in low temperature region. This ferrimagnetism undergoes a first order tran- 
sition to antiferromagnetic ordering with rising temperature [78]. 

3.7.3 Layered Ferromagnets MnAlGe and MnGaGe 

These compounds have an easy axis of magnetization along the c-axis. In these com- 
pounds Mn atoms occupy only the Mn(I)-site shown in Fig. 3.48, and Al (Ga) and 
Ge occupy the Mn(II)-site and the Sb-site, respectively. Mn atoms thus align within 
a two-dimensional plane (001), and the distance between the neighboring planes 
is as long as the c-axis (5.933 A for MnAlGe). These compounds therefore have 
a magnetically two-dimensional character. Probably the RKKY (Ruderman-Kittel- 
Kasuya-Yoshida) interaction and super-exchange interaction via Al, Ga, and Ge are 
responsible for the magnetic interaction between Mn atoms in the adjacent layers. 
Curie temperature Tq and magnetic moment per ?>d atom in Mni-j^M^AlGe are plot- 
ted against X, respectively, in Figs. 3.54 and 3.55 (M = 3 J metals) [79,80]. By mailing 
use of the band of MnAlGe [8 1], we can explain the reduction of moment for the case 
of M = Fe, which has one more electron within the rigid band picture. The case of 
Cr substitution is also explained in the similar way. The prominent increase of Tq by 
80 K due to 15% Cr substitution, which is shown in Fig. 3.54, contrasts to that the 
Cr-substitution decreases the averaged moment of ?>d atom. Such increase of Tq may 
be caused by change of the inter-layer interaction as well as intra-layer interaction. 
MnGaGe shows properties similar to MnAlGe [82]. Band structures and magnetic 
properties of these compounds are also discussed in Part II of this book. 
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Fig. 3.54. Tq of Mni_;cMjcAlGe (M = 3d atom) plotted as functions of composition (Repro- 
duced from [79]) 

3.7.4 Application of the First Order Transition of Mn2-xCr^Sb 

As shown in Fig. 3.53, the hysteresis of the antiferromagnetic-ferrimagnetic transi- 
tion is rather large. A composite material made of compounds Mn2-xCrxSb with two 
slightly different Cr compositions is expected to show a magnetization-temperature 
curve as shown schematically by solid bold lines in Fig. 3.56. 

To explain this figure, let us start from an antiferromagnetic state at the initial 
temperature Tq, which is not shown in the figure, however, between T2 and 7i . When 
the temperature is raised, the initial magnetization Mq is zero at temperatures T <T\, 
and then the nonmagnetic state changes discontinuously to the magnetic state with 
residual magnetization M\ when temperature is raised up to T\. M\ remains un- 
changed in the temperature range of 7i < T < T2. With further heating, residual 
magnetization Mi similarly increases to M2. 

It is thus possible to store three values Mq, Mi, and M2 at one position (bit) 
by different three heating temperatures T's with T2 < T < T\, T\ < T < T2, and 
T2 <T. Thus, a composite material made of the three compounds with different Cr- 
composition can, in principle, store four values. It may be possible to store more 
data using a composite material made of compounds with continuously different 
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Fig. 3.55. Magnetic moment per 3d atom in Mni-j^^Mj^AlGe (M = 3d atom). Pgff is effective 
moment estimated from the Curie-Weiss constant, and /Ia is magnetic moment evaluated from 
spontaneous moment (Reproduced from [79]) 
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Fig. 3.56. Magnetization-temperature curves of composite material made of two compounds 
having different Cr compositions. Ti and T2 (or and are the first order transition tem- 
peratures of the compound in heating and cooling processes, respectively. {Ti — T[) etc. is 
temperature width of hysteresis 



Cr compositions. This mechanism is superior to an usual magnetic memory, which 
stores two values at a position (bit) corresponding to plus/minus magnetization. We 
need to pay attention that the hysteresis curves in Fig. 3.56 are idealized. Actual 
magnetization curves show jump of magnetization within a finite temperature range 
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as shown in Figs. 3.53 and 3.46. It is necessary to reduce this temperature range 
by, for example, homogenization of Cr composition [83]. This compound system 
Mn2-xCrxSb may be applicable easily for temperature controller because large vari- 
ety of transition temperatures can be obtained by changing the Cr composition. 
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Appendix: Magnetic Transition and Free Energy 



Free energy G of a magnetic material can be written in the Landau type expansion 
form as 



where (J is a magnetic moment of a ferromagnet (or a sublattice moment for an an- 
tiferromagnet) normahzed by the value at T = 0 K. When an external magnetic field 
is applied, Zeeman energy —oH is added. This expansion form is irrelevant to the 
nature of electrons, localized or itinerant; however, physical meaning of the coeffi- 
cients depends on the character of electrons and theoretical model. In this Appendix, 
we show an expansion form of magnetic entropy, which appears in the Bean-Rodbell 
theory of magnetic phase transition of MnAs, up to the sixth order (7^. Further, we 
show expression of G{a) and its temperature dependence. 

Let us consider a paramagnetic material consisting of A'^ particles with spin j in 
a magnetic field H. In this case free energy is written as G = —gjjl^aHN — TS{a). 
Replacing S{a) by Nk^S{a), we have G/Nk^T = -{gjii^o /kBT)a - S{a) = 
—aa — S{g). Applying the minimum condition (9G/(9cj = 0 to the final equality, 
we obtain 



On the other hand, a is related to a through the Brillouin function with spin j, 

Bj{a), as 



G = G(0) + a2G^ + (24Ct'^ + aeO^ + . . ., 



(Al) 



a = -dS{G)/dG. 



(A2) 



G = Bj{a) =Aa-\-Ba^-\-Da^-\-... 



(A3) 



where A, 5, Z), 



are given by 



A = (l/3)[(2j+ 1)2 - l]/(2j)^ B = -(l/45)[(2j+ 1)4- l]/(2;r, 
Z)=(l/945)[(27+l)^-l]/(2;f.... 



(A4) 
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By expanding S{a) in (A2) with respect to O", we have 

S{a) = S{0) + aa^ + /?aV ca^ + • • • 

= ln(2 J + 1 ) + acj^ + Z?(j'^ + ccj^ + • • . (A5) 

As the number of microscopic states, in the case of O" = 0, is given by Z = (2j + 1)^ 
and S = kB InZ, it is reasonably understood that S{0) becomes ln(2j + 1 ) . By inserting 

(A5) into (A2), we have a = — {2aa -\-4ba^ -\- 6ca^ H ) • By inserting this relation 

into (A3) and setting the coefficient of each a" term zero, we have the expansion 
coefficients a, /?, c, . . . as 

a = -l/A<0, b = B/4A'^<0, 

c = -{\/2)(B^/A^) +Z)/6A^ < 0, (A6) 

where A, and D are given by (A4). We insert (A5) into the expression of magnetic 
entropy (3.2) of Sect. 3.2.1, and then we put a condition of a minimum energy with 
respect to volume dG/dV = 0 into (3.2) again. Thus, we obtain the free energy G 
also minimized by crystal volume. By setting T/Tq = t (the relative temperature), 
P = 0, and = 0, we have 

G{a)/NkBTo = -dn(2 j + 1) - a{t - l)a^ - b(t -n)a^- cta^ + • • • , (A7) 

with 

T] = -{a^/2b)NkBKToli^ > 0. (A8) 

Figure 3.57 displays [G((j) - G(0)]/A^^b^o calculated by (A7) with parameters t] = 2 
and j = 3/2, which were chosen so as to be applicable to MnAs. 



0.08 




Fig. 3.57. Free energy with parameters j = 3/2 and r] = 2 vs. normalized moment a drawn 
for various relative temperatures 
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On heating, the minimum point of free energy disappears at temperature t = 
1.083 as seen in Fig. 3.57, which means that magnetization a disappears discontinu- 
ously due to the absence of energy barrier. Therefore, this is of the first order and also 
illustrated in the G-T curve with T] = 2 in Fig. 3.5. While in the process of decreas- 
ing temperature from about ^ = 1.083, the peak of free energy disappears at ^ = 1. 
This correspond to sudden appearance of magnetization. We note that the curves in 
Fig. 3.57 are somewhat inaccurate in the region cr 1 because of truncation of the 
expansion formula of entropy. (Numerical calculation is necessary to evaluate en- 
tropy without the expansion form.) The coefficient of the a'^-term of free energy has 
an important meaning. When T] > 1, this coefficient is negative in the temperature 
range 1 < ^ < t] . As coefficients of O"^ and terms are positive for ^ > 1 , the free 
energy does not have minimum at finite G unless the coefficient of the G^ term is 
negative. Therefore, the first order transition takes place only for T] > 1. 



Part II 



Itinerant Electron Theory 



4 



Electronic Band Structure and Magnetism 
of NiAs-Type Compounds 



4.1 Band Calculation 

Extensive studies on band structures of NiAs-type transition-metal compounds 
started in the middle of 1980s by the group of Motizuki, who is one of the present 
authors, of Osaka University and the group of Haas of the University of Groningen. 

For the band calculations, the group of Osaka University employed the self- 
consistent augmented plane wave (APW) method, while the group of the University 
of Groningen used the self-consistent augmented spherical wave (ASW) method. In 
spite of the difference in the method used for calculations, the general features of 
calculated bands are similar to each other. In the APW calculations carried out by 
the group of Osaka University, muffin-tin approximation was used for atomic poten- 
tial. For exchange and correlation interaction, the local (spin) density approximation 
(LDA or LSDA) was adopted within the formula of Gunnarsson and Lundqvist (GL). 
The equation of von Barth and Hedin (HB) for the LDA was also used in some cal- 
culations. In this book, we show mainly the results calculated by the group of Osaka 
University. 

In recent years, we have applied more sophisticated methods. For example, we 
carried out band calculations by using the self-consistent full-potential linearized 
augmented plane wave (FLAPW) method to obtain more accurate band structures. 
Kulatov and his co-workers carried out band calculations by using the self-consistent 
linear-muffin-tin orbital (LMTO) method. Nakada and Yamada recently calculated 
bands by using the LMTO method with atomic sphere approximation (ASA); this 
method which employs localized basis functions is convenient for the calculations 
for a system including many atoms in the unit cell. 

Scalar relativistic effects are usually included, however, spin-orbit interaction 
(SOI) is neglected in many band calculations. The SOI removes the degeneracy in 
the dispersion curves and also induces mixture between up-bands and down-bands, 
resulting in reduction of magnetic moment of ferromagnetic state. Furthermore, the 
SOI in the band calculations is needed to investigate optical processes. Recently, the 
SOI has been included in the calculations by Kulatov et al. 
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Generalized gradient approximation (GGA) has been used to improve the LDA. 
The GGA has enabled us to see difference of total energies between different struc- 
tures. For example, total energy calculation for Co As having the NiAs- and the MnP- 
type structures has been carried out by Nakada, and his results by the full-potential 
LMTO method with the GGA show good agreement with experiments. 



4.2 Band Structures and Optical Properties 
4.2.1 Pnictides: MnAs and MnSb 

Nonmagnetic Band of MnAs and MnSb with NiAs-Type Structure 

Figure 4.1a, b shows the crystal structure and the first Brillouin zone of the NiAs- 
type structure, respectively. Figure 4.2a shows the dispersion curves of nonmagnetic 
MnAs calculated by the self-consistent APW method with the muffin-tin approxi- 
mation and the LDA [1]. The two bands in the low-energy region consist of As-4^ 
orbitals. Above the gap, there are 16 mixing bands consisting of As-4/7 and Mn-SJ 
orbitals. Figure 4.2b shows the density of states (DOS) of the p-d mixing bands. 
It is seen from this figure that there are two remarkable peaks: the low-energy re- 
gion is contributed by bonding orbitals of A^-Ap and Mn-SJ, and the peaks in high- 
energy region by antibonding orbitals mainly consisting of 3 J states of Mn. As seen 
in Fig. 4.2b, the width of the antibonding p-d mixing bands is about 5 eV (or 0.4 Ry), 
which suggests that d states of MnAs are not localized but are of itinerant character. 
Density of states of MnSb has features similar to those of MnAs, however the band 
width is slightly smaller than that of MnAs. 




Fig. 4.1. (a) The NiAs-type crystal structure and (b) the first Brillouin zone 




As shown in Fig. 4.2b, the Fermi level is near the large peak of the density of 
states, which is contributed mainly by the d orbitals of Mn. The density of states at 
E^f, p{Ey^), of MnAs is larger than that of NiAs-type compounds with 3d elements 
other than Mn, which is consistent with the fact that both MnAs and MnSb are fer- 
romagnetic below the Curie temperatures 7c = 318 K and 537 K, respectively. 

These densities of states calculated by the APW method agree well with the 
bands of MnAs and MnSb calculated by Sandratskii by using the Green's function 
method [2]. Recently, Nakada and Yamada have calculated nonmagnetic bands of 
MnAs having the NiAs-type structure by the FLAPW method with the GL formula 
and the LMTO-ASA method with the BH formula. Their results shown in Fig. 4.3a, 
b are rich in fine structure, however the general features of the bands are similar to 
the bands obtained by using the APW method shown in Fig. 4.2b. 

Ferromagnetic Bands of MnAs and MnSb with the NiAs-Type 

Among the NiAs-type compounds, only MnAs, MnSb, and MnBi are ferromagnetic 
as shown in Part I. We have calculated ferromagnetic bands using the self-consistent 
APW method. The muffin-tin approximation for atomic potential and the LSDA 
within the GL formula for the exchange and correlation interaction were used. 

The density of states of the p-d mixing bands of ferromagnetic MnAs is shown 
in Fig. 4.4. Due to spin-polarization, the mixing bands of As-4p and Mn-3d orbitals 
split in the up- spin and the down- spin subbands, whereas no polarization occurs 
in the 4^ bands in the low-energy region. In Fig. 4.4, the broken and dot-dashed 
curves show components of Mn-3J and As-4/? orbitals in the muffin-tin spheres, 
respectively. Comparing Fig. 4.4 with the nonmagnetic bands shown in Fig. 4.2, we 
see that the p-d mixing bands are not shifted rigidly, but deformed significantly due 
to the spin-polarization. 
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a 




Energy (Ry) 

Fig. 4.3. Density of states of nonmagnetic MnAs calculated using the (a) FLAPW method and 
(b) LMTO-ASA method (K. Nakada, Dr. thesis) 

We evaluated total magnetic moment from the calculated band as 3. 1 /iB/formula. 
This value is somewhat smaller than the observed magnetic moment 3. 4 jUb /formula 
(see Part I). Magnetic moment in the muffin-tin spheres of a Mn-site and an As-site 
is 3.12/iB/formula and —0.1 5 /iB/formula, respectively, hence most of the magnetic 
moment is attributed to Mn atoms. In addition to the magnitude of the total magnetic 
moment, the calculated values of magnetic moments at Mn- and As- sites correspond 
well to the results of neutron diffraction measurements, therefore we see that mag- 
netic moments in MnAs are well understood from the band theory. 

Ferromagnetic band of MnSb with the NiAs-type structure is similar to the 
band of MnAs except that the width of p-d mixing band is slightly smaller than 
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Fig. 4.4. Density of states of ferromagnetic MnAs calculated by the APW method. The broken 
and dash-dotted lines are Mn-3J and As-4/? components, respectively 



that of MnAs. According to the band calculation, the total magnetic moment is 
3.2/iB/formula; magnetic moment inside the muffin-tin sphere at Mn- and Sb-sites 
are 3.41/iB/formula and — 0.12/iB/formula, respectively. These values of magnetic 
moment agree fairly well with those measured by Yamaguchi et al. [3]. 

Results of our band calculation for MnSb agree well with the bands calculated 
by Coehoorn et al. by using the self-consistent ASW method [4], namely both dis- 
persion curves and density of states in these two works agree well with each other, 
and the magnetic moments by Coehoorn et al. are 3.3j(iB at a Mn-site and —0.06 jUb 
at an Sb-site. Density of states at the Fermi level p(£'f) is obtained by us to be 37.0 
(states/Ry unit cell), which agrees well with 35.8 (states/Ry unit cell) calculated by 
Coehoorn et al. These values are also in good agreement with 32.6 (states/Ry unit 
cell) evaluated from measurements of specific heat. 

X-ray photoelectron spectroscopy (XPS) measurements for MnSb, CoSb, and 
NiSb were carried out by Liang and Chen [6]. Coehoorn et al. compared the XPS 
spectrum with their ferromagnetic bands (the 5^ band and the p-d mixing band) of 
MnSb. As shown later in Fig. 4.17 in Sect. 4.2.5, there is good correspondence for 
the peaks of the XPS spectrum and density of states (see also Sect. 4.2.5). 
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a 




-0.3 -0.2 -0.1 0 0.1 0.2 0.3 



Energy (Ry) 

Fig. 4.5. Density of state of ferromagnetic MnAs. (a) FLAPW method, (b) LMTO-ASA 
method [5] 

Nakada and Yamada calculated ferromagnetic band of MnAs with the NiAs- 
type structure by using both the FLAPW method and the LMTO-ASA method [5]. 
Calculated densities of states are shown in Fig. 4.5a, b. Compared with the band 
obtained by the APW method, these densities of states are rich in fine structures, 
however, their general features are similar to those calculated by the APW method. 

Nakada and Yamada also evaluated contributions of five Mn-d orbitals and three 
As-p orbitals to the APW wavef unctions. Their results are shown in Fig. 4.6, where 
five J-components for density of states and their values (or numbers of electron be- 
longing to the respective components) are depicted as a function of electron energy. 
Contributions of five J-components to charge and magnetic moment are also shown 
in Table 4. 1, where the values are normalized so that the summation of contributions 
from the five d orbitals is 100%. These calculations revealed that Mn-J orbitals are 
not extended in the c-plane, but have rather spherical distribution. They also plot- 
ted charge distribution in the (0001) plane and (1100) plane which include Mn-Mn 
bonds. The results also confirmed that charge density of Mn-J orbitals has spherical 
character. Furthermore, they clarified by plotting charge density in the (1120) plane 
including Mn and As atoms that the Mn-As bond is the strongest bond. The charge 
distributions are shown in Fig. 4.7. 
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Fig. 4.6. Density of states of ferromagnetic MnAs calculated by the FLAPW method, (a) 
Integrated DOS and (b) components of each d orbitals (K. Nakada, private communication) 



Table 4.1. Contribution rates of d orbitals for charge and magnetic moment, which are inte- 
grated up to the Fermi level, of MnAs 
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They also investigated the effect of crystal volume change on magnitude of mag- 
netic moment. They carried out band calculations in two cases: one is the case of 
fixed-axial ratio c/a (i.e., the bond angle Mn-As-Mn is kept constant) and another 
is the case of unfixed-axial ratio (i.e., the bond angle is variable). They found that 
magnitude of magnetic moment is significantly affected by the change of bond an- 
gle. This result indicates that magnetic moment is sensitive to lattice distortion in the 
direction of the a-axis. On the other hand, effect of lattice distortion along the c-axis 
on magnetic moment is small, though Mn-Mn distance in this direction is smaller 
than the distance along the a-axis. These results of band calculation are consistent 
Wi\h the experimental data seen, for example, in Figs. 3.3 and 3.11. 

Nonmagnetic Band of MnAs with MnP-Type Structure 

MnAs shows a structural transformation to the MnP-type structure at 7^ = 3 1 8 K (see 
also Fig. 3.2). We calculated band structure of MnAs with the MnP-type structure 
by using the self-consistent APW method. The band structure deforms due to the 
lattice distortion accompanied with the structural transformation, which results in 
decrease of p(£^f)- In Sect. 4.3, we will show in detail that anomalous behaviors of 
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Fig. 4.7. Charge distribution of ferromagnetic MnAs (a) (0001) plane, (b) (1100) plane, (c) 
(1120) plane. In each figure, (i), (ii), and (iii) show majority spin, minority spin, and moment, 
respectively [5] 
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paramagnetic susceptibility and thermal expansion observed for MnAs (see Fig. 3.2) 
are well explained by considering spin fluctuation effect and the deformation of the 
band structure. 



4.2.2 FeAs, CoAs, and NiAs 

It is known that Co As shows structural transformation at 7^ = 1 , 250 K from the MnP- 
type structure (low-temperature phase) to the NiAs-type structure (high-temperature 
phase). FeAs has the MnP-type structure, and NiAs has the NiAs-type structure.^ 
FeAs has the double-helical magnetic ordering lower than = 11 K, whereas CoAs 
and NiAs are nonmagnetic at all temperatures (see Table 2.7). 

Nonmagnetic Bands of CoAs and NiAs with the NiAs-Type Structure 

We calculated band structures of CoAs and NiAs with the NiAs-type structure by 
using the self-consistent APW method within the muffin-tin potential approxima- 
tion and the LDA. Density of states of the p-d mixing band of CoAs is shown in 
Fig. 4.8. Density of states at the Fermi level is p(£'f) = 46 (states/Ry unit cell) for 
CoAs and 31 (states/Ry unit cell) for NiAs. These values are much smaller than 156 
(states/Ry unit cell) of MnAs, which is consistent with the fact that CoAs and NiAs 
are nonmagnetic but MnAs is ferromagnetic. 

We also depicted Fermi surfaces of CoAs and NiAs, and then found that CoAs 
has two hole surfaces around the A-axis and these surfaces are expected to have good 
nesting for a displacement of wave vector FM. Since this wave vector FM is relevant 




ENERGY (Ryd) 

Fig. 4.8. Density of states of nonmagnetic CoAs with the NiAs-type structure calculated by 
the APW method 



^ We have related argument in Sect. 4.4. 
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to deformation into the MnP-type structure, the Ni As-type structure of Co As has an 
instabihty to the MnP-type structure. This fact corresponds well to the high 7^ of 
Co As. On the other hand, the structural transformation to the MnP-type structure 
will not occur for NiAs because no nesting of Fermi surfaces is expected [7] . We 
have related discussion in Sects. 4.4 and 4.7. 



Nonmagnetic Bands of FeAs and CoAs with the MnP-Type Structure 

Figure 4.9 shows DOS curves for the p-d mixing band of nonmagnetic Co As 
and FeAs with the MnP-type structure. These compounds have characteristic DOS 
curves, namely, the densities of states show abrupt increase just below the Fermi level 
[7]. We will show in Sect. 4.3 that anomalous temperature dependence of paramag- 
netic susceptibility of CoAs and FeAs can be well explained by the spin fluctuation 
theory including such feature of the densities of states. 

It is known that total energies evaluated from band structures make it possible 
to know which crystal structure is most stable at OK. Recently, Nakada and co- 
workers calculated nonmagnetic bands of CoAs with the NiAs-type and the MnP- 
type structures by using the full-potential LMTO method with the GL formula for 
the exchange and correlation interaction. From the calculated total energies they 
found that the MnP-type structure is more stable than the NiAs-type structure for 
CoAs. This result is consistent with observation that CoAs remains in the MnP- 
type structure up to 1,250 K. They obtained also similar result by calculation of full- 
potential LMTO method with the GGA term. 




Fig. 4.9. DOS curves calculated by the APW method for nonmagnetic CoAs and FeAs having 
MnP-type structure 



4.2 Band Structures and Optical Properties 



85 



4.2.3 CrSb, CrAs, and CrP 

These compounds take the NiAs-type or the MnP-type structure [8]. As shown in 
Table 4.2, the lattice constants of these compounds generally decrease as the an- 
ion changes in the order Sb, As, and P. In particular, change of the a-axis is re- 
markable. The CrX compounds show various magnetic properties with variation of 
anion atom. CrSb is an antiferromagnet with Cr moment of 3.0/iB [9, 10]. CrAs 
takes a double-helical magnetic ordering, and the magnetic moment of Cr is 1.7/iB, 
which is only about a half of Cr-moment of CrSb [11]. CrP is a paramagnetic com- 
pound whose magnetic susceptibility is almost independent on temperature [12]. The 
double-helical magnetic ordering of CrAs has been known to disappear by a few per- 
cent substitution of P atoms for As [13]. The experimental results mentioned earlier 
suggest that there is a close correlation between the crystal structure (or lattice pa- 
rameters) and magnetic properties in the CrX compounds. 

To investigate the correlation between crystal structure and magnetic properties 
of the CrX compounds, we carried out first-principle band calculations by using the 
full potential linear muffin-tin orbital (FP-LMTO) method [14]. The crystallographic 
parameters in Table 4.2 were used in the calculations. As for the local density approx- 
imation (LDA) for exchange and correlation interaction, we have used the formula 
of Vosko, Wilk, and Nussair [15]. We carried out iterative calculations using 125 
A:-points in the 1/24 of the first Brillouin zone so that muffin-tin potential converges 
within 10~^ mRy. 

Figure 4.10 shows densities of states of nonmagnetic CrP, CrAs, and CrSb. From 
these densities of states we see following three features: 

1. Total band width increases in the order CrSb, CrAs, and CrP. 

2. Hybridization between Cr SJ-orbitals and anion orbitals takes place in both high- 
and low-energy regions, but not in the intermediate region. 

3. Width of the large DOS peak near the Fermi level in CrSb, where there is little 
Cr-X mixing broadens and splits roughly in two groups as X goes from Sb to P. 

Table 4.2. Crystallographic parameters of CrX (X = P, As, or Sb) at the room tempera- 
ture. Lattice constants a, b, and c are common to both NiAs- and MnP-type structures (see 
Sect. 1.3). u, V, w, and X denote displacement parameters [8] 
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As shown in Table 4.2, the smallest distance between Cr atoms in the c-direction, 
which is equal to c/2 (~2.7A) in the case of the NiAs-type structure, is almost com- 
mon to all the CrX compounds. On the other hand, length of the a-axis which cor- 
responds to the distance between Cr atoms within the c-plane, depends considerably 
on X of the CrX compounds. This difference in the length of the a-axis influences 
the bonding character in the c-plane, resulting in the characteristics of the density of 
states especially in the vicinity of as mentioned earlier. 

We evaluated contribution of five J-orbitals to the total density of state. The re- 
sults for CrSb and CrAs are shown in Fig. 4.11. As shown in Table 4.2, CrSb and 
CrAs take different crystal structure. The length of the a-axis of CrAs is much smaller 
than that of CrSb. We discuss relation between such difference of the a-axis length 
and magnetic properties. As shown in Fig. 4.11, the density of states of CrSb has a 
sharp peak near the Fermi level. This part consists of mainly d{xy) and d{x^ — y^) 
components which extend in the c-plane and do not mix with p states of Sb as shown 
in Fig. 4.10. On the other hand, states near the Fermi level of CrAs consist of not 
only d{xy) and d{x^ — y^) but also d{xz) and d{yz) orbitals. Because of the short 
a-axis of CrAs compared with that of CrSb, width of d{xy)- and d{x^ — j^)-bands 
become larger for CrAs, and then the Fermi level is pushed down to the low-energy 
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Fig. 4.11. Contribution of 3d orbitals on density of states of CrSb and CrAs 



side of the peak of the density of states. The atomic displacement of As which is 
denoted by v in Table 4.2 due to the MnP-type distortion gives rise to lowering the 
energy of d{xz) and d{yz) components, however, these states remain in the energy 
region higher than the Fermi level. The low-energy (high-energy) peak is bonding 
(antibonding) band consisting of d{xz), d{yz), <i(3z^ — r^) states and p state of anion 
as seen in Fig. 4.10. 

It is important that the change of the a-axis, which induces change of d{xy) and 
d{x^ — y'^) bands, plays an important role to explain the variety of magnetic properties 
seen in CrSb, CrAs, and CrP. 
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Fig. 4.12. Density of states of Cr for one of the two magnetic sublattices of antiferromagnetic 
CrSb (in unit of states/Ry atom) 



From the density of states of CrP shown in Fig. 4.10, p(£'f) is evaluated as 
27.3 (states/Ry f.u.), which is 1/2 and 1/4 of the value of p(£'f) estimated from 
measurements of specific heat [16] and magnetic susceptibility measurement [17], 
respectively. 

Figure 4.12 shows density of states for one of the two magnetic sublattices of an- 
tiferromagnetic CrSb, which gives Cr magnetic moment of 2.7/iB- This value agrees 
well with experimental values 2.7 /Ib [10] and 3.0 /ie [9]. Since there is little p-d 
mixing in the region near the Fermi level as shown in Fig. 4.10, there is strong po- 
larization of Qx-d states; most of the up-spin states lies above the Fermi level. This 
fact indicates that Cr- J states in CrSb has a considerable localized character probably 
because of the large length of a-axis. It is also notable that the down-spin band has a 
gap-like shape near the Fermi level. 

Calculated total energies for paramagnetic, ferromagnetic, and antiferromagnetic 
states of CrSb are plotted in Fig. 4. 13 as a function of crystal volume. As seen in this 
figure, the band calculation indicates that the antiferromagnetic state is most stable. 
This is consistent with observation that CrSb is an antiferromagnet. 

To understand magnetic properties of CrAs (double helix, /icr = 1-7/iB, 7n = 
250 K, see Sect. 3.4 for detail), it is expected to investigate wavevector-dependent 
susceptibility %{ci). This is an important problem left for the future. 



4.2.4 Chalcogenides: CrTe, CrSe, and CrS 

Compounds between Cr and a chalcogen (S, Se, or Te) also take the Ni As-type 
structure. CrTe is a ferromagnet with 7c = 340 K, whereas CrSe and CrS take an- 
tiferromagnetic ordering below 7n = 285 K and 460 K, respectively. Effects of high 
pressure on these compounds have attracted much interest. From electron-spin res- 
onance (ESR) measurements, Shanditsev pointed out that ferromagnetism of CrTe 
disappears at the critical pressure 2.8 GP [18]. Kanomata and his co-workers inves- 
tigated pressure dependence of Curie temperature 7c and magnetic moment [19]. 
They pointed out that the ferromagnetic state is not the ground state under high 
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Fig. 4.13. Total energy of CrSb as a function of unit cell volume. The total energies for para- 
magnetic state, ferromagnetic state, and antiferromagnetic state are plotted. Vexper denotes 
experimental value of unit cell volume 



pressure. According to observation by Ishiduka and Eto, the ferromagnetism of CrTe 
disappears at 7 GP and a structural transition to the MnP-type structure occurs at 
13GP [20]. 

To explain the results of these high-pressure experiments, we have theoretically 
studied the possible magnetic transitions by investigating pressure dependence of 
band structures. We will discuss the pressure effect in detail in Sect. 4.5, and in 
this section we describe the results of band structures calculated by Nakada and 
co-workers; they investigated dependence of bond length (Cr-Cr) and bond angle 
(Cr-Te-Cr) on charge distribution and magnetic moment in ferromagnetic CrTe. 

On the basis of the FLAPW band, calculations for ferromagnetic CrTe, Nakada 
decomposed the APW wavefunction into the components from three p orbitals of Te 
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and five d orbitals of Cr. The results are shown in Fig. 4.14a, b and Table 4.3. The 
density of states and integrated density of states are shown in Fig. 4.14b and 4.14a, 
respectively. 

Contribution from each d orbital to charge and magnetic moment in CrTe is 
shown in Table 4.3, where the values are normalized so that the summation of con- 
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Fig. 4.14. (b) Density of state of ferromagnetic CrTe calculated by the FLAPW method, (a) 
Integrated DOS component from each d orbital (K. Nakada, private communication) 



Table 4.3. Contribution rates of respective d orbitals to charge and magnetic moment in CrTe 
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tributions from five d orbitals is 100%. We see that all d orbitals have almost equal 
contribution to charge and magnetic moment, which means that charge of the d states 
of a Cr atom distributes spherically as also seen in MnAs. This result that charge dis- 
tribution of Cr is spherical has been confirmed by depicting charge distribution in the 
(0001) and (1100) planes which include Cr-Cr bond. They also found strong Cr-Te 
bond from charge distribution in the (1120) plane which include Cr and Te. They 
have clarified that the strongest bond between Cr and Te is the bond in the (1120) 
plane. 

They also investigated the effect of volume change on magnetic moment. They 
calculated magnetic moment as a function of volume for the two cases: (i) bond angle 
of Cr-Te-Cr is fixed (i.e., the axial ratio c/a is fixed) and (ii) bond angle of Cr-Te- 
Cr is unfixed (i.e., c/a is variable). The calculated results indicate that magnetic 
moment is sensitive to lattice change along the a-axis where Cr-Cr distance is large, 
and insensible to lattice change along the c-axis where Cr-Cr distance is small. This 
is a situation similar to that of MnAs. 

4.2.5 Optical Properties 

We can discuss optical properties of compounds on the basis of band calculation 
including the spin-orbit interaction. 

Kulatov and co-workers calculated ferromagnetic bands of MnAs and MnSb by 
using the LMTO method including the spin-orbit interaction [21]. In the calculations, 
they used the muffin-tin approximation and the expression of von Barth and Hedin 
(BH) for the exchange and correlation interactions, and they also took the spin-orbit 
interaction into account. 

Figure 4.15a, b shows the calculated densities of states of MnAs and MnSb, 
respectively. Except the split of the band due to the spin-orbit interaction, overall 



b 




ENERGY, Ryd ENERGY, Ryd 

Fig. 4.15. Densities of states of (a) ferromagnetic MnAs and (b) MnSb. (A) shows total DOS. 
(B) and (C) show partial DOS of Mn-J and As(Sb)-/7, respectively (Reproduced from [21]) 
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features of the bands are similar to the bands calculated by Motizuki and her 
co-workers with the self-consistent APW method. Magnetic moment evaluated here 
by the LMTO method agrees well with experimental value. 

On the basis of these bands, Kulatov investigated optical properties of MnAs and 
MnSb. He obtained a dielectric function by making use of the Kubo formula for the 
linear response theory. Optical absorption coefficient is proportional to the imaginary 
part e2{co) of the dielectric function, and £2{co) consists of two terms: the Drude term 
£2^(0)) arising from intra-band transition, and the (^) term arising from interband 
transition. In the dipole approximation, (^) given by 

^2ico) = V^-27^H m'\i\kX)?fkAl-fkx)mi'-Eki-n(o), (4.1) 

CD 

eh 

where fj^i is the Fermi-Dirac distribution function, 7 = — i — V, and A is a band 

m 

suffix. The matrix elements for the dipole transition {kX'\j\kX) were calculated from 
the LMTO wavef unctions. Interband optical photoconductivity (j(a)) is related to 
£2(^^) as 

e|(a)) = — cj(o)). (4.2) 
ft) 

Figure 4.16a, b shows the photoconductivity a (ft)) calculated for ferromagnetic 
MnAs and MnSb, where magnetic moments are assumed to align along the c-axis 
(z-axis). In these figures, Gxx{(^) and (^^^(ft)) are plotted by the solid and broken 
curves, respectively. Due to anisotropy arising from the hexagonal symmetry of the 
crystal structure, shape and peak positions are different for Oxx and a^z- 

The features of a (ft)) of MnAs and MnSb in Fig. 4.16a, b are similar even in 
the structures. The main peaks of cr(ft)) are attributed to transition of electrons as 
follows: 

• The peak at around /ift) :^ 0.3 eV is attributed to direct interband transition which 
has become allowed owing to mixing between the up- spin and down- spin bands 
due to the spin-orbit interaction. 

• The peaks at around /ift) 1 eV correspond to interband transition between down- 
spin bands of Mn. 

• The broad and strong peak in the energy region of /ift) = 2.5 ~ 4.5 eV is mainly 
due to transition between /?-band of Sb (As) and J-band of Mn. 

In Fig. 4.17, we show XPS spectrum of MnSb measured by Liang and Chen (by 
the dashed curve) [4], and the calculated density of states of ferromagnetic MnSb 
below Fermi level (the solid curve). Since the measurements were carried out with a 
poly-crystal sample, the experimental spectrum corresponds to a^z -\-2axx- There is 
rough correspondence between the experiment and the calculation, however, experi- 
mental data on single-crystal are necessary for detailed comparison. 




^40 



o 



20 



0 f I I I I I I I I I I I I I I I [ 
012345678 
ENERGY, eV 

Fig. 4.16. Optical conductivity for (a) ferromagnetic MnAs and (b) ferromagnetic MnSb cal- 
culated by Kulatov. The solid and dashed curves show Gxx and cr^^, respectively 



4.3 Spin Fluctuations and Anomalous Magnetic 
and Elastic Properties 

4.3.1 Paramagnetic Susceptibility and Anomalous Thermal Effect of MnAs 
and MnAsi_^Px 

With increasing temperature, MnAs shows a first-order transition from ferromag- 
netic state to paramagnetic state at 7c = 3 1 8 K. ^ At this temperature, the structural 
transformation from the Ni As-type to the MnP-type structure also occurs with 2% of 



^ See Sect. 3.2 for experimental detail of MnAs, MnP and their mixed compounds. 
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Energy (eV) 

Fig. 4.17. Density of states for ferromagnetic MnSb (solid curve) and XPS spectrum (dashed 
curve). The dotted curve shows DOS broadened with a Gaussian of full with at half height = 
0.6 eV (Reproduced with permission from [4], Copyright (2009) by the American Physical 
Society) 



volume decrease. In the temperature region above Tq, reverse change to the NiAs- 
type structure takes place at 7^ = 398 K (see Fig. 3.3. The NiAs-type and MnP-type 
structures are also shown in Fig. 4.38). In the temperature region, T > 7^, paramag- 
netic susceptibility obeys approximately the Curie-Weiss law. According to obser- 
vation carried out by Ido et al., the X^-^ curve shows a slight upward convex in this 
temperature range. Magnitude of magnetic moment evaluated from the curve 
in the high-temperature region is 3.8 /Ib- This is larger than saturate moment in the 
ferromagnetic state 3.4 /Ib- 

In the intermediate temperature range (7c <T<Ti) where the crystal structure 
is of the MnP-type, the X^~^ curve shows anomalous behavior that the curve shows 
a broad peak, and crystal volume increases rapidly with increasing temperature (see 
Figs. 3.2 and 3.3, etc.). Crystal volume at T[ is almost equal to the volume just below 
7c. Such anomalous behavior of susceptibility and thermal expansion in the MnP- 
type phase are more remarkable in MnAsi-^Px (0 < x < 0.275). CrSbi-^As^ and 
Cro.4Mno.6As in the MnP-type phase also show such anomalous behavior as shown 
in Part I. According to neutron diffraction measurements by Schwartz [22], magnetic 
moment of MnAs and MnAsi-^Px in the MnP-type phase is smaller than the moment 
in the NiAs-type phase. They also found that magnitude of magnetic moment shows 
smooth change with temperature. 

According to measurements of magnetization process carried out by the group 
of Osaka University (Fig. 3.4 and Fig. 3.15-3.17, etc.), magnetic moment of 
MnAsi_xPx depends little on magnetic fields at temperatures below 7c, whereas 
it depends considerably on magnetic fields at around 7^. In the MnP-type phase in 
the temperature region just above 7c, field-induced metamagnetic magnetization 
curve with hysteresis was observed. As shown in Fig. 3.4, temperature dependences 
of magnetization measured in 20 and 400 kOe revealed that magnetic moment in 
20kOe shows a discontinuous change at 7c. This change is considered due to the 
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structural phase transition at 7c. On the other hand, magnetization in 400 kOe con- 
tinuously changes with temperature, which suggests that MnAs remains to be of 
the NiAs-type structure even in the temperature region above Tq. These results 
mentioned earlier indicate that the metamagnetic behavior is accompanied with the 
structural transition from the MnP-type phase (low-fields) to the NiAs-type phase 
(high-fields). 

Yamaguchi et al. carried out neutron diffraction measurements for the ferromag- 
netic state of MnAs [3]. They reported that a magnetic moment 0.23/iB exists at an 
As- site and it is antiparallel to magnetic moment at Mn-site. Their neutron diffrac- 
tion measurements also revealed that an Sb atom in MnSb has moment 0.3/iB in the 
opposite direction of the moment of Mn atom. 



4.3.2 Spin Fluctuation and Magnetism 

In the itinerant model, correlation between electrons is a very important factor to 
discuss magnetic properties in finite temperature. Moriya and his co-workers de- 
veloped the theory in which the correlation is treated in terms of fluctuation of spin 
density. They succeeded in the unified theory in describing both the extreme case 
of weak ferromagnetism, where the spin fluctuation has a small amplitude and local 
character in the ^-space, and another extreme the case of local moment, where all 
^-modes of spin fluctuation and their mode-mode coupling are important [23-25]. 

Motizuki and Katoh investigated magnetism of MnAs [26, 27] by applying the 
spin-fluctuation theory within the formalism of Moriya-Takahashi [24] and Usami- 
Moriya [25]. 

We start from the following single-band Hubbard Hamiltonian: 



0^ 0^ I 0^ 



U 9 9 



(4.3) 



where the symbol tji denotes a transfer integral, a a creation operator which creates 
an electron with spin a at a site j. nj and Sj in Jifi are operators of charge density and 
spin density, respectively. U and / are an effective intra-atomic Coulomb integral 
and an effective intra-atomic exchange integral, respectively. Assuming that five d 
orbitals degenerate, U and / are related to and / (Coulomb and exchange integrals 
for a single orbital Hubbard Hamiltonian, respectively) as 

U=^{9U'-4f) 
/=^(f/' + 4/). 

By using the Stratonovich-Hubbard transformation, we can transform this many- 
body problem into a single-particle problem under a magnetic field (^) and an elec- 
tric field (t]) fluctuating temporally and spatially. Within the static approximation, 
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a partition function Z (or free energy F) can be written in the functional integral 
form as 



Z = e-^^ = Tr 



-/3(^o-M^) 



and 



(4.4) 



(4.5) 



Here, jU is the chemical potential, A' a total electron number operator, and ^ = \/k^T. 
The functional ^[<^ , T)] = , »?] + ^^i , is given by 



and 



TVexp 



where 



(4.6) 



(4.7) 



(4.8) 



and and T]^ denote local magnetic and electric fields, respectively, operating at the 
jth site. In addition, we introduce a quantity 



1 f 



(4.9) 



where E^q is a Fourier component of the fluctuating magnetic field and Nq is the 
number of lattice points in the system. We note that means that the term ^ = 0 is 
excluded. Xa in (4.9) is a square average of local amplitude of the -field, and related 
to a square average of local amplitude of spin fluctuation as 



(4.10) 



By assuming the mode-mode coupling of spin fluctuation to be local and applying 
a saddle-point approximation for the charge field, Moriya and the co-workers intro- 
duced the following form of the functional: 



?7rT ' 

P q a 



(4.11) 
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The first term of (4.11) represents the nonlocal part, where terms independent 
on relative direction of spin fluctuation are neglected. An appropriate function of q is 
assumed for Xqa [x, ^q] . The second term of (4. 1 1) represents the local part, for which 
Usami and Moriya obtained an expression for L[x, (^o] by using the coherent potential 
approximation (CPA). 

From the functional thus obtained, the partition function of the system is 
calculated and then various quantities such as paramagnetic susceptibility and local 
magnetic moment can be evaluated successfully. 

We show here the results of calculations carried out for paramagnetic MnAs and 
MnAsi_xPx [26, 27]. First, we simplified the density of states calculated by the self- 
consistent APW method by replacing the curves with a set of straight lines and then 
obtained model density of states as shown in the lower part of Fig. 4.18. We use the 
model density of states for the present calculations. For example, the curve denoted 



by (1) in Figs. 4.18 and 4.19 show calculated local magnetic moment J (S^-) and 



inverse paramagnetic susceptibility x~\ respectively. They are calculated with the 
model density of states (1) and in the case that electron number per orbital n is 
0.65 and U and / are equal to 2.8 (in unit of 5/6 eV). As shown in Fig. 4.19 the 
calculated X^'^ curves show the Curie- Weiss behavior with weak upward convex. 
Magnetic moment estimated from the curve is 2.58/iB, which is slightly larger 

than the value 2.30 jUb evaluated from the model density of state and the Hartree- 
Fock approximation at OK. As seen from Figs. 4.19 and 4.20, there is thus good 
correspondence between the present calculations and the experimental results except 
for the numerical values. For example. Curie temperature estimated from an equation 
= 0 is only 1/10 of the observed value. It is needed to point out that the relation 
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Fig. 4.18. Calculated local magnetic moments mioc using the model density of states (l)-(4). 
The model density of states (1) corresponds to the NiAs-type MnAs, (2)-(4) correspond to 
MnAs distorted to the MnP-type structure (Reproduced with permission from [26]) 
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Fig. 4.19. Inverse paramagnetic susceptibility calculated for the model densities of states (1)- 
(4). The dotted and dash-dotted curves represent the intermediate temperature range 

Tc<T <Tt with Tt = ITq and Tt = STc, respectively (Reproduced with permission from [26]) 
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Fig. 4.20. Temperature dependence of (Ap — Ap^) calculated for the model densities of states 
(l)-(4) shown in Fig. 4.18. The dotted and dash-dotted curves represent (Ap — Ap^) in the 
intermediate temperature range Tc<T <T[ with T[ = ITq and T[ = 3Tc, respectively (Repro- 
duced with permission from [27]) 
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X~ gives 7c for the case of a second-order transition and the transition at Tq 
of MnAs is of a first order. We may expect a better estimation of Tq by considering 
effect of degeneracy of orbitals and quantum effects which are neglected in the spin- 
fluctuation theory. 

Next, we investigate paramagnetic susceptibility of MnAs in the intermediate 
temperature range of 7c < T < 7^ where the crystal structure is of the MnP-type. 
As noted in Sect. 4.2, the APW band calculations have revealed that the height of 
the sharp peak of density of state near the Fermi level in the NiAs-type structure 
becomes smaller with lattice distortion to the MnP-type structure. In order to take 
such change of band structure into account, we modified the model density of states. 
The model density of states in the MnP-type structure is shown in the lower part 
of Fig. 4. 18 for the three cases denoted by (2), (3), and (4), which correspond to 
the cases where p{Ey^) is 90%, 85%, and 80% of p{Ey^) of (1) in the NiAs-type 
phase, respectively. The solid curves denoted by (2), (3), and (4) in Figs. 4.18 and 
4.19 display the temperature variations of the local magnetic moment and inverse 
paramagnetic susceptibility calculated for the model densities of states (2), (3), and 
(4), respectively. It is seen from Figs. 4.18 and 4.19 that the local magnetic moment 
decreases in the order (1), (2), (3), and (4), resulting in steeper inclination of the 
X^-T curves. In order to explain the observed X^'^ curve of MnAs, in Fig. 4.19 
we start from a point on the curve (1) in the high-temperature region and consider the 
X~^-T curve in the process of lowering temperature. With decreasing temperature, 
X~^ decreases first along the curve (1), and in the MnP-phase below the transition 
temperature Tt, x^ changes to the values of points on the curves (2), (3), and (4) 
successively because the lattice distortion increases as temperature decreases. The 
dotted- and the dashed-curves show temperature dependences of for two cases 
of (i) Tt = 27c and ii) Tt = 37c, respectively. For the case (i), we assumed that the 
lattice distortion changes monotonically with temperature in the temperature range 
between Tt and 7c. For the case (ii), the monotonical change of distortion in the 
temperature region between 7^ and 1.5 7c, and constant distortion in the temperature 
region between 1.5 Tq and Tq are assumed. Measurements of lattice distortion for 
MnAs and MnAso.8Po.2 has revealed that MnAso.8Po.2 takes the MnP-type structure 
in wider temperature range than MnAs (see, for example. Fig. 3.19 and [28]). It also 
has been found that temperature dependence of the lattice distortion of MnAso.8Po.2 
is large only in a narrow temperature region just below 7^, but almost constant in 
other temperature regions below 7^ [28]. Thus, the X^'^ curve calculated for the 
case (i) corresponds to MnAs and that for the case (ii) to MnAso.8Po.2- The X^'^ 
curves of these compounds show a peak in the temperature region of the MnP-type 
phase. The peak of MnAso.8Po.2 is especially remarkable (see, for example. Figs. 3.2 
and 3.13). 

As mentioned earlier, we have shown that the anomalous temperature depen- 
dences of susceptibility observed for MnAs and MnAsi-^Px compounds are rea- 
sonably understood by the spin- fluctuation theory including the effect of lattice 
distortion of the MnP-type structure. 
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4.3.3 Spin Fluctuation and Elastic Properties 

We describe in this section how to explain anomalous elastic properties of MnAs and 
MnAsi_;cPx in the MnP-type phase from the viewpoint of itinerant electrons taking 
the effect of spin fluctuation into account [27] . 

dF 

Cohesive pressure of electron system is given hy P = — — with the free en- 

dV 

ergy F. For simplicity we assume that electronic energy is proportional to band 
width W and that only W depends on volume. The pressure P in this approximation 
is written as 



P = ?iA/V. 
V dW 

In this equation, A is defined by A = r — and A is given by 

W oV 

A V 



(4.12) 



(4.13) 



From free energy evaluated by using the spin- fluctuation theory described in 
Sect. 4.3.2, we calculated the quantity A for paramagnetic NiAs-phase and para- 
magnetic MnP-phase. The free energy per atomic site is written as a summation 
of Fq and AF, which are terms irrelevant and relevant, respectively, to the spin 
fluctuation. The quantities Fq and AF are given by 



^0^ 



1 



and 



^F 



-L(x), 



(4.14) 



(4.15) 



where the nonlocal term has been omitted. By using an expression of the local term 
L{x) obtained within the coherent potential approximation (CPA), we can calculated 
the quantity A from (4.13). The parameter Ap, which means A of the paramagnetic 
state, is given by 

Ap = -^ y^de/(e + ju) xIm{(e + jU-i:p-7p)Fp}, (4.16) 

where / is the Fermi-Dirac distribution function and e energy of electrons, and other 
symbols in (4.16) are explained next. We note that | T]^ p and x in the expression of AF 
are assumed not to depend on explicitly. The chemical potential /i is determined by 



- / d£/(£ + ja)ImFp: 
n J 



(4.17) 



4.3 Spin Fluctuations and Anomalous Magnetic and Elastic Properties 



101 



where n is electron number. The symbol 7p is nU /2, is a diagonal element of 
single-particle Green's function, and Xp denotes site-independent self-energy. Within 
the Hartree-Fock approximation in which the spin fluctuation is neglected, the local 
magnetic moment is zero and Xp = 0. In this case, Ap becomes Ap ^ given by 

Af = 2|d£/(£ + 7^)£p(£), (4.18) 

where p{e) is the density of states. Contribution of the spin fluctuation to A is given 
by (Ap — Ap^). Change of pressure is related by the compressibility k to change 
of volume AV as 

AV 

— = K^p (4.19) 

We see from (4.12) that AP is written by A A, therefore volume change in (4.19) is 
given by 



^ = A A(Ap -Af)/y + A K^Af/y 
AVi ^V2 

= — ^ + — . (4.20) 

The first and the second terms denote volume changes arising from the spin fluctua- 
tion and that from change of electronic energy, respectively. 

Calculated (Ap — Ap ^) is plotted against temperature in Fig. 4.20. The curves 
(l)-(4) are results calculated for the model densities of states (l)-(4) shown in 
Fig. 4.18, respectively. Since (Ap — Ap^) arises from the spin fluctuation, temper- 
ature dependence of this quantity reflects temperature dependence of local ampli- 
tude of spin fluctuation plotted in Fig. 4.18. As seen in Fig. 4.20 that (Ap-A^F) 
at each temperature decreases in the order (1), (2), (3), and (4), which means that 
spin fluctuation is suppressed as the lattice distortion to the MnP-type increases. To 
explain anomalous thermal expansion observed for MnAs, etc., (see Figs. 3.3, 3.19, 
etc., in Part I) we consider temperature dependence of (Ap — Ap^) in the process of 
decreasing temperature; (Ap — Ap^) decreases first along the curve (1), and in the 
temperature region of 7c < T < 7^ below 7^, it changes successively to the points 
on the curves (1) ^(2) ^(3) ^(4) as the lattice distortion increases. The dotted- 
and dashed-dotted-curves denote temperature dependence of (Ap — Ap^) for the case 
Ti = Tq and 7^ = 37c, respectively. These results clearly indicate that the spin fluctu- 
ation plays an important role on the remarkable volume change with temperature for 
MnP-type MnAs and MnAsi-^Px, because temperature dependence of (Ap — Ap^) 
gives rise to the volume change. 

Next, we investigate the volume change at 7c. Using values /c = 4.5 x 
lO-i^m^/N, y = 34 X lO-^^m^ and A = 5/3 (estimated by Heine [29]), we eval- 
uated at 7c the quantity AVi/V from (Ap - A^^) to be 3.4%, 4.8%, and 5.8% for 
the changes of (Ap — Ap^) between curves (1) and (2), (1) and (3), and (1) and (4), 
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respectively. We also estimated AV2/y comparable to AVi/V from the tight-binding 
calculation for the MnP-type structure, which means that decrease of local magnetic 
moment due to the structural transition and energy change of electron system have 
comparable contributions to the discontinuous volume change at Tq in MnAs (see 
Fig. 3.3 in Part I). Calculated values of AV/V is 7.1% for the change (1) ^(2), 9.5% 
for (1) ^(3), and 12% for (1) ^(4), respectively. Experimental volume change at 
7c is 2% for MnAs and 8 - 10% for MnAsi-;,?;^. 

Einally, we note that we estimated A 1.1 from our band calculation for MnAs. 
If we use this value of A instead of Heine's result A = 5/3, the calculated values of 
AV/y become smaller by about 2/3 than the value mentioned earlier. 

4.3.4 Paramagnetic Susceptibility of CoAs and FeAs 

The compounds EeAs and NiAs crystallize in the MnP-type and the NiAs-type 
structures, respectively, while CoAs shows the structural phase transition at 1,250 K 
from the MnP-type to the NiAs-type structures as temperature rises. Below this tem- 
perature, the crystal structure of CoAs becomes the MnP-type. EeAs takes the MnP- 
type structure at all temperatures, and NiAs remains in the NiAs-type structure ex- 
cept at very temperatures. 

Morifuji and Motizuki calculated paramagnetic susceptibility x of these com- 
pounds by using the spin- fluctuation theory and nonmagnetic bands calculated by 
the APW method [30]. It has been revealed that % of NiAs is almost temperature- 
independent, while x'^ of CoAs and EeAs, which have the MnP-type structure, 
shows anomalous temperature dependence. According to the experimental data [31], 
inverse susceptibilities of CoAs and EeAs obey the Curie-Weiss law in high temper- 
ature region and increase for further decreasing of temperature. 

The temperature dependences of paramagnetic susceptibility of CoAs and EeAs 
are similar to that of a weak ferromagnetic semiconductor EeSi. Takahashi and 
Moriya explained anomalous temperature dependence of susceptibility of EeSi on 
the basis of the spin- fluctuation theory and characteristic shape of density of states 
[32]. According to their theory, negative mode-mode coupling is thermally induced 
when the Eermi level lies at a dip of density of state and the density of states just 
below the Eermi level has a large value. In this situation, the negative mode-mode 
coupling gives rise to a dip of the curve and a rapid increase of ;^~Mn the tem- 

perature region below the dip. Eurthermore, the Curie-Weiss behavior of susceptibil- 
ity occurs in the temperature region above temperature of the dip. The Curie-Weiss 
behavior is brought about by the saturation of local amplitude of spin fluctuation. 

Morifuji and Motizuki calculated X^~^ curves for CoAs and EeAs by taking 
account of the spin-fluctuation effect, as well as the fact that the densities of states 
for EeAs and CoAs shown in Eig. 4.9 are similar to that of EeSi [33, 34]. 

We show in Eig. 4.21 the calculated results for CoAs. We used the model density 
of states, shown in the inset, in which the calculated curve is replaced by fragments of 
straight lines. We see that the temperature dependence of inverse susceptibility x~^ 
agrees qualitatively with that of measurement plotted by x [31]. This result indicates 
that the anomalous temperature dependence of paramagnetic susceptibility of 
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Fig. 4.21. Inverse paramagnetic susceptibility of Co As calculated by using the model density 
of states inserted in the figure is plotted by the solid curve. The observed susceptibilities at 
several temperatures are shown by x 



Co As is well understood in terms of the spin- fluctuation effect except discrepancy in 
magnitude of and the temperature of the dip. We have obtained similar results 
for FeAs. 



4.4 Fermi Surface of NiAs and the de Haas-van Alfen Effect 

We can investigate Fermi surface by observing the de Haas-van Alfen effect. Since 
the cross section normal to an applied magnetic field is related to period of the 
de Haas-van Alfen oscillation, it is possible to investigate shape of Fermi surfaces 
from the de Haas-van Alfen measurements carried out in various directions of ap- 
plied magnetic fields. 

Recently, various techniques and theories for band calculation have also been 
developed to calculate more accurate band structures. As for the correlation between 
electrons, the local density approximation, which is usually used, may be insufficient 
because average electron density is used. Recently, improvement of the LDA such 
as GGA has been carried out. It is thus important to justify calculated bands by 
examining the shape of Fermi surfaces. 

The Group of Kamimura measured the de Haas-van Alfen oscillation using sin- 
gle crystals of NiAs and NiSb with the NiAs-type structure [35]. Since the NiAs-type 
structure has two metal atoms and two anions in the unit cell, the dispersion curves 
near the Fermi level have a complicated dependence on wavevector due to mixing 
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between J-orbitals of Ni and /7-orbitals of As or Sb. Many branches are observed in 
the experiments, and therefore it is difficult to reproduce the shape of Fermi surfaces 
from de Haas-van Alfen measurements. 

Recently, however, Harima developed a computer program to calculate de Haas- 
van Alfen frequencies from Fermi surfaces obtained by band calculations [36]. It 
has enabled us to easily make comparison between calculated band structures and 
de Haas-van Alfen measurements. To make direct comparison, we first constructed 
Fermi surfaces of NiAs from results of band calculation and determined dependences 
of the de Haas-van Alfen frequencies in the direction of magnetic field, and then 
compared them with the experimental data [35, 37]. 

We calculated initially band structure of nonmagnetic NiAs and NiSb by us- 
ing the APW method, and later the FLAPW method to obtain more accurate band 
structures. We applied the LDA with the GL formula for the exchange and corre- 
lation interaction. Scalar relativistic effects were taken into account and the spin- 
orbit interaction was omitted in the calculations. The program codes TSPACE and 
KANSAI-92 were used for these calculations. 

Dispersion curves of NiAs and NiSb are shown in Fig. 4.22a, b, respectively. The 
bands in the low-energy region consist of As-4^ or Sb-5^ orbitals. The 16 bands 
above the gap are p-d mixing bands consisting of Ni-3d and As-4p (or Sb-S/?) 




Fig. 4.22. Dispersion curves of nonmagnetic band of (a) NiAs and (b) NiSb (Fig. 4.22a is 
reproduced form [35]) 
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Fig. 4.23. Density of states of nonmagnetic band of (a) NiAs and (b) NiSb (Fig. 4.23a is 
reproduced form [35]) 



orbitals. We show the total densities of states of the p-d mixing bands p{E) with 
the components of Ni- J orbitals and As- (or Sb-) p orbitals for NiAs and NiSb in 
Fig. 4.23a, b, respectively. The density of states in both compounds can be divided 
into three main parts: p-d bonding band, nonbonding band consisting mainly of d or- 
bitals in the middle region of energy, and p-d antibonding band in the highest region 
of energy. The gross feature of the density of states is similar to that of MnAs. Total 
width of the p-d mixing band is larger than that of MnAs shown in Figs. 4.2 and 4.3. 
The Fermi level is an energy far above the energy region with the large peak of the 
density of states, p(£f) = 29.12 (states/Ry unit cell) for NiAs and 14.53 (states/Ry 
unit cell) for NiSb are obtained. These values are much smaller than 156 (states/Ry 
unit cell) of MnAs, which is consistent with the fact that MnAs is a ferromagnet but 
NiAs and NiSb are paramagnetic. 

The value of p (£^f) of NiAs evaluated from measurement of temperature depen- 
dence of specific heat is 34.54 (states/Ry unit cell) [37], which is somewhat larger 
than the calculated value of 29.12, however, the difference between experimental 
and calculated values in NiAs is smaller than that in MnAs. This fact suggests that 
electron-electron interaction and electron-phonon interaction in NiAs is weaker than 
those in MnAs. Width of the density of states below the Fermi level for both NiAs 
and NiSb is about 6 eV (c^0.44Ry) and there is a sharp peak at 2 eV (c^0.15 Ry) en- 
ergy lower than the Fermi level. Such character of the density of states corresponds 
well to the results of photoemission measurements [38]. 

Figure 4.24 shows calculated Fermi surfaces of NiAs. NiAs has three Fermi sur- 
faces: two hole surfaces with cylindrical shape around the FA-axis and an electron 
surface whose shape is somewhat complicated. 




The de Haas-van Alphen frequencies calculated from the Fermi surfaces are plot- 
ted as a function of angle of direction of magnetic field in Fig. 4.25. The branches 
labeled e, m, and d in Fig. 4.25 seem to correspond generally well to the experimental 
branches labeled £, /i, and 5 as shown in Fig. 4.26 [35]. 

However, the a and j3 branches in Fig. 4.26 cannot be explained by the present 
calculations, and the c- and ^-branches which are seen in the theoretical curves are 
not measured. The absence of observations of the c- and /i-branches may be ex- 
plained by a possible structural transformation due to nesting effect between the 
hole surface around the FA-axis and the electron surface around the KH-axis. Recent 
X-ray diffraction measurements have revealed that NiAs shows a structural transfor- 
mation at a very low temperature [39]. Due to this structural transformation, it is 
considered that the hole surface around the FA-axis and the electron surface around 
the KH-axis disappear. 

Figures 4.27 and 4.28 show the Fermi surfaces and the calculated de Haas- 
van Alfen frequency of NiSb, respectively [37,40]. The results of measurements 
are shown in Fig. 4.29. As seen in NiAs, there is good correspondence between the 
e branch in Fig. 4.28 and the £ branch in Fig. 4.29. Unlike the case of NiAs, the 
a- and the j3-branches which correspond to the c\- and the /ii -branches in Fig. 4.28 
are observed in NiSb. This result corresponds well to the fact that NiSb does not 
show structural transformation, but remains in the NiAs-type structure even in low 
temperature. 
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Fig. 4.25. Calculated de Haas-van Alfen frequencies of nonmagnetic NiAs (Reproduced 
form [35]) 



4.5 Pressure Effect on Magnetic State of CrTe, CrSe, and CrS 

CrTe becomes a ferromagnet below Tq = 340 K. By application of high pressure, 
Tq decreases remarkably and the magnetic moment also decreases, which suggests a 
kind of pressure-induced phase transition. 

To clarify the effects of high pressure on magnetic properties of CrTe, we carried 
out band calculations for nonmagnetic, ferromagnetic, and antiferromagnetic states 
for various values of the unit cell volume. In these calculations, the axial ratio c/a 
is fixed to the value at room temperature. For the antiferromagnetic calculation, the 
sublattice magnetic moments are assumed to be in the c-plane and are antiparallel 



108 4 Electronic Band Structure and Magnetism of Ni As-Type Compounds 



6000 
5000 
4000 

^ 3000 
£)2000 

I 1000 



TA 



200 



TK TM 

I I I I I I I I I M I I I I I I 



TA 



NiAs 



' 9 

u»« £2,3 • ; ; £i,3 £2 



' ® 9 9 9 - 

9 ®«9 9 ®® 



[0001] 30 60 [1210] [1100] 60 30 [ooOl] 
Field angle (Degrees) 



Fig. 4.26. Experimental de Haas-van Alfen frequencies of NiAs (Reproduced form [35]) 
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Fig. 4.28. Calculated de Haas-van Alfen frequencies of NiSb 



between adjacent c-planes. We employed the APW method for the first calculations, 
later we used the FLAPW method to obtain more accurate band structures [41-43]. 
In this section, we show the results obtained by the FLAPW method. The density of 
states of ferromagnetic CrTe calculated by the FLAPW method is shown in Fig. 4.14. 

In Fig. 4.30, we show the calculated total energy as a function of the lattice con- 
stant a for nonmagnetic, ferromagnetic, and antiferromagnetic states in CrTe [41]. 
The total energy for the ferromagnetic state takes the minimum value at a = 3.970 A, 
which agrees well with the experimental value a = 3.981 A. We can express the total 
energy curve of ferromagnetic state by a function 

£'Totai = ^0 + - aof ^b3{a- gq)^ -^b^{a- a^f , (4.21) 

where Eq = -31357.17 (Ry/unit cell), gq = 3.970, Z?2 = 0.471, /?3 = -0.575, and 
/?4 = 1.686. 
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Fig. 4.30. Total energy curves for ferromagnetic, antiferromagnetic, and nonmagnetic CrTe as 
a function of lattice constant a (Reproduced from [41]) 
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Since the APW method gives the minimum energy at the lattice constant a = 
4. 180 A, we see that the FLAPW method brings about improvement of calculated 
results. Generally, a full potential calculation reproduces more correctly the exper- 
imental values of lattice constant than the APW method. As shown in Fig. 4.30, the 
ferromagnetic state is most stable at a = 3.981 A, which is consistent with the fact 
that CrTe is a ferromagnet under atmospheric pressure. 

It is seen in Fig. 4.30 that the antiferromagneic state has smaller energy than that 
of the ferromagnetic state for the lattice constant smaller than a = 3.58 A, which 
indicates that pressure-induced phase transition from the ferromagnetic state to the 
antiferromagnetic state can occur. The inset of Fig. 4.31 shows relation between cal- 
culated electronic pressure and the lattice constant a. From this curve, we see that the 
pressure-induced magnetic transition may occur at a critical pressure of about 40 GP 
which corresponds to the lattice constant 3.58 A. The APW calculation gives critical 
pressure 20 GP. In Fig. 4.31 total magnetic moment of Cr and Te in the ferromagnetic 
CrTe are plotted against the lattice constant a. Large part of the magnetic moment 
arises from Cr atom, however, Te atom has a small negative magnetic moment in the 
opposite direction, (see also Sect. 4.2.4). 

The total magnetic moment of CrTe is calculated to be 3. I/Ib /formula at a = 
3.970 A. This value is not so different from observed moment 2.29/iB/formula. The 
APW calculation gives magnetic moment 3.9/iB/formula. 

We can evaluate electron pressure P from the equation P = — d£^Totai/dV. Since 
the unit cell volume is given by V = 1.35a^, the pressure P is expressed as P = 
— 0.247(2~^d£'Totai/d<3. The bulk modulus, which is given by B = —dP/dlogV 
(V = Veq), is calculated to be 5 = 404kbar, or the compressibility (B~^) 2 A3 




Fig. 4.31. Magnetic moments in ferromagnetic CrTe plotted against the lattice constant a. 
Inset: Electronic pressure as a function of lattice constant (Reproduced from [41]) 
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Fig. 4.32. Total energy for ferromagnetic, antiferromagnetic, and nonmagnetic CrSe (Repro- 
duced from [41]) 

X 10~^ kbar~^ . The calculated compressibility evaluated by the FLAPW calculation 
is slightly smaller than the value 2.6 x 10 ~^ kbar~^ evaluated by the APW calcula- 
tion. Pressure coefficient of magnetic moment is calculated as dM/dP — — 1.33x 
10~^j(iB kbar~^ at (2 = 3.981 A, and then the magnetic moment decreases at the rate 
k = —dlnM/dP = —4.21 x 10~^kbar~\ which is small but in the similar order in 
comparison with the observed value — 12 x 10~^ kbar~^ [42]. 

We also carried out the total energy calculations for CrSe and CrS. As an exam- 
ple, we show the relation between £^Totai and the a-axis of CrSe in Fig. 4.32. The total 
energy calculations show that the antiferromagnetic state is most stable even under 
high pressure both for CrSe and CrS. 



4.6 Magnetic Ordering and Instability of Paramagnetic State 
4.6.1 Double-Helical Magnetic Ordering of MnP-Type Compounds 

The total energy calculation, as described in Sect. 4.5, makes it possible to clarify 
which magnetic state occurs in a material with itinerant electrons. It is, however, a 
hard task to find out the most stable magnetic state from the calculations of total en- 
ergies, especially in a material with a lot of magnetic atoms in a unit cell because we 
need to carry out total energy calculations for many possible magnetic orderings. In 
addition, some of transition metal compounds such as CrAs explained in Part I have 
helical magnetic ordering. If wavevector of the helical ordering and lattice vectors 
are incommensurate, it is impossible to carry out band calculation. 
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An alternative way to search for a possible magnetic state is to study instability of 
paramagnetic state. In this section, we describe how we can investigate the instability 
of paramagnetic state on the basis of paramagnetic band calculations as well as the 
Hubbard model [44]. 

The transition metal pnictides such as FeAs and CrAs with the MnP-type crys- 
tal structure show the double-helical magnetic ordering. Magnetic transition tem- 
perature is 77 K for FeAs and 250 K for CrAs. We describe later the procedure to 
investigate instability of the paramagnetic state for the case of FeAs. 

As shown in Fig. 4.33, the unit cell of FeAs (orthorhombic MnP-type structure) 
contains four Fe atoms and four As atoms. Paying attention only to Fe atoms and 
neglecting the difference in As atoms surrounding the respective Fe atoms, we can 
take a reduced unit cell which contains only two Fe atoms. The thick lines in Fig. 4.33 
shows the reduced unit cell containing two Fe atoms labeled as "a" and "b". 

Experimental studies have revealed that the helical ordering of FeAs is described 
by a wavevector Q = (0, 0, 0.375 x 2k /c) (See also Tables 2.4 and 2.6). We note that 
Q of 0.375 xln/c corresponds to wavelength 2.67b. We also note that the lattice 
vector c here corresponds to b given in Table 2.6. It also has been known that the 
phase difference between magnetic moments at the a- site and the b-site in the same 
unit cell is (j) = 140°,^ and magnitude of the magnetic moment of the Fe atom is 
0.5mb [31]. 




Fig. 4.33. The magnetic ordering of FeAs with the MnP-type structure. The thick lines denote 
the reduced unit cell (see the text), ti and t2 denote transfer integrals between Fe atoms at 
different sublattices. is a transfer integral between Fe atoms at the same sublattices. See 
also Fig. 2. 1 and Table 2.6 



^ (j) is denoted by Ou in Part I. 
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4.6.2 Instability of Paramagnetic State 



We start from a Hubbard Hamiltonian with transfer integrals and intra-atomic 
Coulomb interaction. We assume that each Fe atom has a single orbital. For the 
transfer integral, we consider three types of interaction: between the same sublat- 
tice (a-a or b-b) and t\ and t2 between the different sublattices (a-b) as shown in 
Fig. 4.33. 

By Fourier transformation and the Hartree-Fock approximation for the Coulomb 
interaction, we have the following expression of the Hubbard Hamiltonian: 

^ = S ^1 (^WkG^kG + 7^1* {k)bl^aka + T3 (k) (al^aka + bl^ha) 
ko 



X X V'i-oal+q^ako + Bg-abl+g^bko] -NU^ [Ag+Ag- + Bg+Bg-] 

kq O q 



[^9 +^9 ^1+9-^*+ +^-qVk+q-+B+gbl^bk+q- 

kq 



-NUj^[\A, 

9 



\B 



+ |21 



(4.22) 



where 



and 



(4.23) 



j j 



In (4.24), Ti denotes position of an a- atom in the /th unit cell, and rj + T position 
of a b-atom in the jth unit cell. 

Using this Hamiltonian, we investigate instability of paramagnetic state at 0 K 
against an external magnetic field modulated by a wavevector q. In order to do it, 
first we add Zeeman terms to the Hamiltonian. Next, by solving equations of mo- 
tion of the quantities Aqa, etc., we have the following relation between ^-dependent 
magnetic fields and magnetic moments: 
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rH''{q)\ ^ X"''i<l)\ fM\q)\ ^ /M-{q)\ 

yH\q)j ~ U'n?) X^iq)) \M\q)) ~ \M\q)) ' 



In this equation, M^{q) = jlsiAq- — A^+) and M^{q) = jlB{Bq- — Bq^) are Fourier 
components of spin density at the a- and b-sites, respectively. H^{q) and H^{q) 
denote ^-dependent magnetic fields operating on M^{q) and M^{q), respectively. 
X^^{q), etc., are elements of the inverse magnetic susceptibility tensor. These quan- 
tities are given by 

X-iq) = f'iq) = -^m±![llM^mm (4 26a) 



r''{q)=x'"'*{q) 



1 r2{q)* 
2^imY-\r2{q)r 



where 



(4.26b) 



(4.27a) 



and 



1 ^ T,%k^q)n{k) 
-4Nf\n{k^q)\\n{k)\ 



(4.27b) 



In these expressions, a and j3 denote bands, and XQ^{k^q), etc., are bare electronic 
susceptibilities which are given by 



Xq 



where Ej^ denotes the band energy obtained by diagonalizing the Hamiltonian. When 
paramagnetic state is unstable against an infinitesimal ^-dependent magnetic field, a 
relation 



d^mq)]=r%qy-\r{q)\=o 



(4.28) 



must be satisfied. Substituting (4.26a)-(4.27b) into (4.28), we can evaluate a critical 
value of the Coulomb interaction Uc as a function of the wavevector q.lf U >Uq, 
paramagnetic state is unstable, resulting in occurrence of a magnetic order. When Uq 
takes the minimum value for a wavevector q = Q,we may say that a magnetic order 
described by the wavevector Q is most likely to realize. 

Substituting H''{q) = H^{q) = 0 into (4.25), we have a relation jf^^(^)M^(^) + 
X''^{q)M^{q) = 0. When ^ = g, we see that \x''''{q) \ is equal to Ix^'^iq)]. Therefore, 
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absolute value of M^{Q) is equal to the absolute value of M^{Q). In the real space, 
magnetic moments at the a- and the b-sites are written as M^(r) = M^{Q) and 
M^(r) = M^(e)e-^^•^ respectively. 

From these expressions, the phase difference (j) between the moments of the a-site 
and the b-site is defined by the following equation: 

M%Ll)=M!Me-ie-^.e'^ (4.29) 

Since X^^{Q) is a real number and the phase factor of X^^{Q) is derived from 
12(2)*, (j) the phase difference between moments at a-site and adjacent b-site is 
given by 

(/) = TT + arg [r2{Q)] + arg [e-^^'^j . (4.30) 

We carried out calculations using parameters ti = — O.OTRy, t2/\ti\ = —1.0, 
I = —0.2, and electron number n = 1.1. These values were chosen so that they 
reproduce as nicely as possible the features of Fermi surfaces and density of states 
of FeAs which were obtained by the APW band calculation. 

The calculated results are shown in Figs. 4.34 and 4.35. As shown in Fig. 4.34, 
Uc takes the minimum value for Q = 0.4 x2k/c, which means that the paramagnetic 
state is most unstable against magnetic ordering modulated with this wavevector. We 
see that this value of Q agrees well with the experimental value Q = 0.375 x 27r/c. 
Furthermore, as shown in Fig. 4.35, the phase difference (j) takes 158° at2 = 0.4x 
271 /c. This value of (j) is also close to the experimental value (j) = 140°. 




Fig. 4.34. Critical value of Coulomb interaction Uc evaluated as a function of wavevector q. 
As shown by the arrow, Uc takes the minimum for Q = 0.4 x2k/c 
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Fig. 4.35. Phase difference (j) between spins at the a- site and b-site calculated as a function of 
q. For Q = 0.4 x2n/c, (j) = 158°. This value of (j) is close to the observed value (j) = 140° 



4.6.3 Energy of Double-Helical Spin Density Wave State 

Next, we investigate the double-helical spin density wave (DHSDW) state and eval- 
uate electronic energy with the Hubbard Hamiltonian used in the instability study. In 
order to describe a DHSDW state modulated by a wavevector Q, we neglect ^-terms 
except the one q = Qm (4.22), and have 



je = J^^\k,Q)je{k,Qmk,Q)- 
k 



NUn^ 



-NUl\A+\^ + \B+^^ 



). (4.31) 



The J^{k, Q) and ^(fc, Q) in this equation are expressed by the following equations: 







-UA 


-Q 0 


\ 


T*ik) 


Uk) 
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-UB 


-Q 


-UA+ 


0 


mk+ 


Q) Ti{k + 


Q) 


\ 0 




T*{k^ 


-Q) Uk + 


Q)J 



(4.32) 
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and 



^k+Q- 



(4.33) 



Total energy of the DHSDW state is given by 



£'SDW = X X ^kiif{Ekii) 
k 11=1 



NUn^ 



-NU\ 



{K\'^K\')^ (4.34) 



where E^^ is eigen energies evaluated by diagonalizing the Hamiltonian, f{Ek^) the 
Fermi-Dirac distribution function, Aq and Bq Fourier component of spin density at 
the a- and b-atoms, respectively. The phase difference (/) between spins at the a-site 
and the b-site is given by 



5e-^^^ = e^^ 



(4.35) 



We obtained Aq and Bq by self-consistent calculations, and evaluated the total en- 
ergy £^SDW at 0 K as a function of Q and U . 

In Fig. 4.36, £'sdw is plotted as a function of wavevector Q for the two values 
ofU /\t\\. The total energy £'sdw takes the smallest value at 2 = 0.4 x {2k /c). This 
result corresponds to the result of the instability study described earlier. 




Fig. 4.36. Energy of the DHSDW state calculated as a function of wavevector Q for two values 
ofU /\ti\. Corresponding to the result of the instability study, the energy takes the minimum 
valueat2 = 0.4x2;r/c 



4.7 Phase Transition from the NiAs-Type to the MnP-Type Structure 119 



1 






























-0.2 
-0,1 



Fig. 4.37. Upper part: £^SDW —Epara plotted as a function of U/\ti\. For the curves of each 
value of Q, in the region ofU/\ti \ larger than the value denoted by the arrow, E^jm is smaller 
than £'para- Lower part: \Aq \ plotted as a function of U/\ti \ calculated for Q = 0.4 x2k/c 

We show energy difference between the DHSDW state and paramagnetic state, 
^SDW — Epsim, at 0 K in Fig. 4.37. The vertical arrows indicate the values of U, where 
^SDW is equal to £'para for the respective values of Q shown in the figure. These 
values of U corresponds to Uc shown in (4.28). Since Uc takes the smallest value 
for Q = 0.4 X {2k /c) as shown in Fig. 4.37, the calculated energies of the DHSDW 
state is consistent with the result of instability study. The phase difference (j) between 
spins in the same unit cell is evaluated from (4. 35) as (^ = 153. The value of \Aq \ for 
Q = 0.4 X {271 /c) is shown in lower part of Fig. 4.37 as a function of U. Since the 
magnetic moment per atom is given by 2|A^|, we see that a parameter U/\ti \ =4.8 
reproduces the experimental value of magnetic moment O.S/ie/Fe, therefore we can 
estimate the magnitude of the Coulomb interaction to be [/ = 4.6eV by using ti = 
-0.07 Ry. 

From the discussion described in this section, we see that the double-helical mag- 
netic ordering of FeAs is well explained by the itinerant picture based on the multi- 
band Hubbard model whose parameters (transfer integrals and electron number) are 
obtained so as to reproduce results of the APW band calculations. This theory is 
widely applicable not only to other MnP-type compounds which has the double- 
helical ordering but also to compounds with several magnetic atoms in a unit cell. 

4.7 Phase Transition from the NiAs-Type 
to the MnP-Type Structure 

The structural phase transition between the NiAs-type and the MnP-type struc- 
tures, which is easily induced by temperature change or by composition change of 
compound, is one of the characteristic properties of the transition metal pnictides 
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and chalcogenides with the Ni As-type structure. Correlation between magnetic prop- 
erties and the crystal structure is an interesting problem which has not been fully 
understood yet. In transition metal chalcogenides, VS takes the MnP-type and the 
NiAs-type structures below and above 7^ = 850 K, respectively. Tio.QsSe takes the 
NiAs-type structure, whereas Tii.osSe takes the MnP-type structure at room temper- 
ature. Ti^Vi-xS takes the NiAs-type for x > 0.66 and the MnP-type for x < 0.66. As 
for the pnictides, MnAs, CrAs, and Co As show the structural phase transition be- 
tween the NiAs-type structure above 7^ and the MnP-type crystal structure below 7^ 
with Tt = 398 K, 1,100 K, 1,250 K, respectively. Only MnAs, with decreasing tem- 
perature from Tt = 398 K, shows the second transformation from the MnP-type to 
the NiAs-type structures at 7^ = 3 1 8 K, where magnetic phase transition to the fer- 
romagnetic order also occurs discontinuously (see also Fig. 3.17, etc.). In contrast to 
MnAs, etc., NiAs does not show the transformation to the MnP-type structure, how- 
ever the de Haas-van Alfen effect and the Fermi surface suggest an existence of a 
structural transformation at very low temperatures as noted in Sect. 4.4. 

Various measurements carried out for various mixed compounds revealed that the 
structural transformation between the NiAs-type and the MnP-type structures occurs 
in the following compounds: Cri-j^Mnj^As (all x), Cri-j^Co^As (all x) , Fei-^Co^As 
(0.8 <x< 1.0), Mni_^Ni;cAs (0 < x < 0.58), Mni_;,Co;cAs (0 < x < 0.2, 0.5 < x < 
1.0), MnAsi_xPx (0 < X < 0.28) See also Sect. 3.2. In contrast to these arsenides, all 
the antimonides investigated so far do not show the transformation to the MnP-type 
structure. 

Taking account of the fact that the phase transition occurs in MnAs, CrAs, Co As, 
and VS but not in MnSb, CrSb, NiAs, and TiS, we investigated mechanism of the 
structural phase transition by studying instability of the NiAs-type structure. This is 
a microscopic theory based on band structures obtained by the first principle calcu- 
lations, that is, no fitting parameters are used. 



4.7.1 Electron-Lattice Interaction Coefficient 



Change of energy of an electronic system AF due to lattice deformation is described 
by a phonon normal mode Q^^ (q is a wavevector, A denotes mode) is given by 

^ = -\xiq^)\Qgxf- (4.36) 

X {q^ ) in this equation is a generalized susceptibility, which is given by 

X{q^) = -^I,I,^^j^e-{q^.y)e^\q^y)x''Hvv'.q) (4.37) 



with 



X ^(VV ,q) - ^l^l^lnk^n'k^q \k,n'k+q _ j^O ' (^.38) 



nn k nk n'k+q 
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In these equations, My is mass of the v atom, e{qX^v) is the polarization vector of 
phonon, E^j^ is the band energy in the absence of lattice distortion, and f{E^j^) is the 
Fermi-Dirac distribution function. 

Displacement of a v atom in the a-direction brings about change of crystalline 
potential AV, resulting in coupling between eigen states \nk) and \n'k^q) of the 
system. The coefficient T^^^,^^^ in (4.38) denotes magnitude of the coupling between 
states and ^n'k+q due to the potential change AV , and is called an electron-lattice 
interaction coefficient, which is given by 

i:£,n'k+, = {^nk\m^n'k+,)- (4-39) 

An expression of the electron-lattice interaction coefficient based on the APW for- 
malism is given in [45]. 

4.7.2 Tendency of Structural Transformation from the NiAs-Type 
to the MnP-Type 

The structural transformation to the MnP-type is caused by freezing of a phonon 
mode M4 at the M point of the first Brillouin zone for the Ni As-type structure. The 
normal coordinate of the mode is expressed as 

c\ (zi + Z2) + C2{xi - X2) + C3 (Z3 - Z4) , (4.40) 

where the subscripts 1 and 2 denote metal atoms, and 3 and 4 denote anion atoms. 
The symbol z is the displacement along the c-axis, and x (y) the displacement in the 
c-plane (see Fig. 4.38). As indicated in (4.40), metal atoms shift in the c-plane as well 
as along the c-axis, while anions only displace along the c-axis (see also Chap. 1). 




O metal f J anion 



Fig. 4.38. (a) The NiAs-type structure. The thick lines show the unit cell corresponding to 
the MnP-type structure, (b) Unit cell of the MnP-type structure. The arrows indicate displace- 
ments of atoms which vanish for the NiAs-type structure. We note that magnitude of the 
displacement is exaggerated. The lattice parameters a, b, and c in this figure correspond to A, 
B, C in Fig. 1.4, respectively 



122 4 Electronic Band Structure and Magnetism of Ni As-Type Compounds 



We calculated matrix elements of the electron-lattice interaction I^^^fj^^^ for q = 
FM as a function of k. We note that when a wave vector k is on the X-line (FM), 
A: + ^ is also on the X-line. There are four irreducible representations Xi , X2, X3, and 
X4 in the group of symmetry on the X-line. Among product representations, only 
Xi X X4 and X2 x X3 satisfy compatibility relation with the representation. 

Therefore, matrix elements /'s vanish except for the combination of states Xi- 
and X4 -bands or X2- and X3 -bands. Furthermore, it is clear from (4.38) that the 
electron-lattice interaction reduces the free energy only when one of the states \nk) 
and \n'k') is below the Fermi level and the other is above the Fermi level, and the 
reduction of free energy increases as the energy difference E^j^ — E^,^ decreases. 

MnAs, MnSb, CrAs, and CrSb 

Dispersion curves along the X-line near the Fermi level calculated by using the APW 
method for MnAs, MnSb, CrAs, and CrSb are shown in Fig. 4.39. We see that the 
Xi-band and X4-band lie near the Fermi level, so we calculated I^^^fj^^^ for the states 
k of the Xi-band and A: + ^ of the X4-band. Calculated electron-lattice interaction co- 
efficients are plotted as a function of k along the X-line in Fig. 4.40. The symbols •, 
X , and o denote the results for the displacements in the x-direction of a metal atom, 
in the z-direction of a metal atom, and in z-direction of an anion, respectively. In 
the hatched regions, on the A:-axis, the electron-lattice interaction contributes to the 
reduction of energy of electron system. Figure 4.41 shows the reduction of free en- 
ergy l^l^k^z^k^ql'^/l^^ik - ^^Ak+qV When this quantity is large in the hatched region 
on the ^-axis shown in Fig. 4.41, large reduction of energy of the electron system 
due to the lattice distortions arises, therefore we can say that phase transition to the 
MnP-type structure is likely to occur. The results shown in Fig. 4.41 suggests that the 
phase transition to the MnP-type structure is expected to occur in MnAs and CrAs, 
but not in MnSb and CrSb. 
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Fig. 4.39. Dispersion curves of the nonmagnetic bands of MnAs, MnSb, CrAs, CrSb along 
the X-line. The broken line denotes the Fermi level. The numbers on the curves 1, 2, 3, and 4 
denote the irreducible representations Xi,X2,X3, and 1^4 of the bands (Reproduced form [45] 
with permission) 
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Fig. 4.40. Elecron-lattice matrix elements, |/^^^ E4k-\-q\' MnAs, MnSb, CrAs, CrSb as func- 
tions of the wavevector k. The wavevector q is fixed to FM. The metal x-, metal z-, and anion 
z-components of the electron-lattice matrix elements represented by filled circles, times, and 
open circles, respectively. The hatched regions on the A:-axis denote, where the coupling of 
electronic states due to lattice distortion reduces electronic energy (Reproduced form [45] 
with permission) 
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Fig. 4.41. Energy reduction |/^^^ X4A:+^ 1^/ l^^i^ ~ ^x^k+q I = FM) calculated as functions of 
k. The symbols filled circles, times, and open circles denote the metal x-, metal z-, and anion 
z-components of the MnP-type distortion (Reproduced with permission from [45]) 



We also see from the figure that the energy reduction due to the displacement of 
metal atoms in the x-direction is larger than that due to other two displacements. Ob- 
servation for IMnAs in the IMnP-type structure has revealed that the displacement of 
metal atoms in the x-direction is the largest in the three displacements. 
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From the discussion mentioned earlier, we see that wavevector dependence of the 
electron-lattice interaction coefficients and dispersion curves near the Fermi level 
have a decisive role in the phase transition form the NiAs-type to the MnP-type 
structure. 



CoAs and NiAs 



Using the APW method, we calculated the band structure of CoAs and NiAs near 
the Fermi level for A:-points in the planes = 0, = 1 / SF A, and = 2/ SPA [46] . 
As an example of calculations, we show dispersion curves along the lines parallel 
to the X-line (FM) on the plane of =0 in Fig. 4.42. As we have noted earlier, 
the combination of I^i-band and IJ4-band (or the X2-band and Xs-band) contribute 
to the reduction of electronic energy due to the MnP-type lattice distortion. Only 
when one of the states is above the Fermi level and the other is below the Fermi 
level, the electron-lattice interaction reduces electronic energy. The hatched region 
in Fig. 4.42 indicates where the reduction of electronic energy occurs. 




Fig. 4.42. Dispersion curves of the band of CoAs and NiAs near the Fermi level along the di- 
rection parallel to the FM. the hatched regions on the A:-axis mean where the energy reduction 
due to the electron-lattice interaction occurs (Reproduced with permission from [46]) 
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By calculating |^x|a:X4A:+^I^/I^^i^ —Ej:^k+q\, we found that reduction of elec- 
tronic energy due to the MnP-type lattice distortion is large for CoAs but small for 
NiAs, which indicates that structural transformation to the MnP-type structure is 
likely to occur in CoAs but not in NiAs. This is consistent with observation. Calcu- 
lated electron-lattice interaction coefficients as well as |/^|^ j^^j^j^q\^ /\Ezik — Ez^k^ql 
are depicted in [46] . 

VS and TiSe 

We carried out similar calculations for VS and TiSe. Dispersion curves near the 
Fermi level of these materials show that the electron-lattice interaction between the 
states k of the Xs-band and A: + ^ of the X2-band is important. Calculated electron- 
lattice interaction coefficients show that reduction of free energy due to MnP-type 
lattice distortion is large for VS in comparison with that for TiSe, which suggests that 
VS has a possibility to deform to the MnP-type structure, however TiSe does not. 
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Itinerant Electron Theory of Magnetism 
of Cu2Sb-Type Compounds 



5.1 Crystal Structure and Magnetic Properties 

The Cu2Sb-type compounds that consist of Cr, Mn, Fe, and As or Sb take tetragonal 
crystal structure. Figure 5.1 shows the crystal structure and the first Brillouin zone 
of the Cu2Sb-type structure. The unit cell contains four magnetic atoms denoted by 
(1, 2, 3, 4) and two As or Sb atom denoted by (5, 6). Magnetic atoms occupy the 
sites M(I) and M(II), which have different symmetry. The Cu2Sb-type compounds 
have attracted much interest because they show various magnetic ordering such as 
ferromagnetic, ferrimagnetic, and antiferromagnetic orderings [1]. 

Magnetic structure with many magnetic atoms in a unit cell is described by two 
parameters Q and (^^j. The vector Q (propagation vector) is a wavevector that spec- 
ifies angle between total moment in a unit cell and that in an adjacent unit cell, and 
(j)ij is a relative angle between moments at the ith and the jth sites in a unit cell. The 
structural parameters and magnetic structure (Q, (j)ij, and magnetic moments at the 
respective sites) are shown in Tables 5.1 and 5.2, respectively. 

In CrMnAs and FeMnAs, which are mixed compounds between Mn2As and 
Cr2 As or Fe2 As, a Cr or Fe atom occupies mainly the M(I)-site and a Mn atom occu- 
pies mainly the M(II)-site. In MnAlGe and MnGaGe, which also take the Cu2Sb-type 
structure, Mn atoms occupy only the M(I)-site, and Al and Ga occupy the M(II)- 
site. The compounds MnAlGe and MnGaGe are ferromagnetic below 7c = 503 and 
453 K, respectively. As the Mn atomic layers, as seen in Fig. 5.1, are separated by 
the layers of Al (or Ga) and Ge, these compounds are of interest as a quasi two- 
dimensional ferromagnet. 

5.2 Band Structures of Cu2Sb-Type Compounds and Magnetic 
and Optical Properties 

Since early 1990s, Motizuki and her co-workers have carried out band calculations of 
the Cu2Sb-type compounds for magnetic and nonmagnetic phases to investigate var- 
ious magnetic properties as shown in Table 5.2 [2]. The self-consistent APW method 
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Fig. 5.1. (a) Crystal structure and (b) the first Brillouin zone of the Cu2Sb-type structure 



Tabl e 5.1. Lattice parameters of Cu2Sb-type compo unds 
Lattice parameters 





a (A) 


c(A) 


Uc 




Cr2As 


3.60 


6.34 


0.325 


0.275 


CrMnAs 


3.88 


6.28 


0.327 


0.266 


Mn2As 


3.78 


6.25 


0.33 


0.265 


MnFeAs 


3.735 


6.035 


0.33 


0.25 


Fe2As 


3.627 


5.973 


0.33 


0.265 


Mn2Sb 


4.078 


6.557 


0.295 


0.280 


Cu2Sb 


3.99 


6.09 


0.27 


0.30 


MnAlGe 


3.914 


5.933 


0.273 


0.280 


MnGaGe 


3.963 


5.895 


0.29 


0.29 



Data are taken from [1] 



and the LAPW method are employed [3] together with the muffin- tin approximation 
for crystal potential and the local density approximation for the exchange and corre- 
lation interaction. The formula of Gunnarsson and Lundqvist (GL) was also used for 
the LDA term [4]. As for relativistic effects, we neglected the spin-orbit interaction 
but scalar relativistic effects are taken into account. Densities of states calculated by 
the APW and LAPW methods for the nonmagnetic bands of the Cu2Sb-type com- 
pounds were found to agree well with each other. APW band calculations including 
the spin-orbit interaction were carried out by Kulatov and his co-workers for the 
Cu2Sb-type compounds. They also calculated optical conductivity from the band 
structure [5]. 
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Table 5,2, Magnetic structure of Cu2Sb-type compounds 



Magnetic 
moments 

Magnetic structure (/is /atom) 







Q 


012 


023 


034 


M(I) 


M(II) 


Cr2As 


AF 


(0,0,;r/c) 


n 


K 


0 


0.40 


1.34 


CrMnAs 


AF 


(0,0,0) 


n 


K 


K 


0.41 


3.14 


Mn2As 


AF 


(0,0,:7r/c) 


0 


K 


K 


2.2 


4.1 


MnFeAs 


AF 


(0,0,;r/c) 


0 


0 


K 


0.2 


3.6 


Fe2As 


AF 


(0,0,;r/c) 


0 


0 


K 


1.28 


2.05 


Mn2Sb 


FI 


(0,0,0) 


0 


n 


0 


2.13 


3.87 


Cu2Sb 


P 














MnAlGe 


F 


(0,0,0) 


0 


0 


0 


1.70 




MnGaGe 


F 


(0,0,0) 


0 


0 


0 


1.66 





Data are taken from [1] 



Cr2As Non- magnetic state 




r M X rZ A R ZMAXR 



Fig, 5,2, Dispersion curves of nonmagnetic band of Cr2As 

5.2.1 Nonmagnetic State of Cr2As, Mn2 As, Fe2As, Mn2Sb, CrMnAs, 
and FeMnAs 

We calculated dispersion curves and densities of states for these compounds in non- 
magnetic state [6]. As an example, we show the dispersion curves of Cr2As in 
Fig. 5.2. Two bands in the low energy region consist of As-4^ orbitals. Above the gap, 
there are mixing bands between Cr-3 J orbitals and As-4/? orbitals. Total width of the 
p-d mixing bands is about 3-4 eV, which is large enough to regard the d electrons 
as itinerant electrons. Figure 5.3 shows the total and the partial densities of states 
of M2AS (M = Cr, Mn, and Fe) and Mn2Sb. When we divide the p-d mixing band 
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Fig. 5.3. Densities of states of nonmagnetic (a) Cr2As, (b) Mn2As, (c) Fe2As and (d) Mn2Sb 



shown in Fig. 5.3 into three parts specified by (1), (2), and (3) from low energy side, 
it was revealed by our calculations that (1) and (3) have bonding and antibonding 
characters between M and As (or Sb), respectively, and (2) is regarded as nonbond- 
ing bands consisting almost of d obitals. In the nonbonding bands, J-orbitals at the 
M(I)-sites and M(II)-sites show strong hybridization as seen in Fig. 5.3. 

It is generally said from Fig. 5.3 that the J-band arising from M(II)-sites is nar- 
rower than that from M(I)-sites, which indicates that d electrons at the M(II)-sites 
have weak itinerancy compared with d electrons at the M(I)-sites. Therefore, we ex- 
pect that magnitude of magnetic moment at the M(II)-site is larger than that at the 
M(I)-site. Both calculations and experiments actually indicate that magnetic moment 
of M(II) is larger than that of M(I) as shown in Table 5.3. 
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Ta ble 5.3. Magnetic moments of Cu2Sb-type compounds evaluated by calculati ons 

Magnetic moments (/is) 





Total 


M(I) 


M(II) 


Anion 


Cr2As 






0.34 (0.40) 


1.37 


(1.34) 


0.03 


Mn2As 






1.72 (2.2) 


3.50 


(4.1) 


0.04 


Fe2As 






1.01 (1.28) 


2.11 


(2.05) 


0.03 


Mn2Sb 


1.76 


(1.74) 


-2.11 (-2.13) 


3.67 


(3.87) 


0.01 


MnAlGe 


1.81 


(1.70) 


1.90 


-0.02 




-0.06 


MnGaGe 


2.22 


(1.66) 


2.31 


-0.02 




-0.07 



Observed moments are shown in the parenthesis 
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Fig. 5.4. Density of states of nonmagnetic Cu2Sb 



Fermi level lies in the nonbonding band. As seen in Fig. 5.3, Mn2Sb has larger 
P(£^f), the density of states at Ey^, while compounds other than Mn2Sb have smaller 
P(£^f) than that of Mn2Sb, which is consistent with the fact that Mn2Sb is ferrimag- 
netic and others are antiferromagnetic. 

Next, we show in Fig. 5.4 the density of states of Cu2Sb which shows no magnetic 
ordering [7] . We see that hybridization between J-orbitals of M(I) and M(II) is weak, 
and density of states at Fermi level, p(£^f), is very small. 
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Fig. 5.5. Density of states of ferrimagnetic Mn2Sb 



5.2.2 Ferrimagnetic Band of Mn2Sb 

We show calculated density of states of ferrimagnetic Mn2Sb in Fig. 5.5 [8]. Com- 
ponents of SJ-orbitals of Mn(I) and Mn(II) and 5/?-orbitals of Sb inside muffin-tin 
spheres are also plotted. It is noted that the shape of the densities of states for the 
up-spin and down-spin bands are quite different from each other in the respective 
components, and also different from the densities of states in the nonmagnetic states. 
The band component of Mn(II) shows a large splitting between up- and down- spin 
bands, while the splitting of the band of Mn(I) is small, which results in that the mag- 
netic moment of Mn(II) is about two times larger than that of Mn(I) (see Table 5.3). 
We also see that small magnetic moment is induced at the Sb-site in the direction of 
moment of Mn(II). Calculated total magnetic moment as well as those of Mn(I) and 
Mn(II) agrees well with experimental vales as seen in Table 5.3. 

The ferrimagnetic band of Mn2Sb is also calculated by Haas and his co-workers 
using the ASW method [9, 10]. Kulatov also carried out band calculation using the 
LMTO method. Gross features of density of states obtained by them are similar to 
the result shown in Fig. 5.5. 
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5.2.3 Ferromagnetic State of MnAlGe and MnGaGe 

In these compounds, atomic planes of Mn at M(I)-site are separated by two planes 
consisting of Al (or Ga) and Ge, which suggests a two-dimensional character in mag- 
netic properties. Observed magnetic moment is much smaller than moment expected 
for localized Mn orbitals. Therefore, we can regard to some extent these compounds 
as an itinerant electron ferromagnets with two-dimensional character. 

Figure 5.6 shows calculated density of states for MnAlGe [11]. MnAlGe has 
density of states similar to that of MnGaGe. Hybridization between J-orbitals of Mn 
and /?-orbitals of Al (Ga) and Ge is small. As shown in Table 5.3, calculated moment 
of MnAlGe is in good agreement with the experimental value, while the agreement 
is not good for MnGaGe. According to the band calculations, very small magnetic 
moment is induced at Al (Ga) and Ge atomic site in the opposite direction of the 
Mn-moment. 

We found, as shown in Fig. 5.7, three Fermi surfaces of the down-spin bands 
consisting of two electron surfaces and a cylindrical hole surface around the Z-axis, 
and four Fermi surfaces of the up- spin bands consisting of a small hole pocket around 
the R-point, a big hole pocket, and two electron pockets near the M-point. The Fermi 
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Fig. 5.6. Density of states of ferromagnetic bands of MnAlGe 
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Fig. 5.7. Fermi surfaces of ferromagnetic bands of MnAlGe 



surfaces of the down-spin bands has large contribution to p(£'f). We expect that the 
strong two-dimensional character of the Fermi surfaces of the down- spin bands will 
be observed as anisotropy of transport properties such as resistivity. 

5.2.4 Antiferromagnetic Bands of Cr2 As, Mn2 As, and Fe2 As 

We carried out band calculations for antiferromagnetic state of these compounds. As 
seen in Fig. 2.2 of Part I and Fig. 5. 1 of Part II, magnetic moments of Cr(l) and Cr(2) 
of Cr2As are in the c-plane and antiparallel, and so Cr(l) and Cr(2) are regarded 
as different sites. On the other hand, in Mn2As and Fe2As, the sites 1 and 2 are 
equivalent because moment at these sites are parallel. Magnetic moments at the M(I)- 
site and adjacent M(II)-site are antiparallel in Mn2As and parallel in Fe2As. The spin 
arrangements for the compounds with Cu2Sb-type structure are illustrated in Fig. 2.2 
of Parti. 
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Fig. 5.8. Total and local densities of states for (a) Cr2As, (b) Mn2As, and (c) Fe2As in anti- 
ferromagnetic state 



Figure 5.8 shows total and local densities of states for MAs with M = Cr, Mn, and 
Fe. It is noteworthy that in Cr2As the splitting of Cr(II)-bands is about 1 eV, whereas 
there is little splitting of Cr(I)-bands. The very small splitting in Cr(I)-local density 
of states results in the very small magnetic moment of Cr(I), which is related to the 
antiferromagnetic order in the c-plane. As for Fe2As and Mn2As, mixing between 
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3J-states at the M(I)- and the M(II)-sites in Mn2As is smaller than that in Fe2As. 
Calculated moments are compared with the experimental data in Table 5.3. As shown 
in Table 5.3, magnetic moment of M(II) is larger than that of M(I) in the respective 
compounds. Agreement between calculation and experiment is good for Cr2As and 
Fe2As but not so good for Mn2As. 



5.2.5 Magnetic Ordering of Cu2Sb-Type Compounds 

The Cu2Sb-type compounds except for MnAlGe and MnGaGe have four magnetic 
atoms in a unit cell, and therefore various types of magnetic structures have been 
observed. We have studied such various magnetic structures from the viewpoint of 
itinerant electrons. For this picture we investigated instability of paramagnetic phase 
within the theoretical framework described in Sect. 4.6. In this section we explain 
only the calculated results. 

Figure 5.9 shows the magnetic phase diagram calculated for Q = (0,0,0) and 
Q= (0,0, tt/c), respectively. The dashed-curves express the relation between the 
critical values of intra-atomic Coulomb interaction divided by a transfer integral. 
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Fig. 5.9. A magnetic phase diagram of the Cu2Sb-type compounds. The symbols open circle, 
open triangle, open square, etc. denote critical values of intra-atomic Coulomb interaction 
calculated for each magnetic ordering (see also the text). The magnetic orders expressed by 
CrMnAs(II)-type and Cr2 As-type are hypothetical ones for the theoretical calculation [14] 
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f/c/|^i|. f/c/|^i| is plotted as a function of average electron number n. As we have 
noted in Sect. 4.6, if U is larger than paramagnetic phase is unstable. The magnetic 
structure that corresponds to the smallest value of Uq is most likely to realize. 

From the calculations, it is indicated that Mn2Sb-type, CrMnAs(II)-type, and fer- 
romagnetic ordering are theoretically possible for Q = (0,0,0). For 2 = (0,0, tt/c), 
the Cr2As(II)-type, Mn2 As-type, and Fe2 As-type can theoretically realize. The mag- 
netic orderings actually observed are included in Fig. 5.9. 



5.2.6 Photoemission and Inverse Photoemission 

It is possible to obtain a rough estimation of density of states below and above the 
Fermi level by photoemission and inverse photoemission measurements, respec- 
tively. Suga and Kimura measured photoemission and inverse photoemission spec- 
tra of Cu2As, Mn2As, Fe2As, MnAlGe, and MnGaGe. Their experimental data are 
shown in Fig. 5.10 [12, 13] together with the calculated density of states summed 
over the up-spin and down-spin bands [14]. Overall features of the spectra corre- 
spond well to the calculations. For Mn2Sb, the spectrum has remarkable peaks at 
about 3 eV below the Fermi level and 2 eV above the Fermi level, while in MnAlGe 




Fig. 5.10. Densities of states determined by photoemission and inverse photoemission spectra 
(dots) and calculated densities of states (solid lines) for the Cu2Sb-type compounds 
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there is no clear peak but broad hump with the width about 4.5 eV lies below Fermi 
level. Such difference between the spectra of Mn2Sb and MnAlGe can be attributed 
to the absence of Mn at M(II) site in the case of MnAlGe. 

From the spectra of Cr2As, Mn2As, and Fe2As, we see that energy difference 
between the peaks above and below the Fermi level in Cr2As and Fe2As is smaller 
than the energy difference in Mn2As. 
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-of Fe2As, 130 

- of ferrimagnetic Mn2Sb, 132 
-of Mn2As, 130 
-ofMn2Sb, 130 

- of antiferromagnetic Cr2As, 134 
-of CoAs, 83 

-ofCrSb, 87, 88 
-ofCrX, 86 
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-ofMnAs, 76-78 

- of MnP-type CoAs, 84 
-of MnP-type FeAs, 84 
-ofMnSb, 76, 79 

-of NiAs, 105 
-ofNiSb, 105 

- of antiferromagnetic Fe2As, 134 

- of antiferromagnetic Mn2As, 134 

- of ferromagnetic CrTe, 90 

- of ferromagnetic MnAs, 91 

- of ferromagnetic MnAlGe, 133 

- of ferromagnetic MnAs, 79-81 

- of ferromagnetic MnSb, 91 
-ofCrAs, 87 

-of CrP,56 
de Haas-van Alfen effect, 103-106 
dielectric function, 92 
dipole approximation, 92 
dispersion curves 

-of CoAs, 124 

-of Cr2As, 129 

-of CrAs, 122 

-of CrSb, 122 

-of MnAs, 77, 122 

-ofMnSb, 122 

-of NiAs, 104, 124 

-ofNiSb, 104 
displacement parameter, 7 
distortion parameter, 7 
double helical ordering, 12, 13, 17, 26, 29, 

48, 53, 54, 83, 85, 113 
Drude term, 92 

elastic coefficient, 31, 32, 35 
elastic energy, 21, 31, 32, 49, 52, 59, 60 
electrolytic manganese, 4 
electron-lattice interaction, 1 20- 1 25 
electronic pressure, 100 
entropy, see magnetic entropy, lattice 
entropy 

exchange and correlation interaction, 75, 77, 

84, 85, 91, 104, 128 
exchange energy, 21 
exchange integral, 95 

exchange interaction, 20, 24, 25, 35, 55, 58, 
63 

exchange inversion model, 56-63 



Fermi level, 77, 79, 83-88, 92, 99, 102, 104, 

105, 122, 124, 131, 137, 138 
Fermi surface, 120 

-of CoAs, 83 

-of FeAs, 116 

-of NiAs, 83, 106 

-OfNiSb, 108 

- of ferromagnetic MnAlGe, 134 
Fermi-Dirac distribution function, 92, 100, 

118, 121 
first Brillouin zone 

- of Cu2Sb-type structure, 128 

- of NiAs-type structure, 76 

FLAPW method, 75, 77, 80, 81, 89, 90, 104, 
109, 111 

formula of Gunnarsson and Lundqvist, 75, 

77, 84, 104, 128 
formula of von Barth and Hedin, 75, 77, 

91 

FP-LMTO method, 85 
free energy, 20, 60, 61, 69-71, 96, 100, 122, 
125 

functional integral, 96 

generalized gradient approximation, 76, 84, 
103 

generalized susceptibility, 120 
Green's function, 101 
Green's function method, 77 
Gruneisen relation, 49 

Hartree-Fock approximation, 97, 101, 114 
hexagonal close packed structure, 3, 5 
high spin state, 29, 30 
Hubbard Hamiltonian, 95, 114 
hysteresis, 20, 22, 25, 37, 38, 45, 61, 62, 64, 
65,94 

instability 

- of NiAs-type structure, 120 

- of paramagnetic state, 112-114, 118, 

136 

inverse photoemission, 137 
isothermal process, 43, 44 
itinerant electron, 20, 26, 32, 36, 76, 95, 100, 
112, 129, 136 

Kittel's model, 56-63 
Kubo formula, 92 
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latent heat, 47 
lattice constant 

- of Cr-compounds, 49 
-ofCrAsi_;,P;„ 50 

- of CrX compounds, 48 

- of Cu2Sb-type compounds, 128 
-ofMn2-;cCr;cSb, 57 
-ofMnAs, 19 

-of MnAso.vSbo.s, 34 
-of MnAsi_;cPx, 32 
-ofMnAsi_;,Sb;c, 34 

- of MnX compounds, 6, 27 
lattice entropy, 21, 42, 43 
linear response theory, 92 

LMTO method, 75-78, 80, 84, 91, 92, 132 
local density approximation, 75, 76, 83, 85, 

103, 104, 128 
local density spin approximation, 75 
local magnetic moment, 97, 99, 101, 102 
local moment, 20, 25, 26, 33, 37, 38, 95 
low spin state, 29, 30 

magnetic energy, 49, 52, 53 

magnetic entropy, 21, 41, 42, 44, 52, 60, 70 

magnetic moment 

-of Cr2As, 134 

-ofCrAs, 50, 85 

-ofCrSb, 55, 85, 88 

-of CrTe, 89, 107, 111 

- of CrX compounds, 12, 48 

- of Cu2Sb-type compounds, 131 
-of Fe2As, 134 

-ofFeAs, 113, 119 

- of FeX compounds, 13 

- of M2X compounds, 13 
-ofMn2As, 134 
-ofMn2Sb, 132 

-of Mni_;,M;,AlGe, 65 
-of Mni_;cTixAs, 41 
-of Mn2_xCr;cSb, 62 
-ofMnAlGe, 133 
-ofMnAs, 18, 78, 80,81,94 

- of MnAso.7Sbo.3, 36, 39, 40 

- of MnAso.88Sbo.12, 38 
-of MnAso.gPo. 1,28-30 
-of MnAsi_;,Px, 28 
-of MnBi, 45 
-ofMnGaGe, 133 

-of MnP, 17 



-ofMnSb, 9, 79 

- of MnX compounds, 12, 27 
magnetic refrigeration, 40-45 
magnetic structure 

-of CrSb, 56 

- of CrX compounds, 12 

-of Cu2Sb-type compounds, 13, 14, 127, 
129 

- of FeX compounds, 13 

- of M2X compounds, 13 

- of MnX compounds, 12 
magnetization process, 94 
magnetocaloric effect, 42, 44 
magnetostriction, 20, 31, 49, 53 

- of MnAso.93Po.07, 31 
magnetovolume effect, 20, 21, 24, 31, 48, 49 
matrix element, 92, 121 
metamagnetism, 11, 19, 20, 22, 23, 29-32, 

36, 39,41,94 
MnP-type structure, 3, 5, 121 
mode-mode coupling, 96, 102 
model density of states, 97, 98, 102, 103 
molecular field, 52 

muffin-tin approximation, 75-77, 83, 91, 
128 

muffin-tin potential, 85 
muffin-tin sphere, 77, 79, 132 

Neel temperature 

- of Cr-compounds, 51 
-of CrAs, 46, 50 
-of CrAso.gPo.i, 46 
-of CrAsi_;,Sb;c, 51 
-of CrAsi_;,Sb;c, 54 

- of CrX compounds, 12 

- of FeX compounds, 13 

- of M2X compounds, 13 
nesting, 83, 84, 106 

neutron diffraction, 47, 53-55, 78, 94, 95 
NiAs-type structure, 3, 5, 76, 121 
non-bonding band, 130 
normal coordinate, 120, 121 

optical absorption coefficient, 92 
optical photoconductivity, 92 
orthorhombic symmetry, 5 

partition function, 96 

Pauli paramagnetism, 40, 50, 102 
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peritectic temperature, 45 
phase transition 

first order - 11, 17, 20, 25, 26, 33, 36, 37, 
39, 41, 44-48, 50-54, 56, 58, 60, 63, 
64, 93, 99 

pressure-induced-, 107, 111 

second order -, 33, 37, 39, 99 
phonon, 120, 121 
phosphorus, 4 
photoemission, 137 
polarization vector, 121 

red phosphorus, 4 
refrigerant material, 40 
relativistic effect, 75, 104, 128 
resistivity, 47 

-of CrSb, 55 
RKKY interaction, 63 

self-energy, 101 
shear, 59 

simple hexagonal lattice, 5 
sintering method, 3 
specific heat, 42, 47, 49, 79, 88, 105 
spin density wave, 117, 118 
spin fluctuation, 83, 84, 93-103 
spin glass, 40 

spin-orbit interaction, 75, 91, 92, 104, 128 
spin-orbit interaction, 128 
stiffness constant, 59 



Stratonovich-Hubbard transformation, 95 
structural transformation, 81, 84, 89, 94, 95, 

99, 102, 106, 119-125 
super-exchange interaction, 63 
susceptibility 

-of CoAs, 102 

-of CrAs, 46 

-of CrAso.9Po.i, 46 

- of CrAso.gSbo.i, 53 

-of CrAsi_;,Sb;c, 51 

-ofCrP,56, 88 

-of CrSb, 55 

-ofFeAs, 102 

-ofMnAs, 18, 94, 98 

-ofMnAso.9Po.i,28 

-ofMnBi, 45 

-of MnP, 17 

-of NiAs, 102 

thermal expansion, 20 

thermal expansion coefficient, 18, 31, 32, 49 
3d electron, 20 
transfer integral, 95, 114 

white phosphorus, 4 

X-ray diffraction, 53, 63 
XPS, 79, 92, 94 

Zeeman energy, 21, 69 



